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Author: Bunnjaweht, Dahmmaet 
Degree: Ph.D. 
Year: 2005 

Corporate Source/Institution: University of Colorado at Colorado Springs ( 0892 ) 
Source: Volume 6512B of Dissertations Abstracts International. 
PAGE 6543 . 163 PAGES 
ISBN: 0-496-92246-7 

...quality of service, the system requires efficient resource allocation. In the system based on a code 
division multiple access technique, every user contributes to the interference characteristics of the system. 
This interference power is proportional to the transmitter power. It is clear that power control is the most 
important factor which allows more efficient resource sharing by balancing received signal quality and 
interference power. 

The multiple path connections idea is well-studied in the network layer of a wired infrastructure. Here 
we investigate a similar idea in the physical layer of a wireless access network. There is a major difference 
in the sense of channel existence. In the wired network, where the channels always exist, algorithms are 
based on smart queues and routings. In the connections scheme for resource sharing. 

The proposed scheme is studied under two different types of channels based on wideband code division 
multiple access standard, the dedicated channel and the common packet channel. The dedicated channel 
represents real-time traffic, while the common packet channel represents non-real-time traffic. Wideband 
code division multiple access technology is chosen because it has emerged as the most widely adopted third 
generation interface.. 
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Reduction of NTSC co-channel interference by referencing carrier frequencies to the LORAN-C 
signal 

Rhodes, C W 

IEEE Transactions on Broadcasting , v 41 , n 2 , p 37-43 , June 1995 
Publication Date: 1995 

Publisher: Institute of Electrical and Electronics Engineers, Inc. , 445 Hoes Ln , Piscataway , NJ , 08854- 
1331 

Country Of Publication: USA 
Publisher Url: http://ieee.org 
Publisher Email: inspec@ieee.org 

Document Type: Journal Article 
Record Type: Abstract 
Language: English 
ISSN: 0018-9316 

File Segment: Electronics & Communications Abstracts 

Reduction of NTSC co-channel interference by referencing carrier frequencies to the LORAN-C 
signal 



Abstract: 

Co-channel interference may be noticed when the desired signal at the receiver is 60 dB stronger than the 
undesired signal. The spacing distances between transmitters which would be required to avoid this 
interference are prohibitive. Because FCC rules require much smaller minimum spacings, frequency offset 
is specified for each station to reduce co-channel interference. Hie offsets used in the U.S. are +/-10 kHz 
while the carrier frequency tolerance is +/-1 KHz. Differences between carriers can be 8-12 kHz. 
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Interference minima occur every 59.94 Hz within this range with interference maxima in between these 
minima. By referencing the visual carrier frequency of each transmitter to a common frequency, the 
difference frequency can be maintained at one of those frequencies at which the interference is minimum. 
This paper explains how the LORAN-C signal can be used as the frequency reference for television 
transmitters and how the ideal difference frequency can be maintained to minimize the visibility of co- 
channel interference. Equipment required is commercially available, and no rule changes are required to 
employ LORAN-C referenced carrier operation 

Descriptors: Interference; Loran-C; Transmitters; Carriers; Offsets; Carrier frequencies; Minima; 
Referencing; Visibility; Receivers; Maxima; Television; Reduction; Stations; Tolerances; Broadcasting; 
Visual; Cobalt; Face centered cubic lattice; Noise levels 
Identifiers: 
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Title: Three-dimensional imaging using ambiguity-free interferometry 
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Author Affiliation: Inst, fur Hochstfrequenztech. und Elektron., Karlsruhe Univ., Germany 
Conference Title: AMTA '97. 19th Meeting and Symposium p. 445-9 
Publisher: Antenna Meas. Tech. Assoc , Boston, MA, USA 
Publication Date: 1997 Country of Publication: USA xxii+574pp. 

Material Identity Number: XX- 1999-00546 
Conference Title: AMTA '97. 19th Meeting and Symposium 
Conference Date: 17-21 Nov. 1997 Conference Location: Boston, MA, USA 
Language: English 
Subfile: B 

Copyright 1999, IEE 
Title: Three-dimensional imaging using ambiguity-free interferometry 

Abstract: One of the key problems in S AR interferometry is the determination of the absolute phase of a 
scatterer by unwrapping the calculated phase difference of two S AR images. Since the phase difference is 
a function modulo 2 pi , software algorithms are used to perform a phase unwrapping to obtain an 
unambiguous phase and thus the corresponding height information. This paper presents a procedure using 
two different transmit frequencies to enlarge the unambiguous range for height determination. 
Measurements performed in an anechoic chamber are used to test the processing. Using this procedure it... 
Descriptors: ...radiowave interferometry; 

Identifiers: ...ambiguity-free interferometry; phase difference; transmit frequencies 
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06773897 INSPEC Abstract Number: A9802-9385-086, B9801-7710-034 
Title: An interferometric SAR processor avoiding phase ambiguities 
Author Schmitt, K.; Wiesbeck, W. 

Author Affiliation: Inst, fur Hochstfrequenztech. und Elektronik, Karlsruhe Univ., Germany 
Conference Title: IGARSS'97. 1997 International Geoscience and Remote Sensing Symposium. Remote 
Sensing - A Scientific Vision for Sustainable Development (Cat. No.97CH36042) Part vol.4 p. 1713- 
15 vol.4 

Editor(s): Stein, T.I. 

Publisher: IEEE , New York, NY, USA 

Publication Date: 1997 Country of Publication: USA 4 vol. lxxi+2105 pp. 
ISBN: 0 7803 3836 7 Material Identity Number: XX97-02132 
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U.S. Copyright Clearance Center Code: 0 7803 3836 7/97/$10.00 

Conference Title: IGARSS'97. 1997 IEEE International Geoscience and Remote Sensing Symposium 
Proceedings. Remote Sensing - A Scientific Vision for Sustainable Development 
Conference Sponsor: IEEE Geosci. & Remote Sensing Soc; Centre for Remote Imaging, Sensing & 
Processing, Nat. Univ. Singapore; NASA; NOAA; Office of Naval Res.; URSI 
Conference Date: 3-8 Aug. 1997 Conference Location: Singapore 
Language: English 
Subfile: A B 
Copyright 1997, IEE 
Title: An interferometric SAR processor avoiding phase ambiguities 

Abstract: One of the key problems in SAR interferometry is the determination of the absolute phase of a 
scatterer by unwrapping the calculated phase difference of two SAR images. Since the phase difference is 
a function modulo 2 pi software algorithm are used to perform a phase unwrapping to obtain an 
unambiguous phase and thus the corresponding height information. This paper presents a procedure using 
two different transmit frequencies to enlarge the unambiguous range for height determination. 
Measurements performed in an anechoic chamber are used to test the processing. Using this procedure it... 

Identifiers: ...interferometric SAR processor SAR interferometry; phase difference; 
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Apparatus and method for reducing multiplicative gain variation distortions in data recording and 
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, USA 
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Document Type: Patent 
Record Type: Abstract 
Language: English 

File Segment: ANTE: Abstracts in New Technologies and Engineering 

Apparatus and method for reducing multiplicative gain variation distortions in data recording and 
transmission channels 



Abstract: 

Data channel gain- variation distortions are reduced by the following combination of operations. 
Successively occurring data items are cumulatively added in modulo fashion thereby providing a sequence 
of partial sums. An exponentiating circuit is introduced to provide an exponential function of the partial 
sums. The exponentiated partial sums are then recorded or transmitted on a data channel wherein the 
signals are distorted by relatively slow gain variations thereof. The received distorted signals are applied to 
a logarithmic circuit wherein the multiplicative gain variation distortions are converted to additive 
distortions. The logarithmically converted signals are applied to a subtracter that sequentially provides 
the differences in modulo fashion between successively occurring signals from the logarithmic circuit 
thereby yielding the data items with the distortions substantially reduced. 

Descriptors: Distortion; Circuits; Sums; Channels; Additives; Exponential functions; Data recording 
Identifiers: 
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(c) 2009 INIST/CNRS. All rights reserved. 
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Modeling and experimental verification of the interconnected mesh power system (IMPS) MCM topology 
LOW Y L; SCHAPER L W; ANG S S 

Bell Laboratories, Lucent Technologies, Murray Hill, NJ 07974, United States; High Density 
Electronics Center (HiDEC) and Department of Electrical Engineering, University of Arkansas, 
Fayetteville, AR 72701, United States 

Journal: IEEE transactions on components, packaging, and manufacturing technology. 
Part B, Advanced packaging, 1997, 20 (1 
) 42-49 

Language: English 

Copyright (c) 1997 INIST-CNRS. All rights reserved. 



...by half. To provide a complete MCM on only two metal layers, the IMPS topology provides a 
unique and complex power distribution system and signal transmission environment. This paper 
reports on models of the IMPS structure based on finite-difference/partial-inductance methods 
and PS pice simulation. The models can easily be constructed using a commercial three-dimensional 
(3-D) field solver and the circuit models... 

English Descriptors: Network analysis; Integrated circuit; Multichip module; Modeling; 
Interconnection; Transmission line; Noise; Characteristic impedance; Meshed network; 
Finite difference method; Distribution system 

French Descriptors: Analyse circuit; Circuit integre; Module multipuce; Modelisation; Interconnexion 
; Ligne transmission; Bruit; Impedance caracteristique; Reseau maille; Methode 
difference finie; Systeme distribution 

Spanish Descriptors: Analisis circuito; Circuito integrado; Modulo multipulga; Modelizacion; 
Interconexion; Linea transmision; Ruido; Impedancia caracteristica; Red mallada cerrada; Metodo 
diferencia finita; Sistema distribucion 
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Title: Blind ISI mitigation 

Author Mostofi, Y.; Cox, D.C. 

Author Affiliation: Stanford Univ., CA, USA 

Conference Title: IEEE 54th Vehicular Technology Conference. VTC Fall 2001. Proceedings (Cat. 
No.01CH3721 1) Part vol.3 p. 1903-7 vol.3 
Publisher: IEEE , Piscataway, NJ, USA 

Publication Date: 2001 Country of Publication: USA 4 vol.(lxxiii+xii+2777) pp. 
ISBN: 0 7803 7005 8 Material Identity Number: XX-2001 -02309 
U.S. Copyright Clearance Center Code: 0-7803-7005-8/01/$10.00 

Conference Title: IEEE 54th Vehicular Technology Conference. VTC Fall 2001. Proceedings 
Conference Date: 7-11 Oct. 2001 Conference Location: Atlantic City, NJ, USA 
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Language: English 
Subfile: B 

Copyright 2002, IEE 

Abstract: In an outdoor communication environment, a percentage of bandwidth is wasted sending a 
training sequence for channel estimation and equalization. For instance in GSM, 17.93% of bandwidth is 
dedicated to the transmission of such a known sequence. Therefore, if a blind algorithm, using only the 
knowledge of input constellation and/or correlation, can achieve an acceptable performance... 

Descriptors: ...adaptive estimation; interference suppression intersymbol interference; ... 

...parameter estimation; 

Identifiers: ...channel estimation; channel equalization intersymbol interference; 
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The COMPASS RICH-1 read-out system 

Baum, G.; Birsa, R.; Bradamante, F.; Bressan, A.; Chapiro, A.; Cicuttin, A.; Ciliberti, P.; Colavita, A.; 
Costa, S.; Crespo, M.; Cristaudo, P.; Dalla Torre, S.; Diaz, V.; Fauland, P.; Fratnik, F.; Giorgi, M.; Gobbo, 
B.; Ijaduola, R.; Kalinnikov, V.; Lamanna, M.; Martin, A.; Menon, G.; Pagano, P.; Schiavon, P.; 
Tessarotto, F.; Zanetti, A. 

Corresp. Author/Affil: Crespo, M.: INFN, Sezione di Trieste, ICTP, Via Beirut 31, 34100 Trieste, Italy 
Corresp. Author email: mcrespo@ictp.trieste.it 
Editor(s): Ekelof, T; Resvanis, L.K.; Sequinot, J 
Conference Title: RICH 2002 

Conference Location: Nestor, Pylos Greece Conference Date: 20020605-20020610 
E.I. Conference No.: 63024 

Nuclear Instruments and Methods in Physics Research, Section A: Accelerators, Spectrometers, Detectors 
and Associated Equipment ( Nucl Instrum Methods Phys Res Sect A ) ( Netherlands ) 2003 502/1 (246- 
250) 

Publication Date: 2003041 1 
Publisher: Elsevier 

CODEN: NIMAE ISSN: 0168-9002 
Publisher Item Identifier: SO 168900203002821 
Item Identifier (DOI): i0.1016/SQ168-9002(03)00282-l 

Document Type: Conference Paper; Journal Record Type: Abstract 
Treatment: T; (Theoretical) 

Language: English Summary Language: English 
Number of References: 18 

This paper describes the reconfigurable read-out system for the 82944 RICH-1 channels of the 
COMPASS experiment (NA58) at CERN. The system is based on 192 identical large front-end boards 
(BORA board). BORA was designed for acquiring, digitizing, threshold subtracting and transmitting 
event data. The overall operation of the board is controlled and supervised by a DSP tightly interacting with 
an FPGA that acts as a parallel co-processor. The DSP allows characterizing each analog channel by 
locally calculating noise and pedestal. Each BORA communicates with the outside world through two 
optical fibers and through a dedicated DSP network. One optical fiber is used to receive event triggers, and 
the other one is used to transmit event data to subsequent processing stages of the acquisition system. The 
DSP network allows reconfiguring and reprogramming the DSPs and FPGAs as well as acquiring... 
...DOLINA DSP. The read-out system can steadily work up to a trigger rate of 75 kHz with maximum pixel 
occupancy of 20%, reaching a transmission data rate of 5.13 Gbytes/s. (c) 2003 Elsevier Science B.V. All 
rights reserved. 

Descriptors: 
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Bidirectional access on HFC CATV networks: Design of a CDMA-based communications system 

Author: Wolters, Robert Peter Christina 
Degree: Ph.D. 
Year: 1999 

Corporate Source/Institution: Technische Universiteit Eindhoven (The Netherlands) ( 0426 ) 
Source: Volume 6101C of Dissertations Abstracts International. 
PAGE 302. 214 PAGES 
ISBN: 90-423-0072-8 

Publisher: Shaker Publishing B.V., St. Maartenslaan 26, 6221 AX Maastricht, The Netherlands 

...coaxial access network. In older upgrades, a PON like architecture was applied in the downstream 
direction, requiring that special techniques are applied. The CATV upstream channels are hampered by 
ingress noise and radiation from short-wave radio, and by non-ideal amplitude and phase characteristics. 
The most suitable transmission technique, as concluded in this thesis is Direct Sequence Code Division 
Multiple Access (DS-CDMA). 

The upstream physical design is based on DS-CDMA, employing preferentially phased Gold codes. 
Since the system is operating in synchronous mode except during the initialisation procedures, the 
downstream channel has been defined as the timing reference for all upstream channels. The differences 
in delay resulting from different physical positions of the modems is compensated for by an initialisation 
procedure. 

The Medium Access Control determines which modem is allowed to use a certain spreading code and 
time slot. It is therefore important that a modem is assigned bandwidth when it is required. A distinction 
needs to be made between modems that actively transmitting data, and modems that are connected to the 
network but are not actively sending data. Active modems use piggybacking for sending the bandwidth 
needs to the MAC controller. Idle modems use dedicated cells instead. Based on these requests, the MAC 
algorithm determines for each modem when it is allowed to transmit. The modem determines for which 
traffic type the assigned bandwidth is used. 
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Integrated data processing system with CPU and digital signal processor - integrates functions of 
CPU and DSP modules as parallel, independently operating functions, tuned for voice applications 
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...Original Titles: Integrated circuit having CPU and DSP for executing vector lattice propagation 

instruction and updating values of vector Z in a single instruction cycle Parallel operating CPU core and 

DSP module for executing sequence of vector DSP code instructions to generate decoded constellation 

points in QAM/TCM modem application Integrated circuit having CPU core operable for switching 

between two independent asynchronous clock sources of different frequencies while the CPU continues 
executing instructions. Original Publication Data by Authority ArgentinaPublication No. ...Original 
Abstracts: execution by the DSP module and for access by the GP CPU for storage and retrieval of 
instructions and data. The DSP module includes a lattice filter mechanism for executing a single 
instruction to update propagating values of a vector responsive to the externally provided signal. ... the DSP 
module is independent of and in parallel with execution of GP instructions by the GP CPU. The system 

executes a novel QAMYTCM modem algorithm The present invention is directed to various features 

of an integrated data processing system that includes a general purpose (GP) CPU core for processing data 

in accordance with a GP instruction set and Claims: being independent of and in parallel with execution 

of GP instructions by the GP CPU; and wherein the DSP module executes a dual tone multi-frequency 
(DTMF) detection algorithm by executing a single command-list code instruction VCPOLY) that computes 
ACCEPT frequency detectors, REJECT frequency detectors and total signal energy; and wherein the GP 
CPU includes decision logic that identifies a DTMF signal based on the ACCEPT frequency detectors, 
REJECT frequency detectors and total signal energy computed by the DSP module.An integrated circuit 
structure for use in identifying a parallel disconnect signal provided on an analog transmission medium, 
the integrated circuit structure comprising, a central processing unit that executes instructions to perform 
data processing operations, the data processing operations including a successive approximation analog-to- 
digital conversion operation to provide a digital output signal based upon an input data signal; a pulse 
width modulation (PWM) element that converts the digital output signal generated by the central 

processing unit to a square-wave output signal having a duty cycle circuit structure, the external 

comparator receiving as inputs the output of the external low pass filter and the parallel disconnect signal 
provided by the analog transmission medium and providing said external comparator output a digital 
signal corresponding to a difference between the output of the external low pass filter and the signal 
provided by the analog transmission medium, the input port connected to provide the output of the 

comparator as the input data signal to the central processing unit command-list code instructions and 

vector operand data for utilization by the DSP module, the DSP memory element connected to the DSP 
module by a dedicated DSP bus, but not connected to the shared internal bus; and wherein the DSP module 
includes at least one vector address pointer register that includes an address field that identifies an address 
of vector operand data that can be retrieved from the DSP memory element by the DSP module via the 
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DSP bus for use by the DSP module in executing vector DSP command-list code instructions, an address 
wrap-around field that facilitates wrap... processing means for processing an externally-provided signal 
received by the DSP module by executing DSP command-list code instructions, the signal processing 
means including lattice means for executing a first of said command-list code instructions to update 

propagating values of the vector "Z" responsive to the externally-provided signal said first of said 

command-list code instructions requires only a single instruction fetch/decode cycle of the DSP module to 

update the propagating values of the vector "Z for preventing the DSP module from executing 

command-list code instructions which cause access to the shared external memory, whereby the DSP 
module is temporarily prevented from 
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...uses receiver having data aided multi-user interference canceller with feed forward and feed 
backward filters and subtracting stage ...Original Titles:Direct sequence code division multiple access 

receiver with multiuser interference canceller Recepteur pour acces multiple par multiplexage de 

code a sequence direct avec annuleur d T interference pour utilisation multiples Direct sequence - Code 

division multiple access receiver with multiuser interference canceller Recepteur pour acces multiple 

par multiplexage de code a sequence direct avec annuleur d' interference pour utilisation multiples 
Alerting Abstract ...The system has a transmitter for multiplying a modulating signal with a spreading 
sequence to form a spread spectrum signal and transmitting the signal on a CDMA channel. A receiver 
receives the signals and cancels the multiple access interference. It uses a feedforward filter, feedback 
filter and quantiser for producing quantised decisions from information derived from the feedforward 

filter Circuits make tentative decisions on received information produce an inversion of the equivalent 

CDMA channel in dependence on the cross-correlation values of all users, and adjust the tap weights of the 
feedforward filter... Title Terms .../Index Terms/Additional Words: INTERFERENCE; ... 
...SUBTRACT; Class Codes International Patent Classification IPC Class Level Scope Position Status 
Version Date ...H04K-001/02 Main "Version 7" H04B-0001/707... H04B-0001/707... Original 
Publication Data by Authority ArgentinaPublication No. Original Abstracts: A DS-CDMA system 
comprises at least one transmitter and at least one receiver operating on a dedicated signal channel which 
is used by other DS-CDMA operators in a non- synchronous way. A receiver for use in such a DS-CDMA 
system includes a data aided multiuser interference canceller (38) comprising a feedforward filter (48) 
connected to a subtracting stage (50). A decision stage (52) is connected to an output of the subtracting 
stage (50). Decisions from the decision stage (52) are stored in a memory which also stores preliminary 
decisions made on the signal applied to the feedforward filter (48) by a simple receiver (46). These 
decisions are applied to a feedback filter (58) the output of which is connected to the subtracting stage 
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(50) thereby removing the effects of pre-cursive and post-cursive interference. The input to, and output 
from,the decision stage 52 are compared and the difference is used for adjusting the tap weights of the 
feedforward and feedback filters (48,58). In the case of a single input, single output, the feedforward filter 

(48) comprises A DS-CDMA system comprises at least one transmitter and at least one receiver 

operating on a dedicated signal channel which is used by other DS-CDMA operators in a non- 
synchronous way. A receiver for use in such a DS-CDMA system includes a data aided multiuser 
interference canceller (38) comprising a feedforward filler (48) connected to a subtracting stage (50). A 
decision stage (52) is connected to an output of the subtracting stage (50). Decisions from the decision 
stage (52) are stored in a memory which also stores preliminary decisions made on the signal applied to the 
feedforward filter (48) by a simple receiver (46). These decisions are applied to a feedback filter (58) the 
output of which is connected to the subtracting stage (50) thereby removing the effects of pre-cursive and 
post-cursive interference. The input to, and output from, the decision stage 52 are compared and the 
difference is used for adjusting the tap weights of the feedforward and feedback filters (48,58). In the case 
of a single input, single output, the feedforward filter (48) comprises a fractionally spaced equalizer 
which... Claims: 1. A direct sequence-code division multiple access (DS-CDMA) system comprising a 
transmitter and a receiver, the transmitter having means for multiplying a modulating signal with a 
spreading sequence in order to form a spread spectrum signal and for transmitting the spread spectrum 
signal on a CDMA channel, and the receiver comprising means for receiving signals on said CDMA 
channel, means for cancelling multiple access interference, said means including feedforward filtering 
means, feedback filtering means, quantizing means for producing quantized decisions from information 
derived from the feedforward filtering means, means for making tentative decisions on received 
information,means for producing an inversion of the equivalent CDMA channel in dependence on the 
cross correlation values of all the contemporaneous users and means responsive to the production of the 
inversion of the equivalent CDMA channel for adjusting tap weights of at least the feedforward filtering 
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Signal design to transmit binary information of unequal importance using asymmetric modulation 
scheme to transmit less important information at lower power level Original Titles:Method and device 

for transmitting confirmation signal Signal design to transmit binary information of unequal 

importance Conception de signaux pour transmission d'information binaire d'importance inegale... 

...APPARATUS FOR TRANSMITTING ACKNOWLEDGMENT SIGNAL AND ITS METHOD... 

...APPARATUS FOR TRANSMITTING ACKNOWLEDGEMENT SIGNAL AND ITS METHOD A 

response signal transmission apparatus and its method Method and apparatus for transmitting 

acknowledgement signals Method and apparatus for transmitting acknowledgement signals... 

...Method and apparatus for transmitting acknowledgement signals METHOD AND APPARATUS 

FOR TRANSMITTING ACKNOWLEDGEMENT SIGNALS... Alerting Abstract ...NOVELTY - A 
response signal from a receiver can be transmitted by adjusting the transmission power according to the 
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transmit power obtained when the response signal is ACK or NACK, I.e. of higher or lower importance 
and the detection threshold of the receiver is... DESCRIPTION - INDEPENDENT CLAIMS are included 
for a signal transmission apparatus, for a method of avoiding errors and for a method of transmitting 

ACK/NACK signals USE - Transmitting ACK/NACK signals in mobile communications... 

...ADVANTAGE - Increased reliability and power of transmitted ACK signals... Title Terms .../Index 
Terms/ Additional Words: TRANSMIT; Class Codes International Patent Classification IPC Class Level 
Scope Position Status Version ...Original Abstracts:wireless communication system by using automatic 
repeat pray (ARQ) plan, especially claims a response (ACK/NACK) signal power control method for using 
different powers to transmit ACK/NACK signal according to the necessary allowable probability for 

receiving error of the system. When sending and receiving the ACK/NACK signal in the repeat pray 

(ARQ) system, because emissive powers of the signal can be differentiated according to each ACK/NACK 
signal and the ACK/NACK signal is transmitted according to the emissive power satisfying the necessary 
error rate of each signal, it can reduce the necessary power consumption for transmitting the NACK 

signal A method to construct and detect a binary transmit signal is disclosed. The method is based on 

the assumption that certain decision errors are less severe than other decision errors, as e.g. in transmission 
schemes employing ARQ. To this purpose, the important information (e.g. the acknowledgement ACK) is 
transmitted with high power, whereas the less important information (e.g. the non-acknowledgement 
NACK) is transmitted with reduced power. An asymmetric modulation scheme is obtained. The detection 
threshold to decide at the receiver which of the two signals 
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Node apparatus in universal mobile telecommunication system, determines down link data channel 
power level, when modulation technique applied to downlink channel is high order modulation 
technique Original Titles: Apparatus and method for transmitting/receiving uplink transmission power 

deviation and power level APPARATUS AND METHOD FOR TRANSMITTING AND 

RECEIVING UPLINK TRANSMISSION POWER OFFSET AND HIGH-SPEED DOWNLINK 
SHARED CHANNEL POWER LEVEL IN COMMUNICATION SYSTEM EMPLOYING HIGH-SPEED 

DOWNLINK PACKET ACCESS SYSTEM Apparatus and method for transmitting/receiving uplink 

transmission power offset and HS-DSCH power level in a communication system employing HSDPA... 
...Apparatus and method for transmitting/receiving uplink transmission power offset and HS-DSCH 
power level in a communication system employing HSDPA Alerting Abstract ...NOVELTY - The node 
apparatus has determiner that determines an uplink power offset according to the comparison of different 
between signal-to-interference ratio (SIR) of uplink dedicated channel signal and preset target (SIR). 
Another determiner determines the downlink data channel power level, when the modulation technique 
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applied to the downlink data channel is a high order modulation technique. ... uplink power offset and 
downlink data channel power level transceiver apparatus; node apparatus controlling method; uplink 
power offset and downlink power level transmission/reception method; uplink transmission power 
controlling method; uplink transmission power controlling apparatus; downlink data channel power level 
transmission method; downlink data channel power level transmission apparatus; downlink data channel 
power level transmission/reception method; and downlink data channel power level transceiver 

apparatus USE - Node apparatus for transmitting uplink power offset and downlink data channel 

power level is universal mobile telecommunication system, employing high speed down-link packet access 

and used in code division multiple access (CDMA) communication system ADVANTAGE - Control 

transmission power of uplink high speed dedicated physical control channel (HS-DPCCH) in 

communication system employing high speed downlink packet access (HSDPA DESCRIPTION OF 

DRAWINGS - The figure shows the channel allocation scheme for user equipment (UE) located in a soft 
handover region in the communication system. Original Publication Data by 

Authority ArgentinaPublication No. Original Abstracts: The invention claims a high-speed packet data 
communication system, wherein a signal channel status determiner measures signal-to-interference ratio 
(SIR) of a first special uplink channel signal received from UE, and calculates difference between the SIR 
measured and a pre- set target SIR; a transmission power determiner is to compare the difference and the 
pre- set threshold and determine the uplink power deviation of a second special uplink signal channel 
applied to transmit control information for the packet data received by UE based on the comparison result; 

and a transmitter sends the determined uplink power deviation to the UE by means of a downlink Ein 

Paketdatenkommunikationssystem hoher Geschwindigkeit, wobei ein Kanalzustandsbestimmer ein Signal- 
zu-Interferenz-Verhaltnis (SIR) eines Signals eines ersten gewidmeten Aufwartsverbindungskanals misst, 
welches von einem UE empfangen wird, und eine Differenz zwischen dem gemessenen SIR und einem 

Ziel A high-speed packet data communication system, wherein a channel condition determiner 

measures a signal-to-interference ratio (SIR) of a first uplink dedicated channel signal received from a 
UE, and calculates a difference between the measured SIR and a preset target SIR; a transmission power 
determiner compares the difference with preset thresholds, and determines an uplink power offset to be 
applied to a second uplink dedicated channel transmitting control information for packet data received at 
the UE according to a result of the comparison; and a transmitter transmits over a downlink the 

determined uplink power offset to the UE A high-speed packet data communication system, wherein a 

channel condition determiner measures a signal-to-interference ratio (SIR) of a first uplink dedicated 
channel signal received from a UE, and calculates a difference between the measured SIR and a preset 
target SIR; a transmission power determiner compares the 
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Data processing apparatus for radio communication system, filters output of multipliers, before 
summation and subtraction, to remove noise Alerting Abstract ...NOVELTY - Multipliers (31- 
33,36,37) multiply long pseudo noise (PN) code respectively with in-phase and quadrature pseudo noises 
(Ipn,Qpn), pilot, dedicated control, voice and high speed data transmission channels (PI,DC,FD,SM). 
The output of the multipliers are filtered for removing noise, after which summation and subtraction are 
performed to produce the final result. ...ADVANTAGE - The data processing speed of the radio channel is 

improved and high quality data is output 41,40 Output channels Title Terms .../Index 

Terms/Additional Words: SUBTRACT; NOISE Class Codes International Patent Classification IPC 

Class Level Scope Position Status Version Original Publication Data by Authority ArgentinaPublication 
No. ...Original Abstracts:processing data using the complex quadrature phase shift keying (CQPSK) is 
disclosed. The present invention includes plurality of multiplying units for multiplying a long pseudo noise 
(PN) code with an in-phase pseudo noise and quadrature pseudo noise, and a results of multiplying with a 
pilot channel, a dedicated control channel, a voice channel, a high-speed data transmission channel, a 
result multiplied with channels with gain, filtering units for eliminating noise, subtracting and adding 
units for calculating final result, and in-phase/quadrature channel outputting unit for outputting I- 
channel/Q-channcl. The present invention can improve data processing speed of a radio channel, increase 
capacity of data processing in radio communication and also can be applicable to various kinds of radio 

communication devices transceiving voice, media contents and processing data using the complex 

quadrature phase shift keying (CQPSK) is disclosed. The present invention includes plurality of 
multiplying units for multiplying a long pseudo noise (PN) code with an in-phase pseudo noise and 
quadrature pseudo noise, and a results of multiplying with a pilot channel, a dedicated control channel, a 
voice channel, a high-speed data transmission channel, a result multiplied with channels with gain, 
filtering units for eliminating noise, subtracting and adding units for calculating final result, and in- 
phase/quadrature channel outputting unit for outputting I-channel/Q-channel. The present invention can 
improve data processing speed of a radio channel, increase capacity of data processing in radio 
communication and also can be applicable to various kinds of radio communication devices transceiving 

voice, media contents and processing data using the complex quadrature phase shift keying (CQPSK) 

is disclosed. The present invention includes plurality of multiplying units for multiplying a long pseudo 
noise (PN) code with an in-phase pseudo noise and quadrature pseudo noise, and a results of multiplying 
with a pilot channel, a dedicated control channel, a voice channel, a high-speed data transmission 
channel, a result multiplied with channels with gain, filtering units for eliminating noise, subtracting and 
adding units for calculating final result, and in-phase/quadrature channel outputting unit for outputting I- 
channel/Q-channel. The present invention can improve data processing speed of a radio channel, increase 
capacity of data processing in radio communication and also can be applicable to various kinds of radio 

communication devices transceiving voice, media contents and Claims processing data in a radio 

communication system based on a complex quadrature phase shift keying, comprising: a first multiplying 
means for multiplying a long pseudo noise (PN) code and multiplying with an in-phase pseudo noise (Ipn) 
and quadrature pseudo noise (Qpn) thereby generating a long in-phase PN code and a long quadrature PN 
code; a second multiplying means for multiplying the long in-phase pseudo PN code with a pilot channel 
(PI), a dedicated control channel (DC), a voice channel (FD) and a high speed data transmission channel 
(SM); a third multiplying means for multiplying the long quadrature PN code with a pilot channel (PI), a 
dedicated control channel (DC), a voice channel (FD) and a high speed data transmission channel (SM); 
first and second filtering means for filtering the result from the second and the third multiplying means to 
eliminate a noise controlling gains of filtered signals; fourth and fifth multiplying means for multiplying 

each of filtered signal from the first and second filtering means with a to an eighth gain-controlled 

signals, wherein the first, the second, the seventh, the eighth gain controlled signals are a gain-controlled 
in-phase pilot channel signal, a gain-controlled in-phase dedicated control channel, a gain-controlled in- 
phase voice channel and a gain-controlled in-phase high-speed data transmission channel signals, and the 
third, the fourth, the fifth, the sixth gain controlled signals area gain-controlled quadrature pilot channel 
signal, a gain-controlled quadrature dedicated control channel, a gain-controlled quadrature voice channel 
and a gain-controlled quadrature high-speed data transmission channel signals; a subtracting means for 
calculating a first result by adding the first gain-controlled signal and the second gain controlled signal, 
calculating a second result by adding the third gain-controlled signal and the fourth gain-controlled signal, 
and calculating a difference between the first and the second results; an adding means for calculating a 

third result by adding the fifth and sixth gain-controlled signals, calculating fourth result by adding the 

seventh and eighth gain-controlled signals and calculating a sum of the third and the fourth results; an in- 
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phase channel outputting means for outputting the result outputted from the subtracting means as an in- 
phase channel (I-channel) signal; and a quadrature channel outputting means for outputting the result 

outputted from the adding means as a quadrature channel (Q-channel) signal processing data in a 

radio communication system based on a complex quadrature phase shift keying, comprising: a first 
multiplying means for multiplying a long pseudo noise (PN) code and multiplying with an in-phase pseudo 
noise (Ipn) and quadrature pseudo noise (opn) thereby generating a long in-phase PN code and a long 
quadrature PN code; a second multiplying means for multiplying the long in-phase pseudo PN code with a 
pilot channel 
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Channel estimation method in antenna diversity system, involves determining channel estimate by 
rotating common pilot channel according to estimated phase difference between common and 
dedicated pilot channels ...Original Titles: Antenna phase estimation algorithm for WCDMA closed loop 

transmitter antenna diversity system METHOD FOR INFERRING AND DECIDING CHANNEL 

FOR WCDMA CLOSED LOOP TRANSMITTER ANTENNA DIVERSITY SYSTEM Antenna 

phase estimation algorithm for WCDMA closed loop transmitter antenna diversity system Alerting 
Abstract ...NOVELTY - The phase difference between a common pilot channel and a dedicated pilot 
channel received from an antenna of a diversity system, is computed as a complex coefficient between the 
channels. A channel estimate is determined by rotating a common pilot channel according to the 
estimated phase difference. USE - For enabling a user equipment (UE) e.g. mobile phone, vehicle 
telephone to perform channel estimation for channel received from base station (BS) in antenna diversity 
system, wideband CDMA telecommunication system, transmit antenna diversity system... 
...ADVANTAGE - Does not require an estimation of the noise variable or the relative powers of the 

dedicated and common pilot channels DESCRIPTION OF DRAWINGS - The figure shows the 

block diagram of the 3G- WCDMA base station transmitter arrangement. Title Terms /Index 

Terms/ Additional Words: CHANNEL; Class Codes International Patent Classification IPC Class Level 

Scope Position Status Version Date "Version 7" H04B-0001/707 H04B-0007/06 H04B-0007/08... 

...H04B-0007/26 H04L-0025/02 H04B-0001/707 H04B-0007/04 H04B-0007/08 H04B- 

0007/26 H04L-0025/02 Original Publication Data by Authority ArgentinaPublication No. Original 

Abstracts: This invention teaches a method for use in a wideband CDMA telecommunications system 
wherein the base station (BS) utilizes transmit antenna diversity. The method enables the user equipment 
(UE) to perform a channel estimate for channels received from two antennas of the base station (BS), and 
includes steps of: (a) computing an estimated phase difference between a common pilot channel and a 
dedicated pilot channel received from one antenna of the two antennas; and (b) determining the channel 
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estimate by rotating a common pilot channel estimate according to the estimated phase difference, 

where the phase difference is computed as a complex coefficient between the common pilot channel and 
the dedicated pilot channel. In a second aspect this invention teaches a method for use in a transmit 
antenna diversity system for enabling the UE to verify an antenna phase shift previously signaled by the 

UE to the BS. This method includes a first step of: (a This invention teaches a method for use in a 

wideband CDMA telecommunications system wherein the base station (BS) utilizes transmit antenna 
diversity. The method enables the user equipment (UE) to perform a channel estimate for channels 
received from two antennas of the base station (BS), and includes steps of: (a) computing an estimated 
phase difference between a common pilot channel and a dedicated pilot channel received from one 
antenna of the two antennas; and (b) determining the channel estimate by rotating a common pilot 
channel estimate according to the estimated phase difference, where the phase difference is computed as 
a complex coefficient between the common pilot channel and the dedicated pilot channel. In a second 
aspect this invention teaches a method for use in a transmit antenna diversity system for enabling the UE 
to verify an antenna phase shift previously signaled by the UE to the BS. This method includes a first step 
of: (a) exploiting prior knowledge of BS rotation angle by considering a... Claims:In an antenna diversity 
system, a method for enabling a user equipment (UE) to perform a channel estimate for a channel 
received from a base station (BS), comprising steps ofxomputing an estimated phase difference between a 
common pilot channel and a dedicated pilot channel received from one antenna of the antenna diversity 
system; anddetermining the channel estimate by rotating a common pilot channel estimate according to 
the estimated phase difference, wherein the phase difference is computed as a complex coefficient 

between the common pilot channel and the dedicated pilot channel What is claimed is:l. In an 

antenna diversity system, a method for enabling a user equipment (UE) to perform a channel estimate for a 
channel received from a base station (BS), comprising steps of: computing an estimated phase difference 
between a common pilot channel and a dedicated pilot channel received from one antenna of the antenna 
diversity system; anddetermining the channel estimate by rotating a common pilot channel estimate 
according to the estimated phase difference, wherein the phase difference is computed as a complex 
coefficient between the common pilot channel and the dedicated pilot channel.> 
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Cross talk interference reducing method for transmission system that uses frequency-translated 
signals, involves multiplying complex coupling factor with signal on second line ...Original 

Titles:METHOD AND DEVICE FOR REDUCING CROSSTALK INTERFERENCE Method and 

device for reducing crosstalk interference METHOD AND DEVICE FOR REDUCING 

CROSSTALK INTERFERENCE Alerting Abstract DESCRIPTION - A complex coupling factor for 
cross talk interference is estimated. The estimated factor is multiplied with the signal on the second line to 
estimate the induced cross talk interference, after which the estimated induced cross talk is subtracted 
from the signal on the first line. An INDEPENDENT CLAIM is also included for a transmission system... 

...USE - For transmission system that uses frequency-translated signals ADVANTAGE - Simplifies 

and accelerates computation of the estimated cross talk interference, thus permitting adaptive reduction 
of cross talk interference in real time to accommodate e.g. changing transmitting conditions. Allows, at 

least partly, implementation using dedicated hardware for improved speed DESCRIPTION OF 

DRAWINGS - The figure shows the conceptual diagram of the cross talk reduction in a transmission 

line. Title Terms .../Index Terms/Additional Words: INTERFERENCE; TRANSMISSION; ... 

...FREQUENCY; Class Codes International Patent Classification IPC Class Level Scope Position Status 

Version Date H04B-015/00 H04L-027/26 Main "Version 7" H04B-0001/10 H04B-0003/06... 

...H04B-0003/32 H04B-0003/46 H04B-0001/10 H04B-0003/02 H04B-0003/06 H04B- 

0003/46 Original Publication Data by Authority ArgentinaPublication No. Original Abstracts:In a 
transmission system which makes use of frequency translated signals, said system having a plurality of 

lines, a method for reducing crosstalk induced on a signal on a first line of said plurality of lines by a 

lines, comprising estimating a complex coupling factor for the crosstalk, which when multiplied by the 
signal on the second line estimates the induced crosstalk, and subtracting the estimated induced crosstalk 
from the signal on the first line. The invention comprises an approximation method for multiplying the 

complex coupling factor by the signal on the second In a transmission system which makes use of 

frequency translated signals, said system having a plurality of lines, a method for reducing crosstalk 

induced on a signal on a first line of said plurality of lines by a signal on a second line lines, 

comprising estimating a complex coupling factor for the crosstalk, which when multiplied by the signal on 
the second line estimates the induced crosstalk; and subtracting the estimated induced crosstalk from the 
signal on the first line. The invention comprises an approximation method for multiplying the complex 

coupling factor by the signal on the second line that operates on the In a transmission system which 

makes use of frequency translated signals, said system having a plurality of lines, a method for reducing 
crosstalk induced on a signal on a first line of said plurality of lines by a signal on a second line of said 
plurality of lines, comprising estimating a complex coupling factor for the crosstalk, which when multiplied 
by the signal on the second line estimates the induced crosstalk, and subtracting the estimated induced 
crosstalk from the signal on the first line. The invention comprises an approximation method for 
multiplying the complex coupling factor by the signal on the second line that operates on the signal on the 

second line L'invention concerne un systeme de transmission utilisant des signaux transposes en 

frequence, ledit systeme possedant plusieurs lignes, et un procede permettant d'attenuer les perturbations 
diaphoniques provoquees sur un signal d'une premiere ligne par un signal d'une seconde ligne. Ce procede 
consiste a estimer un coefficient de couplage complexe pour les perturbations diaphoniques qui... 
Claims: What is claimed is: 1. In a transmission system which makes use of frequency translated signals, 
particularly a discrete multi tone (DMT) modulated transmission system or an OFDM (Orthogonal 
Frequency Division Multiplex) transmission system, having a plurality of lines and in which the 
modulation may be effected using a fast inverse Fourier transform (IFFT), a method for reducing crosstalk 
interference induced on a signal on a first line of said plurality of lines by a signal on a second line of said 
plurality of lines, comprising estimating a complex coupling factor for the crosstalk interference, which 
when multiplied by the signal on the second line estimates the induced crosstalk interference; and 
subtracting the estimated induced crosstalk interference from the signal on the first line; further 
comprising multiplying the complex coupling factor by the signal on the second line through an 
approximation method operating on the signal on the second... 
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Controlling power transmission level with a cellular communications system Alerting Abstract 

...NOVELTY - Base stations (12, 14, 16) transmit a pilot channel that includes pseudo-random noise and 
information about the power of the transmitted signal. The mobile (18) measures the ratio of signal to 
noise on the pilot channel and from this makes an estimate of the strength of the received signal. The 
mobile compares this with the strength of its communications channel and transmits the ratio of the 
difference between the two back to the central control. ...ADVANTAGE - Does not need a dedicated 
control channel. Improved assessment of the quality of reception... Title Terms .../Index Terms/ Additional 
Words: TRANSMISSION; Class Codes International Patent Classification IPC Class Level Scope 

Position Status Version Date H04B-007/005 H04B-007/26 Main "Version 7" H04B-001/04 H04B- 

0001/04 H04B-0007/00 H04B-0007/005 H04B-0007/005 H04B-0007/26 H04B- 

0007/26 H04L-0001/00 H04B-0001/04 H04B-0007/00 H04B-0007/005 H04B -0007/26... 

...H04L-0001/00 Original Publication Data by Authority ArgentinaPublication No. Original Abstracts:A 
system and method is taught for controlling the power level of transmissions within a communication 
system having a base station, a mobile station, a communication channel, and a pilot channel. The mobile 
station determines a signal strength value according to a communication signal received by way of the 
communication channel. A pilot channel signal is determined according to a pilot signal transmitted by 
way of the pilot channel. The signal to noise ratio of the communication signal is determined according to 
the determined signal strength value and the pilot channel signal. The power level of a transmission is 
controlled according to the signal to noise ratio. The noise level in a communication channel within the 
communication system is estimated. The pilot channel signal includes pilot energy and pilot noise 
components. The pilot energy component is removed to provide a remaining pilot signal. Communication 
system operations are controlled according to the remaining pilot signal. The power levels of 
transmissions are controlled by determining the signal to noise ratio of a signal within the communication 
channel and determining a difference signal. The difference signal is formed by determining the 
difference between determined and desired signal to noise ratios. The difference signal is transmitted 
between the base station and the mobile station. The pilot channel has at least one frame and the power 
control signal is inserted into the frame. Thus, information representing the strength of the communication 
signal is transmitted to the base by way of the pilot channel within the frame.The pilot channel can have 

two information frames for transmitting the power control signal a plurality of times A system and 

method is taught for controlling a power level of transmissions in a communication system. A first station 
of the communication system transmits signals on a communication channel and a pilot channel. A 
second station of the communication system measures a signal to noise ratio of the received signal, 
determines an expected signal to noise ratio of the received signal, determines an indicator in accordance 
with said measured signal to noise ratio and the expected signal to noise ratio, and transmits the indicator. 
The indicator is received at the first station the indicator and determines the power level of the next 
transmission according to said received indicator. Additionally, methods and systems for the measurement 
of the signal to noise ratio of the received signal, the determination of the expected signal to noise ratio of 

the received signal, and the determination of the indicator are disclosed A system and method is taught 

for controlling a power level of transmissions in a communication system. A first station of the 
communication system transmits signals on a communication channel and a pilot channel. A second 
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station of the communication system measures a signal to noise ratio of the received signal, determines an 
expected signal to noise ratio of the received signal, determines an indicator in accordance with said 
measured signal to noise . . . 
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RF interference mitigation method for multicarrier modulation estimating RF interference in 

accordance with frequency domain model for RF interference and estimated frequency of RF 
interference and removing estimated RF interference from frequency domain data Original 

Titles:SCHALTUNG ZUR DIGIT ALEN HF-INTERFERENZUNTERDRUCKUNG DIGITAL 

RADIO FREQUENCY INTERFERENCE CANCELLER SUPPRESSEUR D INTERFERENCES 

SUR DES TRANSMISSIONS RF NUMERIQUES SCHALTUNG ZUR DIGIT ALEN HF- 
INTERFERENZUNTERDRUCKUNG DIGITAL RADIO FREQUENCY INTERFERENCE 

CANCELLER SUPPRESSEUR D INTERFERENCES SUR DES TRANSMISSIONS RF 

NUMERIQUES DIGITAL RADIO FREQUENCY INTERFERENCE CANCELLER The 

relief Imoderation of the radio-frequency interference in a multicarrier transmission system Digital 

radio frequency interference canceller Digital radio frequency interference canceller Digital 

radio frequency interference canceller. Mitigating radio frequency interference in multi-carrier 

transmission systems Digital radio frequency interference canceller Digital radio frequency 

interference canceller DIGITAL RADIO FREQUENCY INTERFERENCE CANCELLER 

Alerting Abstract ...The method involves obtaining frequency domain data which are associated with a 
frequency band. A restricted frequency sub-band is identified in the frequency band and the frequency of 

the RF interference sub-band is estimated in the restricted frequency sub-band The RF interference is 

estimated in accordance with a frequency domain model for the RF interference and the estimated 
frequency of the RF interference. The estimated RF interference is removed from the frequency domain 

data USE/ADVANTAGE - E.g. for ADSL, VDSL, FTTC. Eliminates or compensates radio 

interference and reduces sidelobe interference. Title Terms .../Index Terms/Additional Words: 

INTERFERENCE; FREQUENCY; Class Codes International Patent Classification IPC Class Level 

Scope Position Status Version Date Original Publication Data by Authority ArgentinaPublication No. 
Original Abstracts:Disclosed are radio frequency (RF) interference cancellation techniques that 
effectively estimate RF interference to the data signals being received using a frequency domain model, 
and then remove the estimated RF interference from the received data signals. Improved techniques for 
digitally filtering multicarrier modulation samples to reduce sidelobe interference due to the RF 
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interference are also disclosed The radio frequency interference used as the seriouslsignificant 

disorderldamagelfailure for a receiver which can receive a transmission-data signal appropriately is 
removed, or it is compensated. The radio-frequency (RF) interference to the data signal received is 
effectively estimated using a frequency-domain model, Subsequently, the radio-frequency (RF) 
interference denial technique of removing this estimated RF interference from a receiving data signal is 
disclosed. In order to reduce the side-lobe interference by RF interference, the improvement technology 
which carries out the wave filter of the multi-careerlcarrier modulationlalteration sample digitally is also 
disclosed.FIG. 7This invention relates to a radio-frequency (RF) interference denial (erasurelelimination) 
device.More Specifically, It is related with RF interference denial in a multi-careerlcarrier transmission 

system Disclosed are radio frequency (RF) interference cancellation techniques that effectively 

estimate RF interference to the data signals being received using a frequency domain model, and then 
remove the estimated RF interference from the received data signals. Improved techniques for digitally 
filtering multicarrier modulation samples to reduce sidelobe interference due to the RF interference are 

also disclosed Disclosed are radio frequency (RF) interference cancellation techniques that 

effectively estimate RF interference to the data signals being received using a frequency domain model, 
and then remove the estimated RF interference from the received data signals. Improved techniques for 
digitally filtering multicarrier modulation samples to reduce sidelobe interference due to the RF 

interference are also disclosed Disclosed are radio frequency (RF) interference cancellation 

techniques that effectively estimate RF interference to the data signals being received using a frequency 
domain model, and then remove the estimated RF interference from the received data signals. Improved 
techniques for digitally filtering multicarrier modulation samples to reduce sidelobe interference due to 

the RF interference are also disclosed Methods and mechanisms for reducing the impact of radio 

frequency interference in multi-carrier data transmission systems are disclosed. In one aspect of the 
invention, dummy tones are used to suppress emissions within restricted frequency bands. In another 
aspect of the invention, signals received on designated quiet tones are detected and used to facilitate 

canceling of RF interference on adjacent active tones Disclosed are radio frequency (RF) 

interference cancellation techniques that effectively estimate RF interference to the data signals being 
received using a frequency domain model, and then remove the estimated RF interference from the 
received data signals. Improved techniques for digitally filtering multicarrier modulation samples to reduce 
sidelobe interference due to the RF interference are also 
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Alerting Abstract ...NOVELTY - Several base stations are communicated via three types of frequency 
bands which utilize three different spreading codes respectively. The spreading codes (502) define 
corresponding channels utilized for transmitting data and control information from mobile station to base 
stations and transmitting information from base station respectively. DESCRIPTION - A mobile station 
sends communication signal which is received via a channel by which location of mobile station based on 
communication is traced. A time difference of arrival (TDOA) or time of arrival (TO A) technique is used 
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for locating the mobile station. The mobile station broadcasts a probe sequence corresponding via a 
Dedicated Physical Location Channel (DPLCH) (501). An INDEPENDENT CLAIM is also included for 

an apparatus for locating mobile station within a communication system ADVANTAGE - A high signal 

to noise ratio at a distant base station without an undesirable mobile station transmit power envelope 
variation is obtained by utilizing spreading codes with orthogonal channelization codes and mobile station 

having Dedicated Physical Located Channel DESCRIPTION OF DRAWINGS - The figure shows 

the block diagram of mobile station spreading technique for wide band CDMA system incorporating a 

location sub-channel 501 Dedicated Physical Location Channel Class Codes International Patent 

Classification IPC Class Level Scope Position Status Version Date H04B-007/26 Main H04B-0007/26... 

...H04J-0011/00 H04J-0013/00 H04B-0007/26 H04J-0011/00 H04J-0013/00 Original 

Publication Data by Authority ArgentinaPublication No. Original Abstracts: A Dedicated Physical 
Location Channel (DPLCH) is utilized by a mobile station (913) to support subscriber location functions. 
The DPLCH is spread by an Orthogonal Variable Spreading Factor (OVSF) code CL (502) of length 256 
which is distinct from those OVSF codes assigned to other channels utilized by the mobile station (913). 
When a power-up function (PUF) is received by the mobile station (913), the DPLCH sub-channel 
amplitude is then modified relative to the other channels being utilized by the mobile station (913) using 
gain module GL (503) prior to combination (504) with the other channels. A Dedicated Physical Location 
Channel (DPLCH) is utilized by a mobile station (913) to support subscriber location functions. The 
DPLCH is spread by an Orthogonal Variable Spreading Factor (OVSF) code CL (502) of length 256 which 
is distinct from those OVSF codes assigned to other channels utilized by the mobile station (913). When a 
power-up function (PUF) is received by the mobile station (913), the DPLCH sub-channel amplitude is 
then modified relative to the other channels being utilized by the mobile station (913) using gain module 
GL (503) prior to combination (504) with the other channels. A Dedicated Physical Location Channel 
(DPLCH) is utilized by a mobile station (913) to support subscriber location functions. The DPLCH is 
spread by an Orthogonal Variable Spreading Factor (OVSF) code CL (502) of length 256 which is 
distinct from those OVSF codes assigned to other channels utilized by the mobile station (913). When a 
power-up function (PUF) is received by the mobile station (913), the DPLCH sub-channel amplitude is 
then modified relative to the other channels being utilized by the mobile station (913) using gain module 
GL (503) prior to combination (504) with the other channels. Selon...Claims:method for mobile station 
location within a communication system, the method comprising the steps of: communicating with a 
plurality of base stations via a first frequency band using a first spreading code, the first spreading code 
defining a first channel utilized for transmitting data from the mobile station to the plurality of base 
stations;communicating with the plurality of base stations via the first frequency band using a second 
spreading code, the second spreading code defining a second channel utilized for transmitting control 
information from the mobile station to the plurality of base stations; andcommunicating with the plurality 
of base stations via the first frequency band using a third spreading code, the third spreading code 
defining a third channel utilized for transmitting information to the plurality of base stations that is 
utilized for mobile station location within the communication system. 
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Communication system has lattice interfered channel that has underlying channel with at least one 
lattice vector on input or output side, and through which mapped input data are sent at rate of not 
more than 0.7 bit per dimension 
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Communication system has lattice interfered channel that has underlying channel with at least one 
lattice vector on input or output side, and through which mapped input data are sent at rate of not 
more than... Original Titles:Practical coding and metric calculation for the lattice interfered channel ... 
...Practical coding and metric calculation for the lattice interfered channel ...Inventor: SUTSKOVER I 
Alerting Abstract ...NOVELTY - An encoder-modulator maps the input data to a signal constellation 
having signal points in predetermined dimensions. A lattice interfered channel (100) has an underlying 
channel (102) with an addition of at least one lattice vector (104,106) at the input or output side. The 
mapped input data are sent through the lattice interfered channel at a rate of not more than 0.7 bit per 
dimension. ... a transmitters data transmission method in communication system;a signal processing 

method in communication system;a receiver; anda communication device ADVANTAGE - Ensures 

effective transmission of predetermined data from a transmitter to a receiver through effective use and 
implementation of lattice interfered channel. DESCRIPTION OF DRAWINGS - The figure shows the 
block diagram illustrating lattice interfered channel. lOOLattice interfered channel 102Underlying 
channelTitle Terms .../Index Terms/Additional Words: INTERFERENCE; CHANNEL; Class Codes 
Original Publication Data by Authority ArgentinaPublication No. ...Inventor name & address: Sutskover, 

Ilan Sutskover, Ilan Original Abstracts: A communication system includes a lattice interfered 

channel to transmit data from a transmitter to a receiver. In one embodiment, an encoding-modulation 
scheme having a rich signal constellation is used to encode data at a total rate of 0.7 bit/dimension or less 
before transmission into the lattice interfered channel. In another embodiment, decoding metric 
approximation techniques are used to process signals received from the lattice interfered channel. In still 
another embodiment, multilevel code (MLC) decoding is used to jointly decode an error correction code 

component and a lattice/MLC component of a received signal. The error A communication system 

includes a lattice interfered channel to transmit data from a transmitter to a receiver. In one 
embodiment, an encoding- modulation scheme having a rich signal constellation is used to encode data at a 
total rate of 0.7 bit/dimension or less before transmission into the lattice interfered channel. In another 
embodiment, decoding metric approximation techniques are used to process signals received from the 
lattice interfered channel. In still another embodiment, multilevel code (MLC) decoding is used to jointly 
decode an error correction code component and a lattice/MLC component of a received signal. The error 

correction code component may then Claims: system comprising:an encoder-modulator to map input 

data into a signal constellation having more than 2k signal points in k dimensions; anda lattice interfered 
channel having an underlying channel together with at least one of the following: addition of at least one 
lattice vector at an input of the underlying channel and addition of at least one lattice vector at an output of 
the underlying channel; wherein said mapped input data is transmitted through said lattice interfered 

channel at a rate of no greater than 0.7 bits/dimension What is claimed is: 13. A method for use in 

processing a received signal in a communication system having a lattice interfered channel, said lattice 
interfered channel having an associated lattice Lambda, comprising: identifying a point in the lattice 
Lambda that is closest to a received vector; obtaining a difference between the identified point and the 
received vector; andapproximating a bit likelihood measure using the difference between the identified 
point and the received vector, wherein approximating includes calculating a sum of probabilities associated 
with a subset of lattice points within lattice Lambda, wherein prior to transmission, input data is mapped 
into a signal constellation having more than 2k signal points in k dimensions for transmission through the 
lattice interfered channel at a rate of no greater than 0.7 bits/dimension 
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Transmission power control method for communication system e.g. CDMA, involves setting power of 
downlink channel for mobile station based on reception quality information when time difference is 
equal or less than preset threshold 
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Alerting Abstract ...that the transmission power is adapted for positional relationships among mobile 
stations and a radio state of each mobile station, hence prevents unnecessary power consumption, reduces 
interference to other mobile stations, and increase a downlink capacity... Class Codes International Patent 
Classification IPC Class Level Scope Position Status Version Date H04B-007/005 Main H04B -000 1/00... 

...H04B-0001/69 H04B-0007/005 H04B-0007/005 H04B-0007/26 H04B-0007/26 H04B- 

0001/00 H04B-0001/69 H04B-0007/005 H04B-0007/005 H04B-0007/005 H04B- 

0007/26 H04B-0007/26 Original Publication Data by Authority ArgentinaPublication No. Original 

Abstracts:The object of the present invention to provide a transmission power control method which 
controls a transmission power of a downlink common channel so that the transmission power is adapted 
to a radio state of each mobile station, and thereby makes it possible to prevent unnecessary power 
consumption, reduce interference to other mobile stations, and increase a downlink capacity. The 
transmission power control method includes (A) acquiring reception quality information of a downlink for 
a mobile station (S 103), (B) measuring a time difference (S 104) between when an uplink signal is received 
from the mobile station and when a downlink signal to be transmitted to the mobile station occurs, (C) 
comparing the measured time difference with a predetermined threshold (S105), and (D) setting a 
transmission power (SI 06) of a downlink common channel for mobile stations in accordance with the 
reception quality information, when the time difference is equal to or less than the predetermined 

threshold The object of the present invention to provide a transmission power control method which 

controls a transmission power of a downlink common channel so that the transmission power is adapted 
to a radio state of each mobile station, and thereby makes it possible to prevent unnecessary power 
consumption, reduce interference to other mobile stations, and increase a downlink capacity. The 
transmission power control method includes (A) acquiring reception quality information of a downlink for 
a mobile station, (B) measuring a time difference between when an uplink signal is received from the 
mobile station and when a downlink signal to be transmitted to the mobile station occurs, (C) comparing 
the measured time difference with a predetermined threshold, and (D) setting a transmission power of a 
downlink common channel for mobile stations in accordance with the reception quality information, when 

the time difference is equal to or less than the predetermined threshold The object of the present 

invention to provide a transmission power control method which controls a transmission power of a 
downlink common channel so that the transmission power is adapted to a radio state of each mobile 
station, and thereby makes it possible to prevent unnecessary power consumption, reduce interference to 
other mobile stations, and increase a downlink capacity. The transmission power control method includes 
(A) acquiring reception quality information of a downlink for a mobile station, (B) measuring a time 
difference between when an uplink signal is received from the mobile station and when a downlink signal 
to be transmitted to the mobile station occurs, (C) comparing the measured time difference with a 
predetermined threshold, and (D) setting a transmission power of a downlink common channel for mobile 
stations in accordance with the reception quality information, when the time difference is equal to or less 
than the predetermined threshold. ...Claims:power of a common channel is calculated based on the 
receiving quality of the said down-link, The transmission... 
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code division multiple access communication system by using high speed downlink packet access 
system and method therefor 
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Author Abstract: High-frequency single- sideband (HF-SSB) radio communication suffers from highly variable 
propagation conditions and high levels of interference noise. One common distortion is frequency shifting of the 
received signal due to differences in carrier frequencies at transmitter and receiver, or to Doppler shifts due to the 
transmission path. This paper presents a novel and efficient technique to compensate for this frequency offset using 
only the received speech signal to determine and correct... 

Text: 
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10892604 Word Count: 5572 

Radio tomorrow 

INTERNATIONAL DEFENSE REVIEW ( IDR ) MARCH 01, 2001 v. 034 no. 003 
Section Heading: FEATURE 
By: Rupert Pengelley 

...may be in sight. Taken to its logical conclusion, the 
advent of "software defined" radio (SDR) technology could well lead 
to the replacement of tactical transceivers by what might more 
appropriately be described as general-purpose RF (radio frequency) 
devices . 

Certainly future tactical voice and data radios are unlikely to look... 
...1980s to the next 

generation of multiband, multimode SDRs. Those cited by the US Army 
include the increased speed and power of embedded processors, 
enabling modem functions to be performed by digital signal 
processors; improved analog-to-digital converters, with higher 
conversion rates and increased dynamic operating ranges; and the 
development ... intelligence (C{4}I) capabilities to the warfighters. 
This family will cover an operating spectrum from 2-2,000MHz 
initially, and will be capable of transmitting voice, video, and 
data. . . JTRS is not a one-size-fits-all system. Rather, it is a 
family of radios that are interoperable, affordable and... 

...I radio is intended to exploit a 

CSMA-based wideband waveform being developed by BAE in association 
with BBN to provide the TOC-to-TOC transmission capability for those 
brigade combat teams belonging to the Army's second digitized 
division (1 CAV) and the third of the US Army's new... 

...five channels in its ground application, six for 
aviation, two for dismounted applications, and one for handheld 
applications.) Consideration is being given to using a common three- 
channel stackable radio for both ground and aviation applications, 
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offering enhanced economy of scale and redundancy. The Block II 
draft RFP was expected before the end... 120km point-to- 
point range) using a vertical-incidence SINCGARS-type whip antenna. 
(Though it could not receive such signals, the IWS could generate 
skywave transmissions for sending messages to base stations over 
even longer ranges, according to ITT.) 

To maintain such links "anywhere, anytime", the SUO SAS 
communication network is... 

. . .provided by embedding serial probes 

in the bit patterns. These retune digitally matched filters in the 
receiver faster than the channel decorrelat ion time, overcoming 
multipath interference from buildings and obstacles. The software 
tunable digital synthesizer automatically changes frequency (up to 
500/sec) throughout the full frequency range, using direct sequence 
spread spectrum rather than frequency hopping. The coding rate, data 
rate, power and frequency may be changed on every packet 
transmission. 

Some partisan engineers are inclined to dismiss JTRS as "a future- 
looking product whose time has passed." Even so, SUO SAS itself is 
billed as... was already a feature of several 
modern tactical radio systems and that SDRs, like their 
predecessors, embody a receiver (or two) as well as a transmitter. 
The thought therefore arises that this receive capacity could 
usefully be exploited for signal intelligence (SIGINT) as well as 
for communication. Pursuing this line of... 

...to the services at the start of FY03. Both 

teams are working on subsystem risk reduction, including co-channel 
and co-site mitigation. A key difference between ACN and standard 
rebroadcast systems is that waveform translation takes place, the 
signal being demodulated and remodulated in another waveform (and on 
another frequency. . . 

...Aerospace concerning vitreous antennas. These 
allow apertures to be stacked like pancakes on substrates, the 
material between them being electrically transparent to the other 
layers /transmitters, and giving coverage from UHF through S-band. 
MEMs are used to reconfigure the array electrically for the 
appropriate frequency, their cantilevered electrostatic design 
approach ... any type of communications, be it 

satellite, hand-held, or wideband. It will have a communications 
network routing architecture that will be independent of the 
transmission media, and will be fully adaptive in terms of quality 
of service versus network conditions and traffic demands. 

MOSAIC is described by David Keetley, director... 

...are two requirements: what type 

of service are you asking for, and what speed of service and time of 
delivery or precedence - as, say, the difference between fire orders 
or conferencing." 
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(57) ABSTRACT 

A multi-carrier transmission system has a plurality of 
transceivers, each having a receiver and a transmitter. Data 
is transmitted between the transceivers by modulating the 
data onto a multiplicity of carrier waves in the form of 
multi-bit symbols, wherein each carrier wave constitutes a 
channel, and wherein the number of bits per symbol (ihe bit 
loading), varies between channels and, within a channel, 
with lime, so that each channel has associated therewith a bit 
loading parameter. In operation, the multi-carrier system is 
adapted to synchronously update, at the transceivers, the bit 
loading parameters associated wilh each channel by trans- 
mission of data over a control channel, in that the control 
channel is established, at system start-up, oa a predeter- 
mined one of the multiplicity of carrier waves whose iden- 
tity is known to the transceivers, and in that the control 
channel is, after start-up, changed from the predetermined 
channel to a further channel, selected by the first transceiver 
on the basis of channel characteristics. 
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In certain types ol transmission system using nQ ?1 M it may 
be desirable to use a control applications in multicarrier 
snrstems with variable bitloading it is particularly desirible 
that the control channel data be separable from user data 
when the value of n. i.e the bitloading. or constellation 
identifier are unknown. Binary data is transmitted using 
nQAM modulation. The binary data comprises both user 
data and cointro,l data. This means that a control channel is 
superimposed on a user data channel within the same nQAM 
constellation. Separation of said control data is lacilitated by 
arranging control data points to have either unique phase, or 
unique amplitude, values. These unique values are not 
employed for user data points. 



Analog front-end 



Digital receiver unit 



— ^ 




Channel 
estimation 
& 

Equalizer 



Symbol L. Deinterleaving 
detection] 1 & decoding 



man ' mujatln} 



Sys tern controller/ PC I 



Symbol 
mapping 



Encoding t ^_Bif 
interleavwQ\\mQnipulQtmg 



tool 



network/ 

application 

interface 



computer 
interface 



Digital transm'iftertinit 



Usb06538986B2 



(12) United States Patent 

Isaksson et al 



(io) Patent No.: US 6,538,986 B2 
(45) Date of Patent: Man 25, 2003 



(54) DATA TRANSMISSION SYSTEM AND 

METHOD USING NQAM CONSTELLATION 
WITH A CONTROL CHANNEL 
SUPERIMPOSED ON A USER DATA 
CHANNEL 

(75) Inventors: Mikaei Isaksson, Lulea (SE); Magnus 
Johansson, Lulea (SE); Harry Tonvall, 
Lulea (SE); Lennart Olsson, Lulea 
(SE); lomas Stefansson, Lulea (SE); 
Hans Ohman, Lulea (SE); Gunnar 
Bahlenberg, Lulea (SE); Anders 
Isaksson, Lulea (SE); Goran Okvist, 
Lulea (SE); Lis-Marie Ljunggren, 
Lulea (SE); Tomas Nordstrom, Lulea 
(SE); Lars-Ake Isaksson, Lulea (SE); 
Daniel Bengtsson, Lulea (SE); Shvert 
Hakansson, Lulea (SE); Ye Wen, Lulea 
(SE) 

(73) Assignee: STMicroelectronics Amsterdam 
(NL) 



( * ) Notice: 

(21) Appl. No,: 

(22) PCT Filed: 
PCT No,: 



(86) 



Subject to any disclaimer, the term of this 
patent is extended or adjusted under 35 
U.S.C. 154(b) by 0 days. 

09/147,745 

Sep. 1, 1997 

PCT/SE97/01458 



§371 (c)(1), 
(2), (4) Date: 



Jun. 1, 1999 
(87) PCT Pub. No,: WO98/10553 
PCT Pub. Date: Man 12, 1998 
(65) Prior Publication Data 

US 2002/0126768 AI Sep. 12, 2002 
(30) Foreign Application Priority Data 

Sep. 2, 1996 (SE) 9603195 



(51) Int. Ck 7 H04J 11/00; H04J 3/12; 

H04L 23/02; H04L 5/12 

(52) U.S. CI 370/207; 375/261; 370/527 

(58) Field of Search 370/204, 206, 

370/207, 208, 464, 465, 203, 498, 522, 
527; 375/259, 261, 262, 295, 298, 260 

(56) References Cited 

U.S. PATENT DOCUMENTS 

4,627,077 A * 12/1986 Armstrong 375/261 

5 ? 022 ? 053 A * 6/1991 Chung et al 375/265 

5 3 105 ? 443 A * 4/1992 Belts et al 375/265 

(List continued on next page,) 

FOREIGN PATENT DOCUMENTS 

DE 44 35 833 5/1995 

Primary Examiner— Kwang Bin Yao 

Assistant Examiner — Kevin C. Harper 

(74) Attorney, Agent, or Firm — Allen, Dyer, Doppelt, 

Milbrath & Gilchrist, PA; Lisa K. Jorgenson 



(57) 



ABSTRACT 



A transmission system using nQAM, A control channel is 
superimposed on a nQAM data channel. It is desirable to 
separate control channel data from user data without using 
a higher level protocol layer. In multi-carrier systems with 
variable bit-loading it is particularly desirable that the con- 
trol channel data be separable from user data when the value 
of n. that is the bit-loading, or constellation identifier are 
unknown. Binary data is transmitted using nQAM modula- 
tion. The binary data includes both user data and control 
data. This means that a control channel is superimposed on 
a user data channel within the same nQAM constellation. 
Separation of the control data is facilitated by arranging 
control data points to have either a unique phase or unique 
amplitude values. These unique values are not employed for 
user data points. 

15 Claims, 14 Drawing Sheets 
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Advanced Signal Processing Algorithms for 3rd Generation Wireless 

Mobile Systems 

Alessandro Nordio 

Sintesi 

Motivata dall'interesse mondiale per i sistemi cellulari di terza generazione, questa 
tesi affronta il problema di migliorare l'affidabilita del livello fisico dei ricevitori di 
terza generazione che utilizzano la tecnica W-CDMA, proponendo algoritmi semplici 
e a bassa complessita adatti ad essere implementati su sistemi esistenti. II lavoro e 
stato sviluppato neirambito di una collaborazione tra Eurecom e EPFL il cui obiettivo 
e quello di realizzare una piattaforma Software Radio capace di convalidare algoritmi 
avanzati per il trattamento del segnale. La ricerca si focalizza su molti aspetti del 
livello fisico, tra cui sincronizzazione, multi-user detection, e il progetto di algoritmi 
per l'interfaccia tra il livello fisico e l'hardware. 

Per quanto riguarda la sincronizzatione consideriamo l'acquisizione iniziale della 
temporizzazione per i sistemi DS-CDMA, considerando canali affetti da fading e mul- 
tipath. Nel nostro modello la Stazione Base diffonde un segnale di sincronizzazione 
composto da un gruppo di chip modulati, ripetuti periodicamente e separati da lunghi 
intervalli di tempo. Dopo opportune ipotesi semplificative deriviamo lo stimatore a 
massima verosimiglianza della temporizzazione risolvendo un problema di massimiz- 
zazione vincolata. II nostro stimatore a massima verosimiglianza ha una complessita 
proporzionale al numero di campioni di segnale trattati. Per mezzo di simulazioni le 
prestazioni dello stimatore a massima verosimiglianza vengono poi confrontate con 
quelle fornite da altri stimatori di temporizzazione proposti in letteratura. Lo stimatore 
proposto fornisce buone prestazioni indipendentemente dal tipo di canale considerate. 
In piu e poco sensibile alle variazioni dello spettro Doppler del canale e fornisce buone 
prestazioni sia in condizioni di fading veloce che di fading lento. 

In modo da aumentare la capacita del sistema, proponiamo successivamente un 
semplice ricevitore iterativo multiutente a cancellazione parallela di interferenza con 
equalizzazione retroazionata per sistemi CDMA. Nel nostro ricevitore, una stima del 
segnale interferente (interferenza da accesso multiplo e interferenxa inter simbolica) e 
ottenuto pesando le decisioni hard prodotte da decodificatori di Viterbi. 

La stima dell' interferenza e poi sottratta dal segnale ricevuto in modo da miglio- 
rare la decodifica nell'iterazione successiva. Utilizzando un'analisi asintotica delle 
prestazioni per sistemi CDMA con spreading casuale, ottimizziamo i pesi della retro- 
azione ad ogni iterazione. Poi consideriamo due effetti (tra loro correlati) che limi- 
tano le prestazioni del ricevitore: la polarizzazione dell' interferenza residua e Teffetto 
ping-pong. Mostriamo che Talgoritmo proposto puo essere migliorato compensando 
la polarizzazione nel calcolo dei pesi della retroazione e proponiamo una modifica 
deiralgoritmo base in modo da evitare effetti ping-pong. Tale modifica permette in- 
oltre di ottenere una piu alta capacita del sistema e una piu veloce convergenza alle 
prestazioni del singolo utente. L'algoritmo proposto e stato convalidato, tramite simu- 
lazioni, in un ambiente compatibile con le specifiche UMTS -TDD, includendo canali 
selettivi in frequenza, fading, utenti asincroni e controllo di potenza. Per tali appli- 
cazioni i decodificatori SISO (implementati ad esempio con Talgoritmo BCJR) non 
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sono dunque necessari per ottenere efficienze spettrali elevate, e, nella maggior parte 
dei casi, semplici decodificatori di Viterbi sono sufficienti. 

In questa tesi presentiamo poi algoritmi a bassa complesssita per Tinterfaccia tra 
il livello fisico e l'hardware adatti ad essere implementati su terminali Software Radio. 
Gli algoritmi proposti rendono la frequenza di campionamento indipendente dalla fre- 
quenza di simbolo e utilizzano la tecnica del campionamento a frequenza intermedia 
per la conversione analogico-digitale e digitale-analogico. Come applicazione parti- 
colare valutiamo l'effetto della distorsione introdotta da tali algoritmi sulla stima di 
canale, in un ambiente UMTS-TDD. 

Descriviamo infine l'architettura della piattaforma Software Radio sviluppata ad 
Eurecom e capace di fornire comunicazioni in tempo reale. La piattaforma e attual- 
mente utilizzata per Timplementazione di protocolli e algortimi di trattamento del seg- 
nale di terza generazione come UMTS e CDMA2000. Viene fornita una descrizione 
deirhardware con esempi di implementazione software sia su DSP ad alte prestazioni 
che su microprocessori standard. 



ii 



Advanced Signal Processing Algorithms for 3rd Generation Wireless 

Mobile Systems 

Alessandro Nordio 

Abstract 

Motivated by the worldwide activity around third generation mobile communication 
systems, this thesis addresses the problem of improving the physical layer reliability 
of W-CDMA third generation transceivers by designing simple and low complexity 
algorithms suited to be implemented on a real system. The work was developed in a 
framework of a joint project between Eurecom and EPFL whose objective is to design a 
real-time Software Defined Radio platform able to validate advanced signal processing 
algorithms. 

The research focuses on several aspects of the physical layer such as synchroniza- 
tion, multi-user detection, and design of the front-end signal processing. 

About the synchronization we consider initial timing acquisition in DS-CDMA 
when propagation is affected by multipath and fading and where the Base Station 
broadcasts a synchronization pilot signal in the form of bursts of modulated chips 
transmitted periodically and separated by long silent intervals. Subject to certain sim- 
plifying assumptions we derive the Maximum Likelihood (ML) estimator by solving a 
constrained maximization problem. Our ML timing estimator has constant complexity 
per observation sample. The relation to other estimation methods is addressed, and 
performance comparisons are provided by simulation. The proposed estimator yields 
good performance independently of the multipath-intensity profile of the channel, pro- 
vided that the delay spread is not larger than a given maximum spread. Moreover, our 
estimator is fairly robust to the mismatch in the fading Doppler spectrum and provides 
good performance for both fast and slow fading. 

In order to increase the system capacity we propose a low complexity multiuser 
joint Parallel Interference Cancellation decoder and Turbo Decision Feedback Equal- 
izer for CDMA. In our scheme, an estimate of the interference signal (both MAI and 
ISI) is formed by weighting the hard decisions produced by conventional (i.e., hard- 
output) Viterbi decoders. The estimated interference is subtracted from the received 
signal in order to improve decoding in the next iteration. By using asymptotic per- 
formance analysis of random-spreading CDMA, we optimize the feedback weights 
at each iteration. Then, we consider two (mutually related) performance limitation 
factors: the bias of residual interference and the ping-pong effect. We show that the 
performance of the proposed algorithm can be improved by compensating for the bias 
in the weight calculation, and we propose a modification of the basic PIC algorithm, 
which prevents the ping-pong effect and allows higher channel load and/or faster con- 
vergence to the single-user performance. The proposed algorithm is validated through 
computer simulation in an environment fully compliant with the specifications of the 
time-division duplex mode of 3rd generation systems, contemplating a combination 
of TDMA and CDMA and including frequency-selective fading channels, user asyn- 
chronism, and power control. For such application, Soft-Input Soft-Output decoders 
(e.g.implemented by the forward-backward BCJR algorithm) are not needed to attain 
very high spectral efficiency, and simple conventional Viterbi decoding suffices for 
most practical settings. 
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Moreover we present low-complexity algorithms for transmitter and receiver front- 
end suited to the implementation of Software Defined radio terminals. The proposed 
algorithms make the processing sampling frequency independent of the symbol rate 
of the digitally modulated signal and use the "IF-sampling" technique for D/A and 
A/D conversion. As a case-study, we consider a training-based joint multiuser channel 
estimation and we show that our front-end algorithms work nicely when coupled with 
an efficient FFT-based joint channel estimator. 

Eventually we describe a software radio architecture developed for providing real- 
time wide-band radio communication capabilities in a form attractive for advanced 3G 
system research. It is currently being used to implement signaling methods and proto- 
col similar, but not limited to, evolving 3G radio standards (e.g. UMTS, CDMA2000). 
An overview of the hardware system is provided along with example software imple- 
mentation on both high performance DSP systems and conventional microprocessors. 
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Introduction 



1.1 Historical perspective 

The first revolution in mobile telephony can be identified with the migration from ana- 
log systems (so called first generation (1G), such as TACS and AMPS [1]) to digital 
systems (second generation (2G), such as GSM, IS54, IS95 [1]). The existing second 
generation (2G) digital mobile standards offer bit rates up to lOkbit/s and have expe- 
rienced phenomenal growth since their launch. Now they are rapidly evolving towards 
the so called 2.5G systems including, for example, the GPRS (General Packed Radio 
Service) which offers air-interface transfer rates up to 115kbit/s. GPRS is a standard 
oriented to packed data communications developed by the European Telecommunica- 
tions Standards Institute (ETSI) and represents an essential step towards third genera- 
tion (3G) systems such as the Universal Mobile Telecommunication System (UMTS) 
currently under development and soon available on the market. UMTS will provide 
both traditional telecommunication services and new Internet-based services over the 
same network with higher bit rates (from 384kbit/s up to 2 Mbit/s) and will use radio 
frequencies around 2 GHz. Two main modes are currently considered for the UMTS 
standard: Time Division Duplex (TDD) and Frequency Division Duplex (FDD). The 
FDD mode is intended for applications in public macro and micro cell environments 
with wide area coverage and typical data transfer rates of up to 384 kbit/s. The TDD 
mode provides advantages for micro and pico cells and is especially suitable for en- 
vironments with high traffic density (e.g. business areas, city centers, airports) and 
indoor coverage where large data rates are required. 

In order to provide this high performance, 3G systems need an improved physi- 
cal layer. A huge amount of work has already been dedicated to the design of third 
generation transceivers, including for example the RAKE receiver, Multi-User Detec- 
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tion (MUD), multiple antennas, joint adaptive spatio-temporal signal processing and 
modern coding techniques such as "Turbo Codes" at both transmitter and receiver. 
Nevertheless, most algorithms are limited to a theoretical stage (a-part from few cases 
such as the RAKE receiver and Turbo Codes) and their practical feasibility remains a 
major issue. 

At present, most algorithms are tested only through computer simulation which 
often introduces a certain simplifications to model a real system. On the other hand, 
studies done by the main telecommunications companies about the feasibility and on 
software implementation of these solutions are not disclosed to the open literature. 

In this context in 1999 the Eurecom Institute and EPFL (Ecole Poli technique 
Federate de Lausanne) started a joint project whose objective is to study and implement 
a real-time Software Defined Radio (SDR) platform able to provide measurements in 
real transmission conditions and aiming at validating advanced mobile communication 
algorithms. 

Motivated by this project, the work presented in this thesis focuses on the analysis 
and design of a Wideband Code Division Multiple Access (W-CDMA) SDR platform, 
an architecture characterized by relatively few integrated circuits or hardware compo- 
nents except for some high performance Digital Signal Processors (DSPs). The plat- 
form has been initially targeted for the UMTS -TDD standard although other standards 
could be implemented. In other words, the UMTD-TDD standard has been consid- 
ered as a reference for the research developed in these years for the platform project. 
The aim was to improve the solutions proposed by the norm with more efficient and 
reliable algorithms, which concern many aspects of the transceiver physical layer sig- 
nal processing chain. In particular, the work of this thesis focuses on the following 
aspects: 

• Design of the transceiver front-end signal processing; 

• Synchronization; 

• Multi-user detection; 

Some considerations that motivated the work on these items are shortly discussed in 
the rest of this chapter. 

1.2 Front-end signal processing and the Software Radio con- 
cept 

In the future the proliferation of wireless mobile communications systems will re- 
quire radio terminals to be able to reconfigure themselves to handle several different 
standards. These terminals represent a really attractive solution to provide universal 
connectivity to the users and must be designed by using Software Defined Radio tech- 
niques. The basic idea of SDR is that some parameters characteristic of the transmis- 
sion system are chosen and programmed in software instead of using analog or digital 
dedicated hardware. Software mainly means generality, flexibility, and robustness, 
properties that nowadays have extreme relevance in the mobile communication con- 
text. SDR is a very broad term involving several levels in the protocol stack (e.g., see 
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[2], [3], [4] and references therein). At the physical layer, SDR requires signal process- 
ing algorithms suited to implementation on a programmable CPU. The goal is to move 
the software toward the antenna, performing operations like channel selection [5, 6], 
up- and down-conversion [7], synchronization and detection in the all-digital domain, 
by using high performance Digital Signal Processors (DSPs). 

In a modern transceiver, for example, the digital baseband signal produced by the 
transmitter is sent to the A/D converter, up-converted to Radio Frequency (RF) and 
then sent to the antenna. At the receiver side the RF signal incoming from the antenna 
is downconverted to baseband and sampled before being processed. The up and down- 
conversion operations are performed in hardware, passing through an Intermediate 
Frequency (IF) usually ranging from tens to hundreds of MHz. The base-band signal 
processing, instead, is constrained to certain parameters (signal bandwidth, chip or 
symbol rate) imposed by the implemented standard. In order to meet these constraints 
the hardware is created "ad hoc", with the disadvantage of not being able to support 
other standards characterized, by different parameters. 

The advantage of Software Radio instead, is that some parameters characteristic 
of the transmission system, such as the IF, the modulation scheme, the IF to baseband 
conversion, the baseband processing and other features are chosen and programmed 
in software. This solution enables the same hardware architecture to support many 
standards providing a software interface between the baseband processing and the 
hardware constraints. In practice, since the Intermediate Frequency is a parameter 
dependent on the analog modulation devices, in a SDR transceiver the software gov- 
erning the system is supposed to process the transmitter baseband signal in order to 
match it to the IF. This processing mainly consists in applying re-sampling techniques 
at both transmitter and receiver ends. In general, resamplers should be designed in 
order to allow resampling by non rational factors. 

In fact, a part of the work presented in this thesis concerns the design of such re- 
sampling techniques. Many front-end algorithms are proposed and compared, with 
particular attention to the computational complexity. For the transmitter side we dis- 
cuss an efficient generation of an IF analog signal from the baseband digital signal, 
while, at the receiver side, we propose efficients methods to generate a digital base- 
band signal from the IF analog signal. 

The proposed resampling techniques for the front-end signal processing, however, 
have an impact on the performance of the physical layer signal processing functions 
as for example channel estimation and data detection. This impact mainly consists of 
a degradation of the receiver performance in terms, for example, of Signal to Interfer- 
ence plus Noise Ratio (SINR) and Bit Error Rate (BER). For this reason the proposed 
techniques have been chosen and validated in a realistic UMTS -TDD system by find- 
ing a trade-off between degradation and computational cost. 

1.3 Synchronization 

Signal processing requires, in general, perfect synchronization between the transmitter 
and the receiver clocks. The synchronization is acquired when the receiver is switched 
on and then tracked while transmitting and receiving data. In UMTS -TDD, in or- 
der to facilitate the acquisition of the synchronization, a particular downlink signal is 
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broadcasted by the Base Station (BS). In order to assure reliable communications, the 
synchronism must be reliably and rapidly acquired. A good synchronization algorithm 
has a dramatic impact on the synchronization time, especially when the radio channel 
introduces severe degradations to the signal. Classical systems use simple algorithms 
although, in general, when the channel is affected by multipath and fading they show 
quite poor performance. 

In this thesis we present a Maximum Likelihood (ML) approach to the synchro- 
nization problem when the BS signal is corrupted by the presence of a multipath fading 
channel. The derived ML detector provides reliable timing estimations and, at the same 
time, it jointly estimates the channel Multipath Intensity Profile (MIP). This approach 
leads to an algorithm providing better performance compared to classical synchroniza- 
tion systems. 

1.4 Multi user detection 

In a typical CDMA system several users access simultaneously the same radio channel. 
This situation is common in a CDMA systems like the UMTS-TDD where the received 
signal can contain the superposition of several user signals. In such systems the Multi- 
ple Access Interference (MAI), caused by the users accessing the radio channel at the 
same time depends on the correlation between their spreading sequences. The solution 
provided by the UMTS-TDD norm consist of assigning orthogonal sequences to the 
users in order to nullify the MAI and by simple matched filtering the received signal 
using filters matched to the spreading sequences of each user. This idea, however, 
finds its limitations in a non perfect synchronization between users which destroys 
the orthogonality of the sequences. In any case, even with good synchronization and 
tracking algorithms, asynchronism between users is introduced by the propagation 
delays and orthogonality is destroyed by the fading multipath channels. In order to 
reduce the MAI effect, Multi User Detection [8] techniques seem very promising, pro- 
viding outstanding performance, superior to the conventional Single User Detectors 
and allowing very high channel loads. Multi User Detection, traditionally regarded 
as an ensemble of techniques to detect uncoded data in a multiple access waveform 
channel, are the subject of intensive studies in the recent years (see for example [8] 
and references therein). The proposed solutions, however, often require computational 
complexities not affordable by the available technology so that their implementation in 
the communications systems under development is still not considered, at least in the 
first system release. On the other hand, Information Theory shows that much larger 
channel loads can be achieved provided that a non-linear multiuser joint detector and 
decoder is employed [9, 10]. This may range from the impractically complex optimal 
joint decoder to practically appealing successive interference cancellation approaches 
(stripping) [11, 12] where decoded and re-encoded users are subtracted from the re- 
ceived signal. 

In this thesis, with the aim of providing simple and practical solutions for the im- 
plementation on the SDR platform, a Chapter is dedicated to the development of a 
low-complexity Iterative Parallel Interference Cancellation Multiuser Detection algo- 
rithm, that, at the same time, performs joint channel decoding and Decision Feedback 
Equalization. This algorithm can be easily implemented and used for real time appli- 
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cations with the available technology. 

1.5 Thesis outline and contributions 

The thesis is organized as follows: Chapter 2 gives an overview on the UMTS-TDD 
standard main features introducing the general system model of the transmitter and the 
receiver along with part of the notation used through this work. 

Chapter 3 presents the derivation of the Joint ML timing estimator for the UMTS- 
TDD mode in the presence of multipath fading channels. Comparisons with the per- 
formance of other timing estimators are provided. The results shown in this Chapter 
have been published in the following papers: 

• G. Caire, P. A. Humblet, G. Montalbano, and A. Nordio,"Slot-timing maximum 
likelihood estimation with bursty pilot signals for DS-CDMA in multipath fad- 
ing channels ," Proc. IEEE Globecom 2000, S. Francisco, CA, November 2000. 

• G. Caire, P. A. Humblet, G. Montalbano, and A. Nordio,"Initial synchronization 
in DS-CDMA via bursty pilot signals," to appear on IEEE Trans, on Communi- 
cations, 2002 

Chapter 4 develops a simple and high performance multi user detector with joint 
decoding and decision feedback equalization. Related publications are 

• A. Nordio, G. Caire, and M.Hernandez, "Design and Performance of a Low- 
Complexity Iterative Multiuser Joint Decoder Based on Viterbi Decoding and 
Parallel Interference Cancellation," Proceedings of ICC 2002, New York, May 
2002. 

• A. Nordio, G. Caire, and M.Hernandez, "Low-Complexity Turbo Equalization 
and Multiuser Decoding for TDMA-CDMA" submitted to IEEE Transactions 
on Wireless Communications, February 2002. 

Chapter 5 addresses some low complexity algorithms for a SDR transceiver front- 
end. Results, show that the proposed techniques represents a valid solutions for imple- 
mentation. The work presented in this Chapter has been partially published in 

• G. Caire, P. A. Humblet, G. Montalbano, and A. Nordio,"Transmission and re- 
ception front-end architectures for software radio," Proc. ECWT 2000, Paris, 
France, September 2000. 

The architecture of the W-CDMA SDR platform developed at Eurecom Institute 
is presented in Chapter 6. A description of the hardware set-up and of the employed 
signal processing algorithms is also given. Parts of this Chapter are contained in the 
following publications 

• C. Bonnet, G. Caire, A. Enout, P. A. Humblet, G. Montalbano, A. Nordio, and D. 
Nussbaum,"A software radio platform for new generation of wireless commu- 
nication systems," Proc. 12th Tyrrhenian Intern. Workshop on Digital Comm.: 
Software radio technologies and services, Porto Ferraio, Italy, September 2000. 
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• C. Bonnet, G. Caire, A. Enout, P. A. Humblet, G. Montalbano, A. Nordio, and 
D. Nussbaum "A software radio testbed for UMTS TDD systems," Proc. of 1ST 
Mobile Communications Summit 2000, October 2000. 

• C. Bonnet, G. Caire, A. Enout, P. A. Humblet, G. Montalbano, A. Nordio, and 
D. Nussbaum,"Open Software Radio Platform for New Generations of Mobile 
Systems," Proc. of 3rd European DSP Education and Research Conf., Paris, 
France, September 2000. 

• C. Bonnet, G. Caire, A. Enout, P. A. Humblet, G. Montalbano, A. Nordio, D. 
Nussbaum, T.Hohne, R. Knopp, B. Rimoldi, "An open software-radio architec- 
ture supporting advanced 3G+ systems," Annales des telecommunications, tome 
56, n.5-6, Paris, France, Mai-Juin 2001. 

Finally in Chapter 7 the general conclusions of this thesis are pointed out. 



Chapter 2 



A system model based on the 
UMTS-TDD physical layer 



The Universal Mobile Communications System [13] consists of a set of standards and 
new technologies conceived a few years ago and it is still under development. The 
aim of this standard is to provide to the users high-quality wireless multimedia broad- 
band communications, by extending todays mobile technologies (like GSM, WAP, 
GPRS [1]). This choice was motivated by the increasing need of multimedia appli- 
cations usually requiring high data rates such as pictures, graphics, video communi- 
cations, video on demand, mobile Internet access and other wide-band information as 
well as voice and data. 

The development of UMTS has the support of many major telecommunications 
operators and manufacturers worldwide. Several organizations take care of the UMTS 
standardization as, for example, the 3 rd Generation Partnership Project (3GPP) which 
is a collaboration agreement that brings together a number of telecommunications stan- 
dards bodies like the European Telecommunications Standards Institute. One of the 
3GPP scopes is to produce globally applicable technical specifications and technical 
reports for 3 rd generation mobile systems, including a new radio interface known as 
UTRA (UMTS Terrestrial Radio Access) based on the W-CDMA technology. Two 
main modes are considered for the UTRA: the Frequency Division Duplex [14] mode 
and the Time Division Duplex [15] mode. Although the UMTS standard considers for 
these modes several levels of the protocol stack, the research developed in this thesis 
only focuses on the UMTS-TDD physical layer, considering, in particular, the standard 
proposed by 3GPP. 

A short overview of the UMTS-TDD physical layer is presented in the following. 
Later, inspired by its main characteristics, we describe the transmitter and receiver 



7 



8 



Chapter 2 A system model based on the UMTS-TDD physical layer 




Figure 2.1: General structure of a cellular wireless system 



models used in this thesis. The interested reader can find more details on the UMTS- 
TDD 3GPP standard in [13]. 

2.1 Overview on the UMTS-TDD physical layer 
2.1.1 The multiple access strategy 

UMTS is a cellular system schematically represented in Figure 2.1. In the Figure, a cell 
is the region where a Mobile Terminal (MT) can access the network through a certain 
Base Station. Usually communications between MTs and BS are bi-directional. In the 
literature the words uplink and downlink [1] refer to the communication links going 
from the MTs to the BS and from the BS to the MTs respectively. 

The BSs, however, are more complicate than the MTs because they organize the 
flows of information transmitted and received by the MTs according to rules, schedul- 
ing, and algorithms defined in the supported protocols. Communication protocols that 
drive these systems are organized in layers and, in particular, the physical layer or 
layer 1 generates the radio signals to be transmitted and applies suitable signal pro- 
cessing algorithms to the received signals in order to assure reliable communications. 
Like for other digital systems, the UMTS baseband processing is performed in the all 
digital domain while the analog devices only take care of up- and down-conversion op- 
erations, acting from baseband or an Intermediate Frequency to the Radio Frequency. 
The UMTS-TDD layer 1 acts from the bit level to the Digital to Analog (D/A) or Ana- 
log to Digital (A/D) converters. It includes signal processing algorithms as synchro- 
nization, channel estimation, data detection, channel encoding and decoding, modula- 
tion, demodulation, and it implements the multiple access technique used to regulate 
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Figure 2.2: The UMTS-TDD slot structure 



the access of the transceivers to the radio resources. 

The UMTS 3GPP standard proposes different multiple access techniques for the 
FDD and TDD modes. The standard for TDD physical layer [13] recommends a com- 
bination of Time Division Multiple Access (TDMA) and Direct Sequence Code Divi- 
sion Multiple Access (DS-CDMA) techniques [1]. 

On the TDMA side, time is organized as a concatenation of frames lasting Tf rame = 
10 ms where each frame is divided in 15 slots of length T s i ot « 666 ^s. The MTs 
and the BS are allowed to transmit only in certain slots of the frame, called uplink 
and downlink slots respectively. These are dynamically assigned according to certain 
requirements as for example the uplink and downlink data rates, quality of service, and 
cell-load. 

The slot allocation strategy includes the possibility for more than one MT to trans- 
mit at the same time and in the same slot, causing Multiple Access Interference (MAI) 
between MT signals at the BS receiver. By using CDMA techniques, interfering sig- 
nals are separated by means of spreading sequences. 



2.1.2 Spreading and scrambling 

In the UMTS-TDD mode, a spreading sequence is a vector of chips belonging to the 
QPSK [1] constellation. The spreading sequence length, L, is also called spreading 
factor. In the standard, the chip time, T c , is defined such that T s i ot = 2560 T c , that is, 
2560 chips take place in a slot, yielding a chip-rate R c = 1/T C = 3.84Mchips/s. A 
simplified slot structure is depicted in Figure 2.2. Spread data and control information 
occupy only two parts of the slot, in the Figure labeled "Data Field 1" and "Data Field 
2", respectively. The field "Midamble" is used for channel estimation purposes and 
the guard time "Guard" is a silent interval separating the slots. 

The standard employs Walsh-Hadamard [16] spreading sequences with spreading 
factors of powers of two and in the range 1 < L < 16. Moreover the system can 
assign to the users different spreading factors in order to provide different data rates. 
The Walsh-Hadamard spreading sequences of length L = 2*, for a positive integer /, 
are the rows of the L x L Walsh-Hadamard matrix generated by the following 
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recursive procedure: 



H 2 l 

with initial condition 



Hl — Hl 



(2.1) 



Hi = 1. (2.2) 



The cross-correlation matrix of the Walsh-Hadamard sequences is proportional to the 
identity matrix, 

R Hi = H£H L = L1 L . (2.3) 

Since the MAI is proportional to the correlation between spreading sequences, the ad- 
vantage of using Walsh-Hadamard sequences is to nullify the MAI at the receiver. 
We have to point out that this consideration is valid only in the case of perfectly 
synchronous users. When users are asynchronous and their relative offsets contain 
a fraction of T c , (2.3) no longer holds. Moreover, since there are only L orthogo- 
nal sequences of length L, this choice provides a maximum load of 1 user per chip. 
In a general case, if the multiple access interference is not negligible as for example 
when using non-orthogonal spreading sequences or with Walsh-Hadamard sequences 
and asynchronous users, the MAI can be reduced by employing Multi User Detection 
(MUD) techniques [8] for the moment still not considered by the norm. This issue is 
broadly discussed in Chapter 4. 

The MAI, however, is not the only source of interference at the receiver. In fact 
a complete model should take into account also the cell-to-cell interference (i.e. the 
interference caused by transceivers of an adjacent cell). In particular in the UMTS- 
TDD mode, MTs of adjacent cells, can be assigned to the same spreading sequences, 
causing cell-to-cell interference between data signals. In order to reduce this effect 
cells are characterized by scrambling sequences that multiply the transmitted chips. 
The scrambling sequences are initially known by the BS and then acquired by the 
MTs during the synchronization process. 



2.1.3 Synchronization 

The frame and slot timings depend on the BS clock which is used as time reference 
by the MTs. This means that every MT has to acquire the BS time reference before 
starting communication in the cell. This process of timing acquisition is called syn- 
chronization. 

In the UMTS-TDD mode synchronization is performed in two steps, called pri- 
mary and secondary synchronization respectively. These are facilitated by the BS 
which broadcasts to the MTs some synchronization signals located in a known posi- 
tion of the frame. For the primary synchronization the BS broadcasts in a downlink 
slot a burst made of 256 modulated chips, 

s ps = [ 5 ps[0] 5 ■ ■ ■ 5 5 ps[255]j 

separated by long silent intervals and repeated indefinitely at a known rate. The pri- 
mary synchronization sequence, s ps , is common to all the cells and is known by the 
MTs. When a MT is switched on, it first detects the presence and the timing of the 
primary synchronization. 
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Figure 2.3: Generation of training sequences for the UMTS-TDD "Midamble type 2" 

The detection of the frame timing is facilitated by the structure of the primary syn- 
chronization sequence, designed in order to have a very peaky aperiodic autocorrela- 
tion function. Moreover the sequence s ps is a hierarchical sequence made of two base 
sequences of length 16 chips, s pSl and s pS2 . The primary synchronization sequence is 
then given by 

s ps = s psi ® s pS2 (2.4) 
where (® denotes Kronecker product [17]) and where 

SpSl = [-1,-1,-1,-1,-1,-1,+1,+1,-1,+1,-1,+1,-1,+1,+1,-lj 
s pS2 = [-1,-1,-1, +1, +1,-1, +1, +1,-1,-1,-1, +1,-1, +1,-1,-1] 

This feature in the design of the sequence is useful to reduce the complexity of the 
synchronization algorithm as shown in Chapter 6. 

When the frame time reference of at least one BS is successfully acquired, the 
MT searches for some secondary synchronization signal broadcasted by the BS and 
carrying additional information that characterizes the BS (e.g. the slot timing and 
the scrambling sequence index). Once all these parameters have been identified, the 
MT can set up a connection with the BS. For more information on the UMTS-TDD 
synchronization procedures refer to [14, 15]. 

2.1.4 Midamble and channel estimation 

In wireless communication the transmitted signals are degraded by the radio channel 
that must be correctly estimated at the receiver in order to improve the reliability of 
the communications. In the literature a huge amount of work has been published on 
channel estimation, and many techniques have been proposed. 

The UMTS-TDD standard defines a training sequence based channel estimation, 
where a sequence of known modulated chips is sent by the transmitter in the field 
"Midamble" of the slot (see Figure 2.2). In a given slot the training sequence chips are 
built from a common periodic base sequence, 

•t= [ 5 t[0], 5t [l],...,5t[Ptr]-l]] T 

of period Ptr by adding a cyclic extension of length C ex t > providing a total midamble 
length of P tr + (7 ex t chips. In an uplink slot, for example, the MTs allowed to transmit 
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Table 2.1: UMTS-TDD midamble parameters 
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build their training sequences using a cyclic shift of st by an integer multiple of C 0 ff as 
shown in Figure 2.3. This technique that allows a joint estimation of all the channels 
is described in [18], [19], and with some modifications in [20], and references therein. 
Usually P t r is an integer multiple of C 0 q such that the ratio P tr / C 0 ff = U indicates the 
maximum number of channels that the system can jointly estimate; C 0 q also represents 
the maximum length of the channel impulse response that can be estimated. 

The TDD standard defines two midambles of sizes 512 and 256 chips which pa- 
rameters are summarized in the Table 2.1 and allowing the joint channel estimation 
of U = 8 and U = 3 users respectively. An example of midamble generation is 
given in Figure 2.3. More details on the channel estimation signal processing for the 
UMTS-TDD mode are given in Chapter 6. 



2.1.5 Pulse shaping filter 

The last item discussed in this short overview of the UMTS-TDD physical layer is the 
transmission filter which shapes the spectrum of the transmitted signal by providing 
a finite output bandwidth. The UMTS standard proposes, as pulse-shaping filter the 
root-raised cosine waveform ip(t) whose spectrum is given by 



VTc 



»(/) = { 



\ 



sin 



1 ~r 0 s < - ,, < 1 +r off 

l/l>i±^ 



(2.5) 



and where 

/oo 
V>(*)e- j27r/ *di 
-OO 

In the standard, the roll-off factor is set to r 0 ff = 0.22 so that the total bandwidth of 
the UMTS signal is given by 

Since the chip-rate is R c = 3.84Mchips/s this corresponds to the bandwidth = 
4.68 MHz. 



2.2 General system model 
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2.2 General system model 

Driven by the main characteristics of the UMTS-TDD 3GPP system described in the 
previous Sections we introduce in the following the block diagram of a TDMA-CDMA 
transceiver on which the work presented in this thesis is based. Along with the descrip- 
tion we introduce part of the notation used in the following Chapters. Because of the 
generality of this architecture the UMTS-TDD physical layer can be seen as a particu- 
lar case and will be used through this thesis as a reference for comparisons. 

2,2,1 Transmitter 

In Figure 2.4 the block diagram of a TDMA-CDMA transmitter for both MT and BS 
is depicted. Information bits, denoted in the following by the vector b, and provided 
by the upper layers of the protocol stack, are channel-encoded, interleaved, modulated, 
and spread. The channel encoder labeled "ENC" in the Figure can be any code with 
coding-rate R. Blocks of encoded bits are interleaved with the interleaver II and then 
modulated using a modulation alphabet A (BPSK or QPSK [16]). Symbols, denoted 
by the vector 

a= [o[0],o[l],...] 

have energy 

£ = (log 2 \A\)RE h 

where \A\ is the cardinality of the modulation alphabet A and where is the energy 
associated to each information bit. In this system we consider complex spreading 
sequences with spreading factor L, given by the complex column vectors 

s u = [s w [0],s w [l],...,s w [L-l]] 

where u is the sequence index and where the elements s u [l] belong to the QPSK con- 
stellation. 

When the BS communicates at the same time with U MTs its output signal is the 
sum of U contributions, each one characterized by different symbols and spreading 
sequences as shown in Figure 2.4(b). 

In both cases, after spreading, complex chips are sent to the block labeled "Digi- 
tal and Analog front-end" that provides a RF signal to send to the antenna. For most 
systems in this block the D/A converter, which separates the digital signal processing 
from the analog devices is located. The position of the D/A converter inside the "Dig- 
ital and Analog front-end" block depends on the platform. An in-depth discussion of 
the baseband to RF signal processing for the transmitter is given in Chapter 5. 

In the theoretical studies presented in Chapters 3 and 4, for the sake of simplicity, 
we only consider the baseband processing, i.e. the "Digital and Analog front-end" 
block in Figure 2.4 is neglected. In this case, assuming an analog pulse-shaping filter 
?p(t), the baseband signal sent to the antenna is given by 

L-l 

x {t) = V£^a[k] s[l]tl>(t - IT C - kT) (2.6) 
k 1=0 
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(b) BS transmitter 



Figure 2.4: General scheme of the BS and MT transmitters 



for a MT transmitter and by 

U L-l 

= E a« [k] s u [lty(t - IT C - kT) (2.7) 

u=i k ;=o 

for the BS transmitter. The symbol period is denoted by T s = LT C . 

Although this system model has been focused on data transmission (which may 
involve coding and spreading) it can also represent the transmission chain for syn- 
chronization and training sequence signals. In fact, these can be seen as special cases 
where synchronization sequence chips and the training sequence chips are considered 
as source symbols while channel encoding, interleaving, and spreading blocks do not 
perform any operation. 



2.2.2 Channel 

The transmitted signal is corrupted by the radio channel which is characterized by 
the environment where radio signals propagate. The propagation conditions are in 
general assumed time-varying owing to the movement of the MT with respect to the 
BS although other effects can be involved (objects moving, changes of the weather, 
etc). Complete channel models, taking into account all these effects, are usually very 
complicated and when applied in the study of advanced signal processing techniques 
often lead to non-tractable solutions. A commonly used and simple channel model of- 
ten used in the literature, is the Wide-Sense Stationary Uncorrelated Scattering model 
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Figure 2.5: Representation of a channel MIP 

(WSS-US [16]) with Rayleigh fading and multipath impulse response, represented by 

p-i 

h(t,r) = ^2c p (t)S(r-r p ) (2.8) 

where the time-varying complex channel impulse response h(t, r) is a function of the 
time t and of the delay t. In this model the echoes of the RF signal reflected by the 
scatterers are grouped into P independent subsets, called paths, characterized by the 
delays r p and a set of time-varying mutually independent complex path gains c p (t) 
modeled as random processes with distribution Ne(0, cr 2 ). Usually these channels are 
defined by a Multipath Intensity Profile (MIP) [1], given by the delays vector 

t = [t 0 ,ti,...,t p _i] (2.9) 

and by the relative path powers vector 

<7 2 = [og, <!?,.■■ (2.10) 

An example of MIP is given in Figure 2.5. The MIP of certain propagation condi- 
tions have been standardized [21] with the purpose of giving a common reference for 
testing signal processing algorithms. 

The solutions proposed in this work have been tested using some of the standard 
channels as well as with other MIP created "ad hoc" in order to show particular behav- 
iors of the algorithms. 

2.2.3 Receiver 

The considerations done for the transmitter front-end hold for the receiver whose gen- 
eral structure is depicted in Figure 2.6. An in-depth discussion on the RF to base-band 
signal processing represented in the Figure by the block "Analog and Digital Front 
End" is given in Chapter 5. 

For the rest of the thesis, for simplicity, we assume baseband transmission. We 
will denote by y(t) the continuous time baseband received signal and by the vector 

y= [y[0],y[l],...] T 

its discrete-time version, sampled at rate R s = N c /T c where N c is a suitable integer 
that satisfies the Nyquist criterion. 
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Figure 2.6: Scheme of a TDMA-CDMA receiver 



In the uplink case where ?7 MTs are transmitting at the same time the baseband 
signal at the BS receiver is given by 

U /'oo 

2/(0 = E / h u (t,t-r)x u (r)dr + u(t) (2.11) 

where 

• v(t) is the noise contribution, a complex Gaussian random process with zero 
mean and power spectral density 7q that takes into account the thermal noise 
and cell-to-cell interference (assumed Gaussian) 

• x u (r) is the continuous-time complex baseband signal produced by the trans- 
mitter of the 11 th MT given by ( 5.1) 

• h u (t, r) is the channel connecting the 11 th MT to the BS 

In a downlink slot instead the only channel involved is the one that links the BS 
antenna to the MT receiver, so that the received signal is 

= £ / h{t,t-T)x u {T)&r + v{t) (2.12) 

Since we restrict this system model to a single cell, considering a single BS, the 
cell-to-cell interference is assumed to be part of the total noise variance. 

The synchronization block, present only in the MT, estimates the timing of the BS 
signal and provides a time reference to the receiver. The synchronization process is 
active when the MT is switched on or when, for some reason, the MT looses the syn- 
chronization and must re-acquire the timing. Using the slot timing informations the 
channel estimation block identifies the midamble into the slot, processes the training 
sequence signal, and provides channel estimates to the detector/decoder. The deriva- 
tion of an efficient signal processing technique providing the channel estimates is given 
in Section 5.5. 



2.2 General system model 
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Eventually, the detector/decoder processes the signal contained in the slot data 
fields, estimates the information bits, and sends it to the upper layers. It performs 
operations like matched filtering, despreading, demodulation, deinterleaving and de- 
coding. 

In this system model we do not specify the structure of the detector/decoder, leav- 
ing it enough general to include the Single User Detector receiver proposed by the 
UMTS -TDD norm as well as the Multi User Detection and decoding techniques de- 
veloped in Chapter 4. 
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Chapter 3 



Maximum likelihood burst timing 
acquisition 



We consider the initial timing acquisition in DS-CDMA when propagation is affected 
by multipath and fading and where the base-station broadcasts a synchronization pilot 
signal in the form of bursts of modulated chips transmitted periodically and separated 
by long silent intervals. Subject to certain simplifying assumptions we derive the max- 
imum likelihood (ML) estimator by solving a constrained maximization problem. Our 
ML timing estimator has constant complexity per observation sample. The relation 
to other estimation methods is addressed, and performance comparisons are provided 
by simulations. The proposed estimator yields good performance independently of the 
multipath-intensity profile of the channel, provided that the delay spread is not larger 
than a given maximum spread. Moreover, our estimator is fairly robust against the 
mismatch in the fading Doppler spectrum and provides good performance for both 
fast and slow fading. 



3.1 Introduction and motivation 

In a wireless mobile system synchronization represents a major issue. As introduced in 
Section 2.1.3 a MT before starting communication with the BS must acquire its timing. 
When data transmission occurs with a slot structure, as for the UMTS-TDD mode, the 
basic time reference is the slot and frame timing. In third generation wireless com- 
munication systems [22] timing acquisition is facilitated by a pilot signal transmitted 
by the BS. In the UMTS standard this signal is called primary synchronization signal. 
In both FDD and TDD [14, 15] modes, the primary synchronization signal (already 



19 



20 



Chapter 3 Maximum likelihood burst timing acquisition 



described in Section 2.1.3) is bursty, i.e., it is non-zero only for a small fraction of 
time compared with the slot duration. The first task of the MT when it is switched 
on consists of detecting the presence and the timing of the primary synchronization 
signal. Motivated by this scheme, we consider the general problem of initial slot tim- 
ing acquisition in a DS-CDMA system with a bursty pilot signal in the presence of an 
unknown multipath fading channel. 

In [23] a maximum likelihood (ML) timing estimator is derived by assuming that 
the channel MIP is known at the receiver, and in [24, 25] the ML criterion is applied by 
modeling the channel as deterministic unknown and constant with time. Unfortunately 
the multipath intensity profile is not known before initial acquisition and, since the ini- 
tial synchronization phase may last several slots, the algorithms based on the constant 
channel assumption might perform poorly in the presence of time-varying fading. 

In this Chapter, we obtain a low-complexity slot timing estimator requiring a min- 
imum amount of prior knowledge and yielding good performance on a wide range 
of channel MIPs and fading Doppler bandwidths. Since both the channel MIP and the 
noise plus interference power spectral density are not known at the receiver, we formu- 
late a joint ML problem where all these parameters have to be estimated. We derive the 
ML estimator by solving a constrained maximization problem via the Karush-Kuhn- 
Tucker (KKT) conditions [26]. In order to obtain a tractable solution we make several 
simplifying assumptions. When these are not satisfied, our estimator is not exactly ML 
and may suffer from mismatch. The proposed estimator is compared with the simple 
estimator for flat fading [16], with the too optimistic estimator of [23], and with the es- 
timator of [24] for two typical channel MIPs under both the assumption of slow fading 
(constant channel) and fast fading. Computer simulations show that the proposed joint 
ML estimator exhibits good performance in all these conditions and it is fairly robust 
to mismatch of the assumed fading statistics. 



3.2 Signal model 

The continuous-time complex baseband received signal is given by 

y(t) = t(t,e) + u(t) (3.1) 

where v(t) represents noise plus interference, modeled as a zero-mean complex cir- 
cularly symmetric Gaussian process with power spectral density J 0 , and 9) is the 
received synchronization signal component, given by 

M-l 
m=0 

where 9 is the slot timing, x(t) is the bursty pilot waveform of duration T ps , T ^> 
T ps is the period of repetition of x(t), h(t,r) is the time-varying multipath channel 
impulse response and M is the number of transmitted pilot bursts. Because of the 
slot periodicity, 9 must be estimated modulo T. The channel is assumed to follow 
the WSS-US [16] model with Rayleigh fading and multipath impulse response given 
by (2.8) where the channel MIP is defined by the delays vector r = [tq, . . . , Tp-{\ 
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Figure 3.1: Bursty synchronization signal at the transmitter 

and by the path variances a = [(7q, . . . , cr^ ^]. Since T ps is much shorter than T, 
we assume that the channel coherence time T CO h [16] satisfies T ps < T CO h < T. 
This implies that the channel is almost constant during each burst, but changes 
independently from burst to burst. Since T ps is very short 1 this condition is referred 
to as fast fading, and holds approximately for ^ < < where is the fading 
Doppler bandwidth. Moreover, we assume that the delays r p are constant over the 
whole observation window of duration MT. In practical systems [19] the period T 
ranges from one slot period to one frame period and the multipath delays vary at a 
much lower rate than the slot rate. Thus this assumption is satisfied. Subject to the 
above assumptions, Q{t, 9) in (3.2) can be rewritten as 

M-1P-1 

C(* 5 «)=EE -0-T p -mT) (3.3) 

where c mp are complex zero-mean circularly-symmetric Gaussian mutually uncorre- 
cted random variables such that E [c mp c* J = <7pS mjn 5 Pjq . 
The pilot waveform is given by 

N-l 

x{t)=Y< s n'4>{t-nT c ) (3.4) 

n=0 

where T c is the chip duration, s n is a sequence of TV chips known at the receiver and 
ift(t) is the chip-shaping pulse band-limited to [— with ^- < W$ < 
In a digital receiver implementation the signal is low-pass filtered and sampled at a 
convenient rate W > W^p. Hence, baseband processing is performed in discrete time. 
We assume W = AT C /T C , where N c > 1 is the number of samples per chip. Let 
Q = WT denote the number of samples per pilot repetition period, and define the 
discrete-time observed signal as 

y= [y[0],...,y[MQ-l]] T 
After a straightforward derivation, it is possible to write y in the compact form 

y = Sc + u (3.5) 

l In the UMTS standard proposal [14, 15] the pilot waveform consists of a sequence of 256 chips 
convolved with a root-raised cosine chip-shaping pulse. Thus with good approximation we may say that 
Tp S « 270 chip periods, whereas T can range from one slot (2560 chips) to one frame period (15 slots). 
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Figure 3.2: Model of the bursty synchronization signal 

where u = [i>[0], . . . , v\MQ — 1]] T is the vector of interference plus noise samples, 

c = [cj, . . . , c^_ x ] , with c m = [c m0 , . . ■ , c m ( F _!)] , contains all MP channel 
path coefficients over the M periods, and where S is the MQ x MP matrix whose 
(mP + p)^ column, for m = 0, . . . , M — 1 and p = 0, . . . , P — 1, is given by 



0 1 ^ : 0,5 P [0],...,5 P [Q-1], 0 1 ^ 1 0 ] T (3.6) 

mQ (M-m-l)Q 



where 



= -= s n?P(i/W -0-T p - nT c ) (3.7) 

^™ n=0 

Notice that the dependence of y on 9 is hidden in the matrix S as shown in Figure 3.2. 
Since the columns of S are obtained by translating the same waveform, they have all 
the same square magnitude, equal to the energy £ of the pilot waveform. Then, without 
loss of generality we can include the term £ into the path variances cr^, and consider 
l^mpl 2 = 1 for all m,p and delay 9. Under our assumptions, y is a zero-mean complex 
circularly-symmetric Gaussian random vector with covariance matrix 

R y = SAS U + I 0 I M Q (3.8) 

where the superscript (-) H denotes the Hermitian transpose, I m indicates the m x m 
identity matrix and where A is the covariance matrix of c, given by A = 1^ <g> A c 
with A c = diag(ag, . . . , (7p_i)- 



3.3 Maximum likelihood problem formulation 

The log-likelihood function for the parameters vector 9 = (0, r, cr, J 0 ) is immediately 
obtained as [27] 

L{y\0) = -logdet(R y )-y u R- 1 y (3.9) 



3.3 Maximum likelihood problem formulation 
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Notice that the number of channel paths P is also unknown and it is implicitly con- 
tained in the parameters r and a. By applying the matrix inversion lemma [17] to 
(3.8) we get 



Ry 1 = jImq - ^SA^ (l M p + Y^S u SA 1 2^j 



-l 



A2S h 



(3.10) 



Now, we make the key assumption that the columns of S are mutually orthogonal. In 
particular, this is approximately verified if the delays r p are sufficiently spaced apart 
(more than T c ) and if the sequence of chips s n has a very peaky aperiodic autocorrela- 
tion function. In practice, the pilot sequence s n is quite long (N = 256 in UMTS [19]) 
and paths spaced by less than one chip interval are substantially treated as a single path 
(they are not resolvable [16]). Thus, our assumption is not very restrictive. If S has 
orthonormal columns, (3.10) reduces to 



(3.11) 



where 



, and where we define the path average signal-to-noise ratio (SNR) £ p = cr p /Io, for 
p = 0, . . . , P — 1. Subject to the same assumption, the determinant of R y is readily 
obtained as 



det(R y ) 



p=0 



M 



jMQ 



(3.12) 



By using (3.11) and (3.12) in (3.9) and by defining the total received signal energy 
E y = \y\ 2 and the vector of path average SNRs £ = (£ 0 , ■ ■ ■ , £p-i)> we obtain the 
log-likelihood function in the form 



Hy\e) 



p-i 



p-1 



-MQlog/o-M£log(l + £ P ), 

(3.13) 



where we define 



M-l 



(3.14) 



m=0 



The term X(0 + r p ) is obtained by summing for m = 0, . . . , M — 1 the squared 
magnitude of the output of a discrete time filter with impulse response matched to the 
delayed pilot waveform x(t — 9 — mT — r p ). Then, this is the output of a sort of 
square-law diversity combiner [16] collecting the signal energy over all the M pilot 
repetition periods, for a given guess of the timing 9 and multipath component with 
delay r p . We shall refer to these terms as the received path energies at delay 9 + r p . 
From the implementation point of view, X(9 + r p ) for all p = 0, . . . , P — 1 can be 
calculated by sampling at delays 9 + r p the output of the same filter matched to x(t). 
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3.4 Approximated ML estimation 

In order to find the ML estimate of 0, we should jointly estimate also the other un- 
known parameters t,£ and 7q. For simplicity of exposition, we first assume that the 
delays r are known. Then, for any given value of 0 we want to maximize the log- 
likelihood function with respect to £ and 7q, i.e., we want to solve the constrained 
maximization problem 

f maximize £(y|0,£,r, J 0 ) n ^ 

\ subject to £ > 0, 7 0 > 0 K } 

This can be solved by using the KKT method [26]. In this section we let X p = X(9 + 
r p ) for the sake of notation simplicity. 
We have the following: 



Proposition 1. The solution of the maximization problem (3.15) is given by 
"1 



-,X p - 1 



where [•]+ denotes positive part, and where we define the set of indexes B = {pG 
[0, P — 1] : X p > /?}, of cardinality |D| = D. Moreover, this solution exists and is 
unique for any set of received signal energies {X p : p = 0, . . . , P — 1}, for all Q > P 
and M > 1 and for all E y such that Y^p=o X p < E y < p Y^p=o X v 

Proof. See Appendix 3.A.I. 

Solution (3.16) has the following intuitive interpretation: /? acts as an adaptive 
threshold level. If X p > (3 then 9 + r p is a good candidate for being a true channel 
path. 

Let 7r be the permutation of P indexes sorting the X p § in non-decreasing order, 
i.e., such that 

^tt(O) > X ir(l) >'"> X n(P-l) 
Then, by substituting the solution of Proposition 1 into (3.13) we obtain the maximized 
log-likelihood function in the form 

Z(y\9, t) = max£(y|0) = -(Q - D) log E * - _ i og JJ Xw{p) 

(3.17) 

The exact ML joint estimator of 0 is obtained by further maximizing £(y|0, r) with 
respect to all possible r and 9. While maximization with respect to 9 involves search- 
ing in a one-dimensional real space, with complexity 0(Q), the maximization with 
respect to r requires a search over P-dimensional real spaces for all P = 1,2,..., 
which is highly impractical. 

In order to circumvent this hurdle, we restrict the search over a set of tentative delay 
vectors r selected according to some criteria. At the same time, the set of tentative 
r should enforce the condition that the delays are separated by more than T c , which 
guarantees that the columns of S are approximately orthonormal. The selection of the 
tentative r is based on the following proposition and heuristic considerations. 
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Proposition 2. For a given fixed 0, let {r p : p = 0, 1, . . .} be a sequence of delays 
such that the sequence of corresponding received energies {X p : p = 0, 1, . . .} is non- 
increasing. Then, the sequence of log-likelihood functions To, - - - , r p ) : p = 
0, 1, . . .} is non-decreasing. Moreover, if for some D > 1 the condition 

X D <f} D ± ^ ~^ Xi (3-18) 

holds, then 

™£(y|0,To,...,T p ) = &(y\0,To,...,T D -i) 

P>\) 

Proof. See Appendix 3.A.2. 

From Proposition 2 it follows that we can construct a set of "nested" tentative delay 
vectors of dimension P = 1, 2, ... in sequence, by appending more and more delays 
corresponding to decreasing received energies. The maximum of the log-likelihood 
function over the (infinite) delay vector sequence is found when, for some (finite) 
integer D the condition (3.18) holds. Furthermore, we observe that: 

1. The log-likelihood function increases with the received energy sum 

D-l 
p=0 

and decreases with the product 

D-l 

B =U x *w 

p=0 

For a given value of A, £(2/|0,t) is maximized by a "peaky" distribution of 
the path energies X p s. In fact, the product B is maximum for the uniform 
distribution = • • • = X^d-i) = A/D while it is small if the received 

energies are spread, i.e., if X n ^ >• X *{1) > • > ^tt(D-i) (notice that 
^tt(D-i) > P > 0> so B cannot be zero). 

2. In practice, systems are designed to handle multipath channels up to a given 

maximum delay spread (see e.g. [19]) (defined here as = max p r p — 
miripTp [16]). Then, for a given 0, the channel paths are expected to be con- 
tained in a time window [0, 9 + TJ, where is known a priori. 

Driven by the above two considerations, a meaningful choice of the sequence of de- 
lays {r p : p = 0, 1, . . .} is to choose for each p the delay r p E [0, T^] for which X p is 
maximum and r p has not appeared before in the sequence. Moreover, since the paths 
should be separated by at least one chip interval, we place around each chosen r p a 
forbidden region of size 2T C . The newly selected delays should not belong to the for- 
bidden regions of already selected delays. The resulting approximated ML algorithm 
is given as follows. 
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Figure 3.3: An example of the Joint Maximum Likelihood algorithm behaviour 



Initial timing estimation. For all 9 E [0, T] we compute the value L{y \9) given by 
the recursion: initialize p = 0, /?n = E y /Q, and the search interval § 0 = [0, T<j], then 

1. Select r p = arg max rG s p X(0 + r), and let X p = X(0 + r p ). 

2. If X p > fa let = Ey ~^l=_\ X \ let § p+ i = § p \ [r p - T c , r p + T c ], 2 then 
increment p by 1 and go to step 1. 



3. If X p < j3 p or if S p+ i is empty, compute 



p-i 



i(y\0) = -(Q-p)iog/3p-iog 1 [[x i 

i=0 

and exit the recursion. 

— 

Finally, the estimated slot timing is given by 9 = 9 m + To, where 

0 m = arg max H(y\9) 

0<E[O,T] 



(3.19) 



and where to is the estimated delay of the first path (also provided by the above algo- 
rithm) corresponding to 9 m . 

As far as implementation is concerned, some considerations are in order. 

• Up to now, we considered the parameter 9 as continuous. However, in a digital 
implementation, the search for the maximum of L{y\9) is done on a discrete set 
of values. The most computationally demanding operation is the computation of 
the matched filter output z mp (9) = s^ p y involved in the calculation of X(9 + 



2 The notation A \ *B denotes the complement of the set £ with respect to the set A. 
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r p ) (see (3.14)). In most practical applications it is sufficient to acquire the slot 
timing with an error of less than one chip. Therefore, a very fine discretization 
of 0 is not needed. In our numerical examples, we discretize 0 with step 1/W, 
so that we need just a matched filter operating at the signal sampling rate W, 
whose output is given by 

*[«1 = £vb>(0' + 0/WT 

3 

The receiver accumulates the squared matched filter outputs in a vector buffer 
v b = [w&[0], ■ ■ ■ > Vb[Q - l]j such that 

M-l 

v b \i]= £ \z[mQ + i}\ 2 

The search of the delays and the search for the maximum in (3.19) is performed 
over the discrete values {tk = k/W : k = 0, . . . , Q — 1} by processing the 
data buffer Assuming 9 uniformly distributed, the minimum possible Mean 
Square Error resulting from discretization is given by the variance of a uniformly 
distributed random variable over [— +^] , i.e., by 12 ^ 2 - This is referred 
to as "discretization error" in Section 3.5. 

• We have implicitly assumed that the pilot bursts fall approximately in the middle 
of the observation intervals [mT, (m + 1)T]. However, the initial timing refer- 
ence of the MT is arbitrary, and the pilot bursts may fall across the boundaries 
of the observation intervals. Since the slot timing is defined modulo T, in order 
to solve this problem it is sufficient to apply the estimation algorithm by treating 
the data buffer as a circular buffer. 

• The complexity of the proposed algorithm is linear in the observation size QM, 
as opposed to other timing algorithms based on Least-Squares (LS) [24] or 
subspace decomposition, which require matrix-vector multiplication or matrix 
eigen-analysis (see e.g. [28] and references therein). 

• As byproducts of initial timing estimation, the proposed method also provides 
estimates for the interference plus noise power spectral density J 0 , for the path 
SNRs £p, for the path delays r p and for the number of paths P. These can be 
used to speed-up some terminal setup procedures [15, 14]. For example, the 
delays corresponding to the largest path SNRs can be used to initialize the fine 
delay search of a rake receiver, and the knowledge of J 0 can be exploited to 
initialize the power control loop. 



3.5 Results 

We compare the performance of our estimator, hereafter denoted by joint ML (JML), 
with the following estimators: 
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Figure 3.4: RMSE with constant CHI vs. £/7o for M = 10 accumulated bursts 

• The least computationally intensive estimator that selects as estimated delay the 
position of the maximum element in the buffer in the sequel denoted as the 
"MAX" estimator. 

• The approximated least-squares estimator (denoted by "LS") which is an ap- 
proximation of the exact LS estimator that coincides with the ML estimator sub- 
ject to the assumption of constant channel. In fact, when the channel is constant 
during the observation interval the ML slot timing estimator is given by [24, 25] 

6 = argmaxy H S(S H S)- 1 S H y (3.20) 
9 

where y = ± £m=i y m , y m = [y[(m - 1)Q], y[mQ - 1]] , and where 
S is the Q x L convolution matrix whose I th column is given by 

8 l =[s l [0],...,s l [Q-l}] T 

where 

i N-l 

and L is the channel delay spread expressed in sample periods. From the im- 
plementation point of view the algorithm corresponds to filtering the received 
signal y with a filter matched to the pilot waveform, windowing the filtered sig- 
nal by the L x L matrix (S H S) _1 , taking the magnitude square of the output 

— H — 1 

signal and finding the maximum. Even though the matrix (S S)~ can be pre- 
computed this approach is extremely computational intensive since it requires 
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Figure 3.5: RMSE with fast fading CHI vs. £/7o for M = 10 accumulated bursts 

the computation of L scalar products per output sample (i.e. for each possible 
0). Therefore we consider the more practical approximated LS estimator based 
on the assumption S U S « I, also proposed in [24]. 

• The ML estimator with perfect knowledge of both the channel MIP and the 
interference power spectral density 7 0 (denoted by "ML/known MIP") proposed 
in [23]. Since all the parameters but 9 are known the log-likelihood function 
(3.13) for this estimator reduces to 

P_1 P 

^y\9) = ^-f-X(9 + r p ) (3.21) 

As performance measure we use the root-mean square error (RMSE) of the esti- 
mated slot timing normalized with respect to the chip interval (i.e., it is expressed in 
fraction of T c ) and the cumulative density function (CDF) F e {E) of the sum of the 

energies of the paths falling in the estimated window W e = 9 + T^j , normalized 
with respect to the total channel energy. The CDF gives the probability that the frac- 
tion of the total channel energy "captured" inside W e is below a given level E. Since 
in general the performance of the subsequent synchronization phase, involving BS 
identification, frame synchronization etc., depends on the fraction of captured chan- 
nel energy (no matter which algorithm is used), the CDF is able to characterize the 
performance of the whole synchronization procedure, independently of the particular 
algorithm employed in the second phase. 

In our examples the chip-shaping pulse is root-raised cosine with roll-off factor 
a = 0.22. As pilot sequence, we used a PN sequence of length N = 255 and the 
UMTS primary synchronization sequence of length N = 256 defined in [13]. The 
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Figure 3.6: RMSE with constant CH3 vs. £/7q for M = 10 accumulated bursts 



receiver sampling rate is W = 4/T c and both the pilot repetition interval and the slot 
duration are equal to T = 625 T c , corresponding to Q = 2500 samples (in reality Q 
may be much larger but we were limited by the simulation time). In our simulations 
we considered four typical channel MIPs described in table 3.1 and denoted by CHI, 
CH2, CH3, and CH4. CHI is specified as a standard channel test in [21] and it has 
been chosen because it satisfies our working assumptions. CH3 was created "ad hoc" 
in order to test the algorithm on a channel which does not satisfy our assumptions. In 
fact CH3 has non-separable paths. Eventually, CH2 and CH4 were created "ad hoc" 
in order to put in evidence the energy-capture performance and to test the outlier error 
performance of the estimator, i.e., the capability of locating 9 within a window. 

Figures 3.4-3.7 show the timing RMSE versus the pilot energy to interference 
plus noise ratio £/7o for channel MIPs CHI and CH3 under both constant and fast 
fading conditions. We refer to fast fading conditions when the channel coherence time 
is much smaller than the period of repetition of the pilot (but still larger than the pilot 
duration). Notice that the "JML" and the "ML/known MIP" estimators are mismatched 
for constant channel, while the "LS" estimator is mismatched for fast fading. The 
number of accumulated bursts is fixed to M = 10. 

The "LS" estimator performs very poorly in the presence of fast fading because is 
mismatched (see figures 3.5 and 3.7). Indeed in the presence of fast Rayleigh fading, 
as in our scenario, the amplitude of the averaged received signal y = ^ Ylim de- 
creases as the number of accumulated bursts M increases. On the contrary the "JML" 
estimator shows good performance also in the presence of a constant channel despite 
it is mismatched in this case. Moreover it always performs closely to the "ML/known 
MIP" estimator. 



3.6 Conclusions 
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Figure 3.7: RMSE with fast fading CH3 vs. £/7o for M = 10 accumulated bursts 

Extensive simulations show that increasing M reduces the £/7o needed to achieve 
a given RMSE, without affecting the relative behaviors of the estimators. 

Figures 3.8-3.11 show the CDF F e (E) of the fraction of channel energy captured 
in the window W e for channels CH2 and CH4. In these figures the CDF of the "MAX" 
estimator shows abrupt transitions each corresponding to the loss of the channel energy 
associated with a path. For instance, in the case of channel MIP CH2, the transition 
occurs at E = — 3dB that corresponds to the loss of half the channel energy. This 
effect is due to a false lock of the "MAX" estimator on the the second path of the 
channel. 

In the presence of fast fading the "JML" estimator generally outperforms all other 
estimators but the "ML/known MIP" that is based on non-realistic assumptions. With 
slow fading our estimator performs better than the "MAX" and the "LS" estimators 
for high £/7o (i.e. greater than 10 dB when M = 10 as evidenced by figures 3.4 and 
3.6) even though it is mismatched. 

Finally figures 3.12 and 3.13 compare the performance of the "JML" estimator in 
terms of RMSE versus £/J 0 for fast fading channel CHI and CH3 when using the 
UMTS and the PN sequence, for M = 10 accumulations. Performance are worse with 
the UMTS sequence due to its worse acyclic autocorrelation properties. 



3.6 Conclusions 

Motivated by the initial BS acquisition procedure of UMTS, we considered the prob- 
lem of slot timing estimation based on bursty pilot signals. Subject to some simplifying 
working assumptions, we derived a low-complexity algorithm based on joint ML es- 
timation of the slot timing, of the multipath intensity profile and of the interference 
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Figure 3.8: CDF with constant CH2 for £/7 0 = 4dB and M = 10 accumulated bursts 

plus noise power spectral density and we solved the likelihood function maximization 
problem by using the KKT conditions. Comparisons with other algorithms such as 
the simple peak detection, the ideal algorithm that exploits perfect knowledge of the 
channel multipath intensity profile and an approximated least-squares are provided. 
The proposed algorithm shows good performance in all the considered conditions, es- 
pecially for very sparse channels and under fast fading conditions. Moreover it shows 
significant robustness to mismatch namely to the presence of constant channels and 
non-ideal acyclic autocorrelation properties of the pilot sequence. 



3.A Appendix 



3.A.1 Proof of Proposition 1 

We restate the maximization problem (3.15) more in general as follows. For arbitrary 
P, let {X p : p = 0, . . . , P — 1} be a non-increasing non-negative sequence, let Q > P 

and M > 1 be arbitrary integers and let Ylp=o X p < E v < p Y,p=o X v We want to 
maximize the function 



P-l 

n 1 + 



-X p — Ey 



P-l 

-MQlog/o-Mlog "[I X p 

p=0 



subject to £ > 0 and J 0 > 0. 

The KKT necessary conditions for a (local) maximum point of the objective func- 
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Figure 3.9: CDF with fast fading CH2 for £/J 0 = 4dB and M = 10 accumulated 
bursts 



tion in the constraint set yield the system of inequalities 

9 /(£,/o)<0, for i = 0,...,P-l 



56 

^/(^/o)<0 



(3.22) 



where strict equality must hold if ^ > 0 or 7 0 > 0. 

By solving the inequality for ^ we obtain & > JQ/(MJo) — 1. Since the RHS of 
the inequality for & is negative if & = 0 and JQ/(MJ 0 ) < 1, we conclude that the 
solution is given by 



6 = 



- 1 



(3.23) 



where we let /? = M7q. By substituting this into the inequality for Jo, we obtain 



p = 



Q-D 



(3.24) 



where D < P is the minimum integer for which < /?. Notice that since E y > 
T,p=o x P > then P > 0 thus > 0 and (3.23), (3.24) satisfy the KKT necessary 
conditions. 

Next, we prove that this solution exists and is unique. Finally, we prove that this is 
actually the global maximizer for our problem. 
Existence. For d = 1, . . . , P 9 let 

77! 1 y 



Pd = 



Q-d 



(3.25) 
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and define the set D d = {X p : X p > /3 d }. The solution (3.23), (3.24) exists if the 
equation |D^| — d has a solution for some D E {1, . . . , P}. 

By definition of we have that Xq > /?]_. If for 1 < p < P — 1 the condition 
X p > /3 P is never verified, then the solution is obviously D = 1. Otherwise, there 
must exist some 1 < p < P — 1 such that X p > /3 p . Let cf be the maximum of such 
indexes. We can write 



7j y ~ £^p=o y ^p 




which implies Xd > Pd+i- Since by construction d 
Xd+i < ftd+i> which implies that D = d + 1 is a solution. 

Uniqueness. Suppose that there exist 1 < D < D f < P such that \T) D \ = D and 
=ZT,i.e., 

X p > /3 D for 0 < p < D - 1 
X p <j3 D for D < p < P - 1 



{ 
{ 



and 

X p > for 0 < i < £>' — 1 
X p < for D' < i < P — 1 



(3.27) 
(3.28) 



By following the same steps of (3.26) we can show that 

j3 d > X d 0 d+1 > X d (3.29) 

Then, starting from (3.27) we can write the chain of inequalities /3d > Xd Pd+i > 
X D > X D +i =^ Pd+2 > -*Cd+i > -*Cd+2 > X D t_i, which contradicts 

(3.28). 

We conclude that the solution must be unique. 

Global maximum. The function 7q) is not concave in However, it is 

continuous on the (convex) constraint and /(£, Jo) — > — oo if any of its variables 
grows without bound. By continuity, the global maximum of /(£, J 0 ) is finite and the 
maximizer must have all finite components, therefore it must satisfy the necessary 
KKT conditions. Since the solution (3.23), (3.24) exists and is unique, it must be the 
global maximizer. 



3.A.2 Proof of Proposition 2 

Let {X p : p = 0, 1, . . .} denote a sequence of non-increasing path energies. If for some 
D the condition (3.18) holds, then (3d > X p for all p > D and, by Proposition 1, the 
function £(y|0, fo? ■ ■ ■ ? fp 5 to, ■ ■ ■ ? Tp, Jo) for all p > D is maximized with respect to 
the fc's and J 0 by letting (3 = (3 D = (E y - Y^?^ X i)/(Q ~ D )> 7 o = 0/M and & = 
[Xi/f3 — 1]+. Since /?d is determined only by the Xj's for i = 0, . . . , D — 1, by adding 
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Figure 3.10: CDF with constant CH4 for £/7 0 = 4dB and M = 10 accumulated 
bursts 
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Figure 3.12: RMSE with fast fading CHI vs. £/J 0 for UMTS and PN sequences and 
M = 10 accumulated bursts 




Figure 3.13: RMSE with fast fading CH3 vs. £/J 0 for UMTS and PN sequences and 
M = 10 accumulated bursts 
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Table 3.1: Channel MIPs 



MIP 


Delays (T c ) 


Variances (dB) 


CHI 


(0,1,2,3) 


(0, -3,-6, -9) 


CH2 


(0,7.5) 


(0,0) 


CH3 


(0,1.12,1.34, 4.85,5.35,5.62) 


(-3,0, -3,-5,-7, -8) 


CH4 


(0,7.5, 15) 


(0,0,0) 



more delays to the sequence To, ... , td-i with received energy X p < (3d is not going 
to affect neither the value of (3d nor the value of the maximum £(y|0, Ti, . . . , r p ). 

Therefore, the only situation that must be taken into account is when, for given p, 
(3 p < X p -i and we add a delay r p with received energy X p > (3 p . From the proof of 
Proposition 1 we have that 

X p > /3p=> X p > /3 p +i 

i.e., £(y|£o 5 ■ ■ ■ •> To, ■ ■ ■ •> t v -> ^o) is actually maximized by all & > 0, implying that 
T o> ■ ■ ■ ■> t p-i) T o> ■ ■ ■ ? t p ). Proposition 2 is proved by showing that 

£(y|0,T O ,...,T p _i) < H{y\6,TQ,...,T p ) 

By using the fact that 

and by letting Q —p = kwe have 

C(»|0, T 0 , ... , T p ) - £(y|0, T 0 , ... , Tp_l) = 

P P-l 

= -(* - 1) log &+i - log n ^ + a log & + log n x< 

*[k-l\X p k)_ 

=(fc-i)io g ; 7 . +*bg^ 

>o 

since the arguments of both logarithms in the third line above are larger than 1. This 
concludes the proof. 



= -Hog 



- 1 V 
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Chapter 4 



Iterative multiuser detection 



We propose a low complexity multiuser joint Parallel Interference Cancellation (PIC) 
decoder and Turbo Decision Feedback Equalizer for CDMA. In our scheme, an es- 
timate of the interference signal (both MAI and I SI) is formed by weighting the hard 
decisions produced by conventional (i.e., hard-output) Viterbi decoders. The estimated 
interference is subtracted from the received signal in order to improve decoding in 
the next iteration. By using asymptotic performance analysis of random- spreading 
CDMA, we optimize the feedback weights at each iteration. Then, we consider two 
(mutually related) performance limitation factors: the bias of residual interference 
and the ping-pong effect. We show that the performance of the proposed algorithm 
can be improved by compensating for the bias in the weight calculation, and we pro- 
pose a modification of the basic PIC algorithm, which prevents the ping-pong effect 
and allows higher channel load and/or faster convergence to the single-user perfor- 
mance. The proposed algorithm is validated through computer simulation in an envi- 
ronment fully compliant with the specifications of the time-division duplex mode of 3rd 
generation systems, contemplating a combination ofTDMA and CDMA and including 
frequency-selective fading channels, user asynchronism, and power control. The main 
conclusion of this work is that, for such application, Soft-Input Soft-Output (SISO) de- 
coders (e.g., implemented by the forward-backward BCJR algorithm) are not needed to 
attain very high spectral efficiency, and simple conventional Viterbi decoding suffices 
for most practical settings. 



4.1 Introduction and motivation 



In both FDD and TDD modes of the UMTS standard, the basic receiver scheme 
contemplates the use of conventional Single-User Matched Filtering (SUMF). Since 
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Multiple- Access Interference is treated as additional background noise , powerful 
and high-complexity channel coding such as 256-states convolutional codes and turbo 
codes [29] are envisaged in order to attain low Bit Error Rates (BER) at low decoder 
input signal-to-interference plus noise ratio (SINR). In any case, channel loads larger 
than 1 user/chip are practically very difficult if not impossible to attain by the SUMF 
front-end and single-user decoding [9]. 

On the other hand, Information Theory shows that much larger channel loads can 
be achieved provided that a non-linear multiuser joint detector and decoder is em- 
ployed [9, 10]. This may range from the impractically complex optimal joint decoder 
to practically appealing successive interference cancellation approaches [11, 12]. 

In practice, successive interference cancellation must cope with decision errors, 
which prevent perfect cancellation of already decoded users. Then, several iterative 
schemes have been proposed, which limit the deleterious effect of decision errors by 
feeding back soft-estimates of the detected symbols (see for example [30, 31, 32]). 
These schemes require Soft-Input Soft-Output (SISO) decoders, usually implemented 
by the forward-backward BCJR algorithm [33]. However, such SISO decoders repre- 
sent a non-negligible factor in the complexity of whole receiver. 

In real CDMA applications for either TDD and FDD modes, the maximum achiev- 
able channel load is often limited by synchronization and channel estimation issues, 
rather than by the ultimate capability of the decoder itself [34]. Hence, it makes sense 
to investigate simpler joint detection and decoding schemes, which outperform the 
conventional linear SUMF, MMSE and decorrelator, and non-linear Parallel Interfer- 
ence Cancellation (PIC) or Serial Interference Cancellation (SIC) receivers [35], and 
nevertheless yield performance similar to the SISO-based schemes at lower decoding 
complexity. 

Driven by this consideration, this Chapter proposes a low complexity iterative 
multiuser receiver, where SISO decoders are replaced by simpler standard (i.e. hard- 
output) Viterbi decoders. The Viterbi hard decisions are weighted and fed back to the 
interference cancellation stage. By using large-system analysis of random CDMA [36] 
we optimize the feedback weights at each iteration such that the SINR at the decoders 
input of the next iteration is maximized. 

We address the problem of bias of residual interference [37, 38] and of the ping- 
pong effect [35, 39] and we show that the performance of the proposed algorithm can 
be improved by compensating for the bias in the weight calculation and by modifying 
the basic PIC algorithm in order to prevent the ping-pong effect. These modifications 
allow higher channel load and/or faster convergence to the single-user performance at 
almost no additional computational cost. 

We validate the proposed receiver algorithm in UMTS-TDD realistic scenarios, 
including asynchronous transmission, frequency selective fading channels and power 
control. In this regime, the proposed receiver performs very close to the single user 
(i.e., MAI-free) Matched Filter Bound (i.e., ISI-free) performance, even for large chan- 
nel load. 

The remainder of this Chapter is organized as follows: Section 4.2 gives a descrip- 
tion of the system model. Section 4.3 describes the large system asymptotic analysis. 



Notice that the only difference between the SUMF and the linear MMSE fl Iter is that the latter treats 
MAI as colored noise while the former treats it as white noise [8] . 
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Based on the latter, the feedback optimal weights are derived in Section 4.4 for syn- 
chronous and flat channels, taking into account the bias on the residual interference and 
ping-pong effect for large channel load. Section 4.5 deals with asynchronous transmis- 
sion and frequency selective fading channels. Finally in Section 4.6 conclusions are 
pointed out. 

4.2 System model 

We consider the uplink of a DS-CDMA system where U users send encoded infor- 
mation to a common receiver. The baseband transmission chain for the IX th user is 
depicted in Figure 2.4. Source bits are channel-encoded and organized in codewords 
of length TV bits. The codewords are then interleaved and modulated with transmitted 
energy per symbol £ w . For simplicity we assume that all users make use of convolu- 
tional coding and BPSK modulation. Let a u [n], the 71 th symbol generated by the 11 th 
user, be in the set {—1, +1}, and 

a u = [a u [0],a w [l], . . . ,a u [N - 1]] 

represent the code word of user u after interleaving. The symbols are then spread by 
the spreading sequence 

s u = [s w [0],5 u [l],...,5 w [L-l]] 

We assume all the users have the same spreading factor (number of chips per symbol) 
L and unitary energy sequences, that is, 

Vu = 1, . . . , U. If the chip pulse-shaping filter is ip(t) and is common to all the users, 
then the 11 th baseband continuous-time transmitted signal is 

N-l L-1 

x u (t) = a ^ s u[lW - IT C - nT) (4.1) 

n=0 1=0 

where T c is the chip period and T = LT C is the symbol period. 

We focus on UMTS -TDD, where users are synchronous at the slot level but asyn- 
chronous at the chip level. The channels are considered random but slowly varying 
and constant over one slot. Moreover, for the sake of simplicity, we assume that the 
convolutional codewords span a single slot. 

Thus, the 11 th user signal is sent through the channel h u (t) with impulse response 
given by 

Pu-l 

h u (t) = ]T) Cu, p 5(t-T UiP ) (4.2) 

p=0 

where P u is the number of resolvable paths, r U:P is the p ih path delay, and c U:P is the p ih 
path coefficient assumed a Gaussian complex circularly symmetric random variable 
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Figure 4.1: Transmission channel 



with distribution Ne(0,a 2 p) where = E[|c UjP | 2 ]. This model takes into account 
also users' asynchronous transmission, where the relative delays between users are 
included into the channel path delays. 
The signal at the receiver is given by 



/+oo 
X u {T)h u {t-T)dT + u{t) 
u=i -°° 
U N-l 

= ^2^2a u [n]g u {t-nT) + u{t) 



(4-3) 
(4.4) 



u—1 n—0 

where v(t) is the Gaussian noise and g u (t) is defined by 

L-l P u -1 



9u(t) = V&u ^2 Su ^ ^2 hu # " lT ° ~ Tu >v> 

1=0 p=0 



(4.5) 



and represents the overall channel impulse response for the 11 th user as depicted in 
Figure 4.1. 

The signal y(t) is sampled by the receiver at rate N c /T c , an integer multiple of the 
chip-rate, where N c is a suitable integer chosen in order to satisfy the Nyquist criterion. 



The discrete time version of g u {t), sampled at rate 7V C /T C , may have an infinite 
support but most of its energy is concentrated in a finite interval. This interval de- 
pends on the channel maximum delay spread, on the pulse-shaping filter ip(t), and on 
the spreading factor L. After a suitable truncation we can represent the discrete-time 
channel impulse response as a vector 

g u = [g u [0],...,g u [N g -l]\ T 

Vu = 1, . . . , U where N g is the maximum discrete-time channel length of all user 
channels. Subject to these assumptions, the received discrete-time baseband signal 
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Figure 4.2: Structure of the matrix G u 



can be written as follows 



U 



y = ^2 G uQ>u + v = Ga + v 

u=l 



where 



(4.6) 



• V — [y[0], ■ ■ ■ — 1)LN C + 7V ff ]J is the vector of the received signal 
samples. 

• G is a ((TV - 1)LJV C + 7V ff ) x ?77V matrix given by G = [Gi, . . . , G^] where 
the ((TV — 1)LN C + N g ) x N matrices G u contain, in each column, a shift of 
the vector g u as shown in Figure 4.2. 

• a = ^c*i , . . . , a J- J is the concatenation of the U user codewords. 

• u = [^[0]) ■ ■ ■ j ^[(^ ~~ ])LN C + 7V ff ]J contains the noise samples where the 
i/[i]'s are i.i.d. circularly symmetric complex Gaussian random variables with 
distribution u ~ Ne(0, Jol). 

For the sake of simplicity, the receiver front-end is constrained to be the conven- 
tional Single User Matched Filter (SUMF) [8] given by the matrix G H . 2 The matched 
filter performs at the same time pulse-shape matched filtering, channel matched filter- 
ing, and despreading. The output of the bank of SUMFs is given by [8] 



z = Re{G H y} 

= Re{G u Ga + G u u} 
= Ra + v 



(4.7) 



where 



2 Some works as for example [31] and [32] consider linear MMSE front-end. 
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Z — •••f Z U 

concatenation of 



the 



and z u = 
U matched 



-iT 

z u [0], . . . , z u [N — 1] for u = 1, . . . , U is the 
filter outputs. 



• R is the UN x UN cross-correlation matrix given by R = Re {G U G} where 
(•) H denotes the Hermitian operator. 

• the noise term, v = jv^, . . . , vjj is given by v = Re {G H i/} and has distri- 
bution v ~ N(0, § Jt). The vector 

v u = Re {G^u} = [v u [0], . . .,v u [N - 1]] T 

represents the additive colored noise contribution after the 11 th matched filtering. 
In particular, the SUMF output for the 71 th symbol of the IX th user is given by 

UN-l 

z u [n] = z[uN + n ]= ^ R[uN + n, k]a[k] + v[UN + n] (4.8) 



The matrix R is Hermitian and contains U 2 blocks of size N x N each. In the follow- 
ing, Ru^ indicates the (u, v)^ block given by 



Ru,v — Re {Gy]G v } 



(4.9) 



The [i, j] entry of R UjV is given by 



i* Uj ^ = i*[uA^ (4.10) 



From the structure of G u , depicted in Figure 4.2, and from (4.10) it is clear that the 
matrices R UjV are banded with 2D + 1 non-zero diagonals where 



D = 



3l 

LNr 



(4.11) 



Notice that since R u ,v[hj] depends only on the difference (j — i), R UjV is uniquely 
defined by the vector r UjV = [r Ujt ,[-D], . . . , r Ujt7 [0], . . . , r W)U [+D]J where 
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Figure 4.3: Scheme of the proposed turbo multiuser receiver 



Hence, using (4.7) and the above considerations it is possible to rewrite (4.8) as 

U-1N-1 

z u[n] = ^2^2 R ™A n i k ] a v[k] + v u [n] 

v=0 k=0 
U-l D 

v=0 d=-D 
U-l D 

= ^2 ^2 r «,v[ d K[ n + <*] + v u[n] 

v=0 d=-D 

D U-l D 

= \9u\ 2a u[n]+ ^2 r ^[d]a u [n + cf] + ^ ^ r UfV [d\a v [n + d\ +v u [n] 

d=-D v=0 d=-D 

d^O 



s > 



ISI 



MAI 



(4.12) 



where the first term is the useful symbol, the second term denotes the Inter Symbol 
Interference (ISI), the third term is the MAI, and the last term is the colored noise. 

The SUMF output feeds the proposed iterative multiuser detector-decoder depicted 
in Figure 4.3. After SUMF, the signal z passes through an Interference Cancellation 
(IC) stage that uses the estimates ai m ^[n] of the symbols a u [n] to remove ISI and 
MAI (the superscript ( m ) denotes ra 1 * 1 iteration). The hard decisions ai m ^[n] E {±1} 
are provided by a bank of Viterbi decoders, whose outputs are weighted by a factor 
/?£ m ^ E [0, 1]. Thus, the signal at the output of the IC stage in the m th iteration is given 



4»)[n] = z u [n] 



( 



D 



U-l D 



\ 



£ r u , u [d\p^aW[n + d\ + Y, £ r u , v [d\^a^[n + d\ 

. d=-D 

= \9u\ 2 au[n) + Ci m) [n]+v u [n] 



v=0 d=-D 



J 

(4.13) 
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where 

D 



(i m) [n] = E r «,«H (a„[n + d\ - /3^a^[n + d]) 

+ E E r vH (<M« + 4 " ^ m) 4 m) [« + <fl) 



(4.14) 



u=0 d=-£> 

is the residual MAI+ISI at iteration ra, and v u [n] ~ DST (o, f °^ u ^ ) • 

Equation (4.13) shows that the IC stage performs at the same time MAI suppres- 
sion and also non-causal Decision Feedback Equalization (DFE) by removing from 
the 71 th symbol the contributions due to the past and future symbols of the same user. 

At the first iteration, the initial estimated symbols are set to zero, a^[n] = 0, 
u = 0, 1, . . . , U — 1 and n = 0, . . . , N — 1 so that zi^fn] = z u [n]. In the case of 
perfect symbol estimates and (3^ = 1, equation (4.13) reduces to 

zi m) [n] = \g u \ 2 a u [n]+v u [n] (4.15) 

where the MAI and ISI are completely removed and the single-user Matched Filter 
Bound (MFB) performance for user u is attained. For later use, we define the nor- 

\Q I 2 

malized received instantaneous channel energy for user was 7 U = such that 

E[j u ] = 1, and the single-user MFB instantaneous SNR for user u as 

SNRr B = ^7. (4.16) 

Intuitively, the weighting factors should depend on the reliability of the esti- 
mated symbols [n] and be equal to 1 or to 0 in the case of completely reliable or 
completely unreliable symbol estimates, respectively. In the next Section we make this 
statement precise by studying the performance of an idealized synchronous system in 
the large-system limit. 



4.3 Large system asymptotic analysis 

Rigorous asymptotic analysis of this iterative scheme in the large-system limit (i.e., 
for TV, J7, L — >> oo with finite ratio U/L = a and U/N — > 0), for single-path chan- 
nels and synchronous users is proposed in [37]. This analysis is based on the general 
approach of density evolution over graphs: a standard analysis tool for evaluating the 
asymptotic performance of message-passing iterative decoders [40], and makes use of 
results from the theory of large random matrices developed in [36] for the analysis of 
linear receivers with randomly spread CDMA in the large-system limit. In [37] it is 
also shown that this analysis holds only if the symbol estimates ai m ^[n] are functions 
of the decoder extrinsic information [41]. The decoder extrinsic information is defined 
for a SISO decoder based on the sum-product algorithm [42], such as the BCJR algo- 
rithm, but it is not defined for an ML decoder such as the Viterbi algorithm. Hence, 
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we shall optimize the weights (3^ for a. fictitious receiver where the Viterbi decoders 
are replaced by BC JR decoders and the hard decisions [n] are obtained by one-bit 
quantization of the extrinsic likelihood ratios produced by the latter. 

We hasten to say that the fictitious receiver is used as a reference but has no practi- 
cal relevance, for the obvious reason that if BCJR decoders are used, then much more 
efficient soft-estimation of the interfering symbols as in [30, 31, 32] is possible. How- 
ever, as it will be clear from the rest of this section, the weight optimization based on 
the asymptotic analysis of the fictitious receiver allows us to derive a very simple ex- 
pression for the optimal weights, independent of the user sequences and their mutual 
correlations, and a very simple practical algorithm for calculating these weights. 

We make the following assumptions: 

• complex random spreading sequences formed by i.i.d. chips uniformly dis- 
tributed over the QPSK constellation; 

• users are slot and chip synchronous; 

• equal received expected energy per symbol £ for each user; 

• frequency-flat identically distributed propagation channels; 

• we let TV, [7, L — >> oo with U/L = a (a denotes the channel load); 

• the empirical distribution of the user received normalized powers, defined by [36] 

F f ) (^) = ^E 1 {7.<^} (4.17) 

converges weakly to a given non-random cumulative distribution function F 7 (z), 
for U — >> oo. 

With the above assumptions, for a user u randomly selected with uniform proba- 
bility in the user population, we have SNR^ 6 = j^Ju where j u is a random variable 
distributed (in the limit of large U) according to F 7 (z). 

Under these conditions, the SINR at the decoders input in the 771 th iteration for the 
fictitious receiver can be written as follows [37] 

op 

SINR£"> = ?/ m )SNR^ FB = t/™) -f 7u (4.18) 

^0 

where 

rf m) = i— ^ (4-19) 



1 + a^( m ) 



is the degradation factor paid by any user with respect to its single-user MFB perfor- 
mance, and where 

a _£(m) a (m) 2 (4> 20) 



/j( m ) = E 



is the average variance of the residual interfering symbols. In (4.20), a denotes a 
coded symbol of a generic user received at instantaneous power 7, and a( m ) denotes 
the hard decision relative to a at the decoder output. Expectation in (4.20) is with 
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respect to a ~ uniform over {±1}, a( m ) distributed as the hard decision about a at the 
decoder output, and 7 ~ F 1 {z). We refer to rf m ^ as the Multiuser Efficiency (ME) [8] 
at iteration m. Clearly, the single-user MFB performance is achieved for all users if 
lim m _ 00 ^"») = l. 

We can optimize the weighting factor /?( m ) as a function of the channel energy 7 
in order to minimize the expected interference variance /i( m ) at every iteration. There- 
fore, for each iteration m we seek the solution of the simple optimization problem 



min E 

0=0(7) 



j3a 



7 



(4.21) 



By expanding the above conditional expectation, we get 

1 + /3 2 -2/3E[aaW\j] 
l + j3 2 -2j3 [(+1)(1 - e ("%)) + (-l)e(m)( 7 )] 
1 + /3 2 -20(l-2e( m )( 7 )) 



(4.22) 



where e( m ) (7) is the Symbol Error Rate (SER) of a decoder with an input signal-to- 
noise ratio SINR( m_1 ) = j/" 1-1 )^. The solution of (4.21) is easily obtained as 



/?*(7) = l-2e^)( 7 ) 



(4.23) 



Assuming that the residual interference plus noise at iteration m is Gaussian 3 the SER 
is a known function /(•) of the decoder input SINR. We shall refer to /(•) as the SER 
characteristic of the user code (see Appendix 4.A). Thus, for user u at iteration m the 
optimal weighting factor is given by 

= F{lu) = 1-2/ (V™" 1 ^) (4-24) 
The resulting minimized average interference variance is given by 

^ = 4E [ 7 / (y m_i) f 7) (1 - / (y m_i) f 7)) 



(4.25) 



where expectation is with respect to 7 ~ F 7 (z). 

The behavior of the asymptotic system is described by the evolution of the ME 
and it can be represented by the one-dimensional non-linear dynamical system = 
^r(7j( m_1 )) with initial condition ?/ 0 ) = 1+a g/j 0 , and where mapping function *(•) is 
defined by 

= — / »\ / 7 ~ v\T (4-26) 



l+4aAE[ 7 /(,f 7 ) (l-/(i,g7))] 



It can be easily checked that, since /(•) has range [0, 1/2] and it is non-increasing, then 
^(rj) has range [rj(°\ 1] and it is non-decreasing. Therefore, the dynamical system 



3 In the large-system limit and under mild technical conditions this assumption is valid, as shown 
rigorously in [43]. 



4.3 Large system asymptotic analysis 
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Figure 4.4: Evolution of the ME given by asymptotic analysis for CC(5, 7), a = 2, 
and E b /I 0 = 5dB 



defined by ^(rj) with initial condition rj^ has at least one stable fixed point in the 
interval [rj( 0 \ l] , and the sequence {f/ m ^}^=o * s non-decreasing and upper bounded 
by 1. The limit rf = lim m _ KX3 r/ m ) < 1 exists and is equal to the left-most stable 
fixed point of the system in the interval [r]^\ 1]. 

As an example, Figure 4.4 shows the function for channel load a = 2, 

Eb/Io = 5dB, constant instantaneous received power (i.e. j u = 1 for all users), and 
the 4-states convolutional code of rate 1/2 and octal generators {5, 7} (denoted in the 
following by CC(5,7)). For the sake of comparison, we show also the evolution of 
the same system when the BCJR decoder provides soft extrinsic estimates as proposed 
in [30]. The *(ry) function in this case is derived in [37]. The limit jf in this case is 
very close to 1, meaning that all users attain near-single-user performance. Notice that 
soft feedback assures a faster decoder convergence to the single-user performance. 

When the channel load is increased, the * curves are modified so that for a certain 
threshold channel load a, the curve corresponding to weighted hard decisions is tan- 
gent to the diagonal, as shown in Figure 4.5. This means that the system has reached 
its maximum load and is not able to converge to near-single-user performance. On the 
contrary, the system using soft decisions still converges to single user performance. 
This shows that SISO decoding and soft feedback also provides a higher threshold 
channel load, i.e., a larger overall maximum achievable spectral efficiency of the sys- 
tem. 
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Figure 4.5: Evolution of the ME given by asymptotic analysis for CC(5, 7), a = 2.35, 
and E b /I 0 = 5dB 

4.4 Implementation of the proposed receiver 



4.4.1 The basic algorithm 

The performance of the proposed receiver depends on the computation of the weight- 
ing factors /?£ m ^ which do depend on a reliable SINR estimation. Driven by the asymp- 
totic analysis of the previous section, we propose to compute the weighting factor for 
the 11 th user at iteration m as 



= ! -2/ (SINR(™- x )) 



(4.27) 



where /(•) is the (known) SER code characteristics, and SINR^ is the estimated 
SINR at the it* decoder input of the iteration. In order to estimate the input SINR, 
we can use the estimator proposed in [44], given by 



SINR*? 1 ) = 



1 



N Sn=0 ^ M 



(4.28) 



- 1 



;(m) _ zi m) 



where = §j2- 

In [44] it is shown that 



£ m) = ^E|4 m) M| a -l*J 



n=0 
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is an unbiased estimator of the residual MAI plus ISI plus noise variance at the de- 
coder input if Cu^i 71 } i s uncorrected with the desired variable a u [n]. Remarkably, in 
this case, the MSE of this estimator is very close to that of the ML estimator assum- 

ing known useful symbols, given by J2n=o z v™\ n ] ~ \9u\ 2(l u[ n ] (which is not 
applicable here, since the symbols a u [n] are unknown). 

When the symbol estimates [n] are provided by a decision statistics contain- 
ing the current observation interval, such as in the Viterbi decoder, then the residual 
interference term given by (4.14) is conditionally biased given a u [n], that is, 



E 



ti m) [n]\a u [n]\ =6^a u [n 



(4.29) 



where the bias coefficient 5^ is non-positive and depends on the system parameters 
and on the user and iteration indexes as shown in [37, 38]. The bias is not negligible 
especially for high channel load a even in the absence of ISI and the bias reduces the 

/ Am) \ 2 

energy of the useful signal at iteration m by a factor f 1 + ^ p i . 

On the contrary, if the symbol estimates ai m ^[n] are provided by estimates not 
containing the current observation interval, i.e., they are based on the decoder extrinsic 
information [37, 45], then in the limit for large block length (i.e., TV oo) and random 
interleaving, the residual interference is conditionally unbiased, i.e. 



E 



= 0 



Ci m) [n]\a u [n] 

We can rewrite the input of the IX th decoder at iteration m given in (4.13) as 

4 ro) W = (iff J 2 + *£ m) ) ««W + Cl m) M + ««M (4-30) 



where Cu^M * s uncorrected with a u [n]. 

e presence 

(>j 2 +*£ m) ) 



Hence, the true SINR in the presence of bias is given by 

i\oJ z + sr') 2 

SINR u ( m ) = 



(>j 2 +^ m) y 



E |d m) M+Mn]| 2 ] 



(4.31) 



where we define = E 



CL m \n] +v u [n] 



. Now, the SINR estimator proposed in 



(4.28), in the presence of bias and for large N, converges in probability to 



SINR^ 



1 



E 



Ci m) \n]+vu\n] 



1 



Since 5^ < 0 (i.e., the bias tends to decrease the useful signal term), we conclude 
that the estimator (4.28) tends to overestimate the SINR at the decoder input. As a 
consequence, the weights computed according to (4.27) are mismatched in the 
presence of bias. 
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4,4,2 Estimation of the bias 

In order to overcome this problem, a better SINR estimation taking into account the 
bias of residual interference is required. The term cr| in (4.31) appears also in the 
variance of the decoder input signal, given by 

= (l»J 2 + 4 m) ) 2 + < 7 C (4-32) 

and the bias term 6^ is contained in the expression of the correlation between the 
symbol estimates and the decoder input, that, recalling (4.30), is given by 

<f>™ = E[zW[n]a u [n] 

= (l-2/(siNR u ( m - 1 ))+r(siNR u ( m " 1 ))) (|ffj 2 + 4 m) ) 

(4.33) 

where r(-) is a known characteristic function of the convolutional code, proportional to 
the correlation between the symbol estimates and the interference (see Appendix 4.A). 
Using (4.31) and (4.32) the correlation in equation (4.33) can be rewritten as 




(4.34) 



In practice, an approximation of a\ u ' and ^„ ; can be computed as follows 

"L"" 2 - ^E(4 m| N) 2 

and (4.34) can be solved numerically for 6^. Let 6^ be the estimated bias, solution 
of (4.34). Then the estimated SINR is given by 




(4.35) 



Figure 4.6 refers to a flat non-fading system, (i.e. j u = 1 for all users), with 
U = 32 users and spreading factor L — 16 (corresponding to the channel load a = 2), 
Eb/Io = 5dB using the convolutional code CC(5, 7). Users are chip-asynchronous and 
their relative delays are uniformly distributed over an interval spanning one symbol. 
The true and estimated SINR vs. the iterations for a given user are shown. The curve 



y (m) 



[n] 
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10 15 
Iterations (m) 



Figure 4.6: True and estimated SINR for "Viterbi (SINR est.)" and "Viterbi (SINR- 
BIAS est.)" decoders, using CC(5, 7), U = 32, L = 16, and E b /I 0 = 5dB 



labeled "Viterbi (SINR est.)" refers to the proposed iterative receiver that estimates the 
SINR using equation (4.28) while the curve labeled "Viterbi (SINR-BIAS est.)" refers 
to the system that estimates both the bias and the SINR, using (4.35). Notice that in 
the first few iterations the SINR estimator given by (4.28) overestimates the true SINR 
up to several dBs. Instead, the SINR estimation given by (4.35) is very close to the 
true SINR. 



4.4.3 The ping-pong effect and its compensation 

As it has been shown in the previous Sections, the proposed receiver allows system 
loads up to a certain threshold above which the system cannot approach the single user 
performance. In such high load situations, the system parameters as BER, SER, ME, 
and bias tend to oscillate between two convergence patterns [35]. This phenomenon 
is called ping-pong and it is related to the bias in the residual interference term. In 
fact, it does not appear when feedback is obtained from a SISO decoder extrinsic 
information [37]. 4 

A further investigation reported in [39] showed that such a bistable situation is due 
to a fixed subset of the estimated symbols that flip when passing from one iteration to 
the next, while the estimated symbols in the complementary subset do not change. A 
countermeasure to this problem proposed in [39] consists of introducing a perturbation 
into the bistable situation, by feeding back to the IC stage the average of the estimates 



4 For the extrinsic-based schemes, there obviously exists a threshold load above which single-user 
performance cannot be achieved, but no oscillatory behavior appears. 
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Figure 4.7: Bias in the statistic of z£ m) using CC(5, 7), U = 40, L = 16, and E b /I 0 = 
5dB. The figure shows the oscillations due to the ping-pong effect 

obtained from the two previous iterations. Thus, the contribution of the flipping sym- 
bols (considered as not reliable) is mitigated. 

By considering this idea, (4.20) is modified as 



^ (m) = E 



7 



(4.36) 



where the two previous estimates are now weighted with the coefficients /?j m ^ and 
. Now, we can optimize the weighting factors ft[ m ^ and as functions of the 
channel energy 7 in order to minimize the expected interference variance at every 
iteration. In analogy with what done before, we seek the solution of the optimization 
problem 



min E 

/?l=/?l(7),/fe=/?2(7) 



2 






7 



(4.37) 



After straightforward algebra (we skip the details for the sake of space limitations), we 
obtain the solution 



A*(7) 
&*(7) 



A-BC 

1-C 2 
B-AC 

1-C 2 



(4.38) 
(4.39) 
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where 

A = E[aa< m )] =l-2e( m )( 7 ) 

B = Ejaa^" 1 )] = 1 - 2e( m " 1 )( 7 ) 

C = E [iWal" 1 - 1 '] 

In practice, the weights f}^ and of user u at iteration m are computed as follows: 

• for m = 1, jSgJ is irrelevant (since fi(°> = 0) and p[]l = 1 - 2/(SINRi 1) ) 
where SINR^ is computed according to (4.35). 

• for m = 2, 3, . . we jet C = i ^ M^T"^ = 1-2 /(SINRL m) ), 
and B = 1 - 2/(SINRl m_1) ). Thus, /?{™ } and 0$ are given by (4.38) and 
(4.39), respectively. 

Figures 4.7 and 4.8 illustrate the bias and the true SINR plotted versus the iterations 
for a system using Viterbi decoding, with unfaded chip asynchronous flat channel, 
(i.e. j u = 1 Vu) U = 40, a = 2.5 and in the same conditions of Figure 4.6. The 
"Viterbi (SINR est.)" and the "Viterbi (SINR-BIAS est.)" decoders do not converge 
to near-single-user performance and show the ping-pong effect in the bias and in the 
SINR. The receiver that weights the two previous iterations, denoted in the Figures by 
"Viterbi (2-feedback)", converges faster to near-single-user performance and does not 
show oscillations. 

This behavior can also be seen in Figure 4.9 where the BER provided by the 
proposed receivers is shown, for an unfaded chip asynchronous flat channel, and the 
same conditions of Figure 4.8. The curve labeled "Viterbi (ft = 1)" denotes the re- 
ceiver that feeds-back hard decisions (i.e. assuming = 1 for u = 1, . . . , U and 
m = 0, 1, . . .) and the bold horizontal line represents the single user performance. The 
receiver "Viterbi (2-feedback)" converges to single user performance in 6 iterations. 
In these conditions, it converges even faster than employing BCJR decoders providing 
soft extrinsic output, denoted by "BCJR (soft EXT)". 

Finally, Figure 4.10 compares the performance of several receivers in terms of 
speed of convergence to near-single-user performance for E^/Io = 5dB, L = 16, and 
for 7 W = 1 Vu. As test of convergence, we consider the number of iterations required 
to reach an average ME larger than — O.ldB. The basic hard feedback Viterbi receiver 
has maximum channel load a = 1.5. The weighted hard feedback Viterbi with SINR 
and bias estimation improves this limit to a = 2.4. Eventually, the receiver "Viterbi 
(2-feedback)" outperforms any other considered receiver and has limit load a = 3.7 
(corresponding to a spectral efficiency of 1.85 bit/s/Hz). The dashed curves refer to 
the receivers employing BCJR decoders and extrinsic information feedback and are 
shown as a reference. 

We should remark that higher spectral efficiency could be achieved by employing 
more powerful codes (e.g. turbo codes) in the proposed iterative decoder. However, 
since our aim was to find a low complexity solution, we only considered very simple 
convolutional codes as the CC(5,7). So, despite of the good performance shown in 
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Figure 4.8: True SINR of several multiuser receivers using CC(5, 7), U = 40, L = 16, 
and E b /I 0 = 5dB 

terms of spectral efficiency, the proposed receiver works very far from capacity. For 
this reason it makes sense to compare the performance of our receiver to the "single- 
user performance", the performance that the system could achieve if no interference 
were present. 



4.5 Performance in frequency selective fading channels 

In this Section we consider frequency selective channels and asynchronous users. 
Thus, the ISI term appears in (4.12). In presence of ISI, the system converges to 
near single-user Matched Filter Bound (MFB) performance [10] for a channel load a 
generally lower than the limit for flat synchronous channels. Indeed, the presence of 
ISI slows down the convergence to single user performance and allows lower channel 
loads. Two different cases are considered: 

• constant instantaneous power for all users; 

• constant average power for all users. 



4.5.1 Constant instantaneous power 

This case is representative of perfect fast power control that fully compensates for the 
instantaneous (slot-by-slot) effect of fading, so that j u = 1. Thus, 

| ff J 2 = £, Vu = 1,...,U (4.40) 



4.5 Performance in frequency selective fading channels 
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Figure 4.9: BER of several multiuser receivers using CC(5, 7), U = 40, L = 16, and 
E b /I 0 = 5dB 




Figure 4.10: Number of iterations required to attain an ME> — O.ldB for several 
multiuser decoders using CC(5, 7), L = 16, and E^/Iq = 5dB 
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Table 4.1: MIP of the considered channels 



MIP 


Delays (T c ) 


Variances (dB) 


CHI 


(0, 1.2, 34.2, 49.5, 65.7, 76.8) 


(-2.5, 0, -12.8, -10.0, -25.2, -16.0 


CH2 


(0, 16) 


(0, 0) 


CH3 


(0, 16, 32) 


(0, 0, 0) 



Under these conditions, we evaluated by simulation the degradation due to the ISI with 
respect to the results shown in the previous Section. Figure 4. 1 1 shows the number of 
iterations needed to reach ME larger than — O.ldB, as a function of the channel load, 
for Eff/Io = 5dB, L = 16, the "Viterbi (2-feedback)" receiver, and for three different 
channels whose Multipath Intensity Profile (MIP) [16] is shown in Table 4.1. CHI is 
specified as a standard UMTS test channel in [21] while CH2 and CH3 were created 
"ad hoc" in order to test the system in severe ISI conditions. 

CHI shows a negligible degradation with respect to the single-path case while, 
for the channels CH2 and CH3, the maximum achievable channel load is decreased 
to a = 3.3 and a = 2.9, respectively. Notice also that, in the case of strong ISI and 
low channel loads, more iterations are required for the convergence with respect to the 
single-path case. 

4.5.2 Constant average power 

This case is representative of a slow power control system that cannot compensate for 
the instantaneous fading, but maintains a fixed average received power for each user. 
The average received power is assumed to be the same for all users, that is 

E[\ 9u \ 2 ] = £ ; Vu=l,...tf (4.41) 

Figure 4.11 shows the number of iterations needed to reach ME larger than —O.ldB 
as a function of the channel load and for the channels given in Table 4.1. In this case, 
two contrasting phenomena characterize the convergence to single-user MFB perfor- 
mance. On one hand, faded users have little impact onto unfaded users, since they are 
received at much lower power. On the other hand, even if unfaded users can be reliably 
estimated and subtracted from the received signal, faded users have an instantaneous 
SINR that is too low for achieving a small SER, therefore, their contribution cannot 
be perfectly eliminated even in the absence of MAI. The fact that the performance 
for slow power-control are worse than their fast power-control counterparts (see Fig- 
ure 4.1 1) indicates that the effect of uncanceled faded users dominates the performance 
of the receiver. Moreover, Figure 4.11 shows also another interesting fact. Namely, 
the degradation of the receiver performance due to uncompensated fading is larger for 
channels with little multipath diversity (e.g., single-path, CH2 and CH3), while it is 
smaller for channels with rich multipath diversity (e.g., the UMTS test channel CHI). 
This is intuitively clear, since when the multipath diversity is rich, then the random 
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fluctuations of the instantaneous received power \g u \ 2 are reduced, and the fraction of 
users that can be reliably estimated and subtracted is larger. 

In general, the average bit error rate (BER) for a finite-dimensional system can 
be computed by Monte Carlo simulation, by averaging over a large number of frames 
and channel realizations. In order to reduce the complexity of BER computation, we 
propose a semi-analytic approach. Assuming symmetric users, with the same channel 
statistics and average received power, the average BER (averaged with respect to both 
the channel statistics and over the user population) can be bounded by 



where /&(•) is the bit error probability vs. decoder input SNR function of the employed 
convolutional code, where (a) follows from Jensen's inequality applied to the convex 
function /&(•) and where we define the average ME ?j( m ) as the arithmetic mean of the 
ME of all users. Notice that, under mild conditions, for randomly spread CDMA in 
the large-system limit the ME converges to a constant independent of the user index, 
and inequality in (4.42) holds with equality. However, for a finite-dimensional system 
and given spreading sequences the above provides only a lower bound. 

The expectation in the last line of (4.42) is with respect to the normalized channel 
energy 7 (assumed identically distributed for all users). The evaluation of the above 
lower bound is much less computationally intensive than full Monte Carlo simulation 
of the whole system. In fact, the average ME ?/ m ) can be obtained by running short 
simulations, since it requires a much smaller statistical sample to converge than the 
average BER. Then, the expectation with respect to 7 can be obtained by either another 
Monte Carlo simulation or, if the distribution of 7 is known, by numerical integration. 

Figure 4.12 shows the comparison between the full Monte Carlo simulation and the 
semi-analytic approach for a system with U = 32, L = 16, and channel CH3. Solid 
lines show the BER obtained by Monte Carlo simulation for m = 1, 2, 3, 6 iterations 
respectively, plotted versus E^/Iq. The dashed lines show the results obtained using 
the semi-analytic approach. The thick solid line corresponds to the single-user MFB. 
For low SNR the system does not converge to the single user MFB even for a large 
number of iterations while for high SNR the convergence is obtained in a few itera- 
tions. We notice that the semi-analytic approach yields fairly accurate results already 
for such a small system. 




u=l 




(4.42) 



4.6 Conclusions 



We proposed a low complexity iterative turbo equalizer and multiuser decoder/detector 
for TD-CDMA systems, characterized by convolutional coding, hard-output Viterbi 
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Figure 4.1 1: Performance of the "Viterbi (2-feedbacks)" multiuser decoder in the pres- 
ence of multipath fading channels using CC(5, 7), L = 16, and E^/Iq = 5dB. FPC 
and SPC denote fast and slow power control conditions, respectively 



decoding and weighted feedback. From a rigorous large-system analysis of synchronous 
users and flat channels we gained the rationale for the optimization of the feedback 
weights. In the proposed iterative scheme, however, the presence of Viterbi decoders 
(not providing extrinsic information) produces biased statistics after the IC stage. 
Hence, we proposed a method for estimating the bias and improving the performance 
of the basic iterative decoder. Moreover, in order to cope with the ping-pong effect, 
we proposed a modified algorithm where a weighted sum of the decisions made in the 
two previous iterations is used for interference cancellation. The modified algorithm 
outperforms the the basic PIC receiver and even the receiver based on SISO decoding 
and soft feedback, with lower complexity. 

The proposed receiver was validated in a realistic scenario, including asynchronous 
users, frequency selective fading channels and power control. Simulation results show 
that the receiver is robust to severe ISI conditions, even though its performance is de- 
graded by uncompensated fading. However, this degradation is not very evident in the 
presence of sufficient multipath diversity. Eventually, a simple and fast semi-analytic 
approach to compute the BER was proposed and its behavior has been compared with 
the results obtained via Monte Carlo simulation. 

As a concluding remark, we would like to point out that the proposed receiver 
structure is fully suited to be implemented in a UMTS-TDD base station, as an al- 
ternative to conventional single user decoding and high-complexity convolutional or 
turbo codes proposed in the current standard. 



4.A Definition of /(•), f b (-) and T( ) 
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Figure 4.12: BER provided by the "Viterbi (2-feedbacks)" multiuser decoder for chan- 
nel CH3 using CC(5, 7), U = 32, and L = 16 

4.A Definition of /(•), /&(■) and T( ) 

Consider a transmitter that maps the information bit sequence {&[?]} onto a sequence 
of coded binary BPSK symbols {a[n]} and sends over the AWGN channel defined by 

y[n] = a[n] + v[n] 

where v[n] is Gaussian noise with distribution N(0, cr 2 ). The signal-to-noise ratio is 
SNR = ^ • The receiver decodes the signal y[n] providing hard decisions a[n] on 
the transmitted symbols a[n] and hard decisions b\j] on the information bits b[j]. We 
define 

1. The SER vs. SNR characteristic 

/(SNR) = ^EftM'M] (4.43) 



2. The BER vs. SNR characteristic 



1 - E [b\j]b\j]] 



/ 6 (SNR) = l - i- (4.44) 

3. The correlation characteristic between decisions and noise 

T(SNR) = E[o[n]i/[n]] (4.45) 
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The above functions depend uniquely on the code considered (i.e., they are "charac- 
teristics" of the code) and, although generally unknown in closed form, they can be 
pre-computed by Monte Carlo simulation and stored as look-up tables in the receiver 
memory. Therefore, the impact of real-time evaluation of these functions on the overall 
receiver complexity is negligible. 



Chapter 5 



Front-end signal processing 



We present low-complexity algorithms for transmitter and receiver front-end suited 
to the implementation of software defined radio (SDR) terminals. The proposed al- 
gorithms make the processing sampling frequency independent of the symbol rate of 
the digitally modulated signal and use the "IF -sampling " technique for D/A and AID 
conversion. As a case-study, we consider a training-based joint multiuser channel es- 
timation and we show that our front-end algorithms work nicely when coupled with an 
efficient FFT-based joint channel estimator. 



5.1 Introduction and motivations 

Because of the proliferation of mobile communication systems, the implementation 
of Software Defined Radio terminals able to reconfigure themselves to handle several 
different standards is a key issue in order to provide universal connection to the users. 
At the physical layer, the careful design of algorithms for the transceiver front-end 
signal processing seems to provide the required flexibility. The basic idea is to make 
the standard-dependent base-band signal processing independent on the Intermediate 
Frequency which is characterized by the hardware. Moreover these algorithms must 
be suited to the implementation on a programmable CPU. The goal is to perform op- 
erations like channel selection [5, 6], up- and down-conversion [7], synchronization 
and detection in the all-digital domain, by using high performance Digital Signal Pro- 
cessors (DSPs) while, traditionally, these operations are performed by dedicated hard- 
ware. 

In this Chapter, we propose some front-end algorithms for the class of linearly- 
modulated digital signals, including most today's and future mobile communications 
standards, like GSM, IS-54, IS-136, IS-95, and DECT (see e.g. [1] and references 
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therein), UMTS (both FDD and TDD modes) [13], CDMA2000 (see e.g. [22], [46]) 
and EDGE (see e.g. [47] and references therein) l . 

We assume that a multi-band RF section takes care of translating the desired signal 
from a fixed Intermediate Frequency carrier to the required Radio Frequency carrier, 
and vice-versa. Then, on the transmitter side we are concerned with the efficient gener- 
ation of an IF analog signal from the base-band digital signal while on the receiver side 
we are concerned with the efficient generation of a digital signal from the IF analog 
signal. 

Different options are compared in terms of their complexity and performance. 
Complexity is measured in real operations per sample. This can be translated into 
actual processing speed limits once the DSP number of instructions per second is 
known. Performance is given in terms of the Signal-to-interference plus Noise Ra- 
tio (SINR) at the output of the receiving filter matched to the modulation elementary 
pulse, assuming that the transmission channel is pure additive white (not necessar- 
ily Gaussian) noise. In this case the interference should be intended as the result of 
the whole signal distortion generated through the transmitter and receiver chain. The 
SINR is expressed as a function of £/7 0 , where £ is the average symbol energy and 
7o denotes the (frequency-flat) noise power spectral density, in the complex base-band 
equivalent model. By modeling interference plus noise at the output of the receiver 
matched filter as Gaussian, the output SINR can be easily translated into actual error 
probabilities for different modulation schemes with and without channel coding. In 
this way, both complexity and performance measures are independent of the consid- 
ered particular modulation scheme and the comparison between different solutions is 
consistent. 

In order to validate the proposed algorithms on an actual system, we consider 
the training-based multiuser joint channel estimation scheme of UMTS-TDD (see 
e.g. [18], [19], and [20]). We recall briefly an efficient Fast Fourier Transform (FFT)- 
based implementation for this scheme and we apply it to the signal generated by our 
transmit/receive front-end. Performance for a particular multipath test channel are 
given in terms of the actual measured Matched Filter Bound (MFB) at the output of 
the estimated channel matched filter versus the MFB of an ideal system. 



5.2 Transmitter front-end 
5.2.1 Linearly modulated signals 

The continuous-time complex envelope of a linearly modulated signal is given by 



where d[k] is a sequence of modulation symbols belonging to some complex alpha- 
bet (e.g., PSK, QAM [16]), T s is the symbol interval, and ip(t) is the symbol-shaping 

^SM and DECT make use of GMSK modulation, which is non-linear [1]. However, in several 
current implementations a linear approximation of GMSK is used. 




(5.1) 



k 
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pulse, bandlimited over [— Despite its simplicity, several cases of in- 
terest can be expressed by (5.1): 

• In a Direct-Sequence CDMA system [16], T s should be interpreted as the chip 
interval and d[k] = a[|_A;/LJ]s[A;], where a[m] is the modulation symbol, 
s[k] is the chip and L is the spreading gain (this generalizes trivially to sys- 
tems with several spreading layers, like IS-95). 

• In a multi-antenna multi-user system (e.g., downlink with beamforming) with Q 
antennas and U users the signal sent to the antenna is 

x^(t) = ^2d^[k]i;(t-kT s ) 
k 

where 

u 

d^[k]=Y,Wi q) du[k] 

Wu^ is the q ih antenna weight coefficient for the u ih user and d u [k] is the 11 th 
user symbol (or chip). 

In a digital transmitter, the signal x(t) is the output of a D/A converter which takes 
as input the discrete-time signal 

x[n] = x(n/f s ) 

with sampling frequency f s > W^. In classical I-Q modulators, the continuous base- 
band components xi(t) = Ke{x(t)} mdxQ(t) = Im{x(t)} are generated by low-pass 
filtering the output of two separate D/A converters, and the IF signal 

xi F {t) = Re {x{t)et 2nflFt } (5.2) 

is produced by mixing xi(t) and XQ(t) with in-phase and quadrature IF carrier signals, 
and by summing the modulated real signals [16]. This approach requires two D/A 
converters, two low-pass filters, two analog mixers and one adder. 

Another approach consists of producing a sampled version of the IF modulated sig- 
nal xip(t) by using a sampling rate f s greater than 2/if- The continuous-time signal 
is obtained by band-pass filtering the output of a single D/A converter. This solution 
(for the receiver front-end) is described in [7] and [48]. However, since intermediate 
frequencies usually range between tens of MHz up to 100 MHz, this approach is ex- 
tremely computationally intensive as it requires the generation of signal samples and 
the multiplication by the carrier signal at extremely large sampling frequency. 

In the following, we propose a transmit front-end algorithm allowing a sampling 
frequency of the order of the base-band signal bandwidth (and not of the order of the 
IF carrier), no explicit multiplication by the carrier signal and a single D/A converter 
and analog filter centered at /ip. 
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5,2,2 IF-sampling and up-conversion 

We need to generate the IF modulated discrete-time real signal 

x'[n] = Re \x[n\^ 2 < h *l^ n } (5.3) 

from the base-band signal x[n]. We choose throughout this Chapter the sampling rate 
f s according to the expression 

/IF 

fs — » x + , A f° r a positive integer I (5.4) 
t ± 1/4 

With this choice (5.3) reduces to modulating the base-band signal by / 5 /4 and taking 
the real part that in turn corresponds to 

x f [n] = Re{i ±n x[n]} (5.5) 

In this way, the periodic spectrum of x f [n] has a spectrum replica centered at /ip. After 
ideal D/A conversion, a pass-band filter centered at /if removes the other replicas, 
generating the desired IF modulated analog signal. The discrete-time modulation by 
/ s /4 in (5.5) has a negligible computational cost since it corresponds to alternately 
changing the signs of the real and imaginary part of x[n\. Finally, in order to avoid 
aliasing the sampling rate must satisfy the condition /, > 2W^. 



5.2.3 D/A conversion 

In the above description we assumed an ideal D/A converter with flat frequency re- 
sponse. Unfortunately, actual D/A converters show a low pass behavior that can be 
approximated with a frequency response of the form sinc(/// 5 ). The spectrum of the 
continuous-time signal at the D/A converter output corresponding to the discrete-time 
input signal x f [n] takes the form 

X w (f) oc ^X' (^ if/f '- i) ) sinc(/// a ) (5.6) 

% 

where 

X'(^flf*) = ^2x f [n]^ 2nn f^ s 

n 

is the Discrete-time Fourier Transform (DFT) of x f [n]. Hence, the spectrum replica 
located at frequency /if is attenuated and distorted by the D/A impulse response as 
shown in figure 5.1(a). 

A way to reduce the attenuation consists of using a D/A converter working at rate 
fd = ^d/a/s* where L^/A i s a suitable integer, and up-sampling x f [n] by the factor 
L D /A- By choosing L B / A such that 3> /if, the spectrum replica around /if falls 
inside the first lobe of the D/A frequency response since its first zero is located at /^ 
(see figure 5.1(b)). Unfortunately this approach reduces the average signal energy per 
sample by a factor £d/a- Thus, in order to compensate the signal attenuation result- 
ing from high values of Ld/A> high amplification gains may become necessary in the 
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(a) D/A frequency response without up-sampling 
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(b) D/A frequency response with up-sampling 
Figure 5.1: D/A frequency response and IF signal spectrum 

IF/RF conversion stage which can lead to significant non-linear distortion. To over- 
come this problem one may pre-compensate the linear response of the D/A converter 
by introducing a pass-band FIR filter between the up-sampler and the D/A converter. 
The filter must be designed in order to enhance the spectrum replica at IF while at- 
tenuating the other replicas. This would avoid high amplification gains in the IF/RF 
stage yielding less signal distortion. In the following we address a low complexity 
implementation of such a filter. 

Low Complexity Implementation. Denote by x n [n] the up-sampled version of x'[n] , 
given by 

x „ [n]= \Ak] forn = L D/A fc ( . ?) 

^0 otherwise 



and by /&d/aM the pass-band FIR filter impulse response. The filter output is given 
by 

x"'[n] = h D/A [m]x"[n - m] (5.8) 

m 

If the filter impulse response has length Ld/a> there is only a single non-zero term in 
the sum in the RHS of (5.8). Then, x m [n] = /id/aM^W where k = L n /^D/Aj 
and m = n modulo Ld/A- Therefore, the computational cost of the filtering operation 
consists of one real product per output sample at rate Since fd is large, this might 
still be too complex for practical applications. Then, the filter coefficients can be 
constrained to be in the set {0, ±1}. In this way products are avoided (change of 
sign comes basically at no cost). The problem of optimizing the FIR filter reduces 
to exhaustive search for the best frequency response over all possible filter vectors of 
length Ld/a with ternary elements {0, ±1} (a total of 3 L °/ A possibilities). 
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5,2,4 Signal resampling 

In the scheme described above, the sampling frequency f s depends on /if and on the 
signal bandwidth W^p. In a SDR system, we can think of as the maximum signal 
bandwidth over all handled signal formats, while /if is fixed and depends on the analog 
hardware. In any case, f s should be independent of the symbol rate R s = 1/T S of the 
particular signal format. When f s ^ R s , the transmitter must produce the signal x[n] 
at rate f s from the discrete-time symbol sequence d[n] at rate R s . We refer to this 
general problem as "resampling". 

In [49, 50, 51, 52] resampling by rational factors is discussed. A general sampling 
rate conversion scheme by a rational factor L u jL^ consists of an up-sampler by a 
factor L w , a low-pass filter and a down-sampler by a factor L^. The low-pass filter can 
be efficiently implemented by a polyphase filter bank [49] composed by L u phases, 
where each filter phase is sampled at rate R s . Unfortunately, in general the ratio f s /R s 
is not a rational number, therefore f s can be only approximated by j^R s by a careful 
choice of the intes ers L u and L^. 

Since the shaping filter ip(t) of the modulation scheme is bandlimited, it can be 
used as low-pass filter in the resampler. In this way, shaping and resampling are inte- 
grated in a single step with considerable saving in computational complexity. Assume 
that ip(t) is approximated by a finite-support impulse response spanning N^p sym- 
bol intervals (e.g., by windowing [51]) and let ^[i] = ift((iL u + £)/(L u R s )) with 
I = 0, . . . , L u — 1 and i = 0, . . . , — 1, denote the i ih sample of the I th phase of 
the polyphase filter bank implementation of the shaping discrete-time filter. Then, the 
output at rate L U R S of the polyphase pulse-shaping filter bank is given by 



where m = [k/L u \ and £ = k modulo L u . 

A low cost approach for achieving the desired non-rational sampling rate conver- 
sion consists of using a suboptimal interpolation method (e.g., nearest neighbor, linear, 
cubic or cubic spline interpolation) at the output of the polyphase filter bank. In gen- 
eral, the interpolator computes the time epoch t n of the desired 71 th output sample and 
compares it to the input-sample time-instants. The cost of this operation does not de- 
pend on the interpolator type and can be assumed to be (roughly) 1 addition per output 
sample. On the contrary, the computation of the output sample value depends on the 
interpolation method. Therefore, when comparing the complexity of the interpolation 
schemes, we shall always refer to this cost, neglecting the part common to all meth- 
ods. Notice that the polyphase filter-bank coefficients can be pre-computed thus the 
up-sampling factor L u only affects the memory occupancy and not the complexity of 
the interpolation algorithm. In general the higher L u the higher the interpolation accu- 
racy will be. In particular the more the complexity of the interpolation algorithm the 
lower values for L u are required to achieve the same interpolation accuracy. 

Nearest neighbor (NN) interpolation is the simplest technique, since it approxi- 
mates the desired output sample to the time-nearest input sample. In fact using nearest 
neighbor interpolation the overall cost of the resampler reduces to storing in memory 
the shaping filter bank. The nearest neighbor interpolator produces the output sequence 




(5.9) 
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x[n] at sampling rate f s as 

x[n] = v[k] 

where 

and where t n = to + n/f s is the desired sampling epoch (to is a fixed time-offset). 
The computational complexity is then the cost of the inner product given by (5.9) per 
output sample. Normally the shaping filter coefficients are real while, in general, the 
symbols d[k] are complex. Then, the resulting computational complexity is 2N$ mul- 
tiplications and 2(N^p — 1) additions per output complex sample. However, as indicated 
in (5.3), only the real signal x f [n] is needed in our transmit front-end. This is equiv- 
alent to computing alternately xi[n] and a;Q[n]. Therefore, the actual computational 
complexity per output sample reduces to multiplications and — 1 additions. 

It is clear from (5.10) that the accuracy of nearest neighbor interpolation criti- 
cally depends on L u . However, since memory is a very cost-effective resource, in the 
transmitter it is advisable to choose L u large enough and insist on nearest neighbor 
interpolation rather than consider a lower L u and implement a more complex interpo- 
lation method (e.g., linear, cubic, etc.). We shall see in next section that this conclusion 
does not apply to the receiver front-end. 



5.3 Receiver front-end 

5.3.1 IF-sampling and down-conversion 

The IF received analog signal yi^(t) is sampled by an A/D converter at rate f s . If 
fs > ZWjjj is chosen according to (5.4), because of the periodicity of the discrete-time 
signal spectrum, the resulting real sampled signal r[n] = yw(n/ f s ) is pass-band with 
a spectrum replica centered at / 5 /4 (although /if and f s at the receiver can be different 
from / IF and f s at the transmitter, for simplicity we use the same notation). 



5.3.2 Resampling 

To simplify synchronization and data detection, the received signal should be sampled 
at rate R f s , integer multiple of the symbol rate (or chip rate, in the case of CDMA). 
As in the transmitter, the system should be able to choose R f s independently of f s 
(within a certain range of supported signal formats). In addition, R f s /f s is non-rational 
in general and resampling is needed. A general approach to have a perfect transparent 
resampling system consists of resampling the I and Q components of the real pass- 
band signal r[n]. Once the I and Q components sampled at rate R f s are generated, we 
can remodulate the signal at frequency i^/4 by simple sign change (as done in (5.5)) 
in order to have a real pass-band signal sampled at the desired rate R f s . We consider 
the following two resampling approaches. 
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Figure 5.2: Receive resampling approach 1 



Approach 1 



This is a fairly similar resampling scheme to the one proposed for the transmitter (see 
section 5.2.4). Notation refer to the scheme of figure 5.2. The low-pass real inter- 
polation filter h is implemented as a polyphase FIR filter bank with L u phases (ht[i] 
denotes the I th phase 2 th coefficient, for t = 0, . . . , L u — 1 and i = 0, . . . , — 1). 
Let r'"[n] = r'£[n] + jr"'[n] denote the output of the polyphase filter bank, given 

by 

N h -1 

r"'[n]= £ htWHr-'rim-i] (5.11) 

where m = [n/L u \ and £ = n modulo L w . Using nearest neighbor interpolation 
produces the output sequence z[k] at sampling rate R f s as z[k] = r"'[n] where 



(5.12) 



and where ^ = to + k/R f s is the desired sampling epoch (to is a fixed time-offset). 
Finally, the resampled pass-band signal modulated at frequency i?' 5 /4 is given by 
y[k] = Re{j fc r"'[n]}. Since the interpolation filter is used to remove the high fre- 
quency noise components it can be chosen to be real. In this case y[k] can be directly 
written as 



y[k] = { 



N h -l 


~V[m 


i=0 
N h -l 




- £ Mi]Im{(-j) m - 


"V[m 


i=0 
N h -l 




!=0 


"V[m 




i=0 


"V[m 



for integer / 



(5.13) 

where m and £ depend on n, which is a function of k through the resampling epoch 
relation (5.12). 

After some algebra it is possible to show that the computation of y[k] reduces to 



5.3 Receiver front-end 



71 



r[»] 



r'[n] 



rate / s 



Re 




@ 


1 — 


Interp 


— H 





Im 




@ 




Interp 



rate 



Hi 



Re 



Figure 5.3: Receive resampling approach 2 



computing alternately the values 

N h /2-l 

Vo[k] = Yl ko\i]r[m-2i] (5.14) 

N h /2-l 

Vi[k] = E hi\i}r[m-2i-l} (5.15) 

where the modified (real) filter coefficients 

ht,o\i] = i-Whm 
ht+li] = (-lj^pi + l] 

can be precomputed. 2 

The overall complexity of this resampling approach amounts to an inner product 
of length (i.e., N^/2 multiplications and N^/2 — 1 additions) per output sample 
at rate R f s . 

One may notice that we could use as interpolation filter the same shaping filter 
tf) used in the transmitter. In this way, resampling and pulse shaping matched filter- 
ing could be implemented jointly. However, while this approach gives a complexity 
saving in the transmitter, it is not suited to the receiver because better low-pass FIR 
responses h, with the same number of taps, can be synthesized by using appropriate 
filter design methods (see e.g. [49]). Moreover, after filtering by ip(t) (which is nor- 
mally a real Nyquist pulse) and sampling at R f s > R s the resulting interference plus 
noise sequence is correlated even if the interference plus noise process at the receiver 
input is white. Correlation in the noise sequence might prevent some blind channel 
estimation techniques based on second-order statistics to work [53, 54] (or they might 
require the correlated noise to be whitened, with added complexity), and should be 
avoided. Then, we conclude that the interpolation filter h must be as close as possible 
to an ideal low-pass filter. 

Moreover, when the received signal is affected by linear distortion caused by a mul- 
tipath propagation channel h(t), a channel estimation scheme applied to the discrete- 
time resampled signal y[k] directly provides an estimate of the whole convolution 
g(t) = ip(t) * h(t). In this way the pulse-shaping matched filter is included in the 
channel matched filter. 



2 The actual value of y[k] is given either by ±yo[&] or by ±yi [k], where the sign depends on the parity 
of k and m. 
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Approach 2 

An alternative approach consists of exploiting the fact that in our system the base- 
band signal r f [n] is already oversampled with respect to the Nyquist rate. Thus the 
interpolation filter h can be avoided with a significant reduction of complexity for the 
same type of interpolation. Figure 5.3 shows the block diagram of this re-sampling 
system. 

The received signal is first converted to base-band, then the real and imaginary 
parts of r'[n] are treated separately and eventually recombined. Since r[n] is real, 
the real and imaginary parts of r f [n] = (— j) n r[n] have zeros in every other sample. 
This can be seen as the effect of high-frequency components in the spectrum of r'[n] 
originated by the spectrum replicas around ±/ 5 /2. Hence, interpolation is obtained by 
eliminating the zero samples (by downsampling by a factor 2) and by applying some 
interpolation algorithm to the non-zero samples. After recombination, the complex 
signal z[k] is remodulated and the real part is taken, thus obtaining y[k]. Notice that 
the whole processing of remodulating z[k] at frequency and taking the real part 
reduces to computing alternately the real and the imaginary part of z[k] and alternately 
changing the sign to obtain y[k]. Therefore, the overall complexity amounts to one 
interpolated real value per output sample at rate R f s . In this work we consider nearest- 
neighbor, linear, cubic and cubic spline interpolators (the details are given in [55] and 
references therein) whose computational cost is given by [55]: 

• 0 for nearest neighbor interpolation. 

• 2 additions and 1 multiplication per output sample for linear interpolation. 

• 1 1 additions and 8 multiplications per output sample for cubic interpolation. 

• 5 additions and 9 multiplications per output sample plus 4 additions and 3 mul- 
tiplications per input sample for cubic spline interpolation. When the ratio be- 
tween the interpolator input and output rates is approximately 1 the overall com- 
plexity can be considered roughly 9 additions and 12 multiplications per output 
sample. 

We shall show in the next section that approach 2 yields significantly better perfor- 
mance with respect to approach 1 even with lower complexity. 



5.4 Performance over additive white noise channel 

In this section we illustrate the impact of the various system parameters on the perfor- 
mance of the proposed front-end algorithms over an additive white noise channel. As 
performance measure we use the SINR measured at the output of receiver front-end 
versus the nominal £/Jq ratio, where the interference should be intended as the conse- 
quence of the signal distortion. In practice the SINR is defined in this case as the ratio 
between the mean square and the variance of the received signal at the output of the 
pulse-shaping matched filter. The SINR loss with respect to the nominal £/Jo is due to 
two effects: Inter-Symbol Interference and energy decrease of the useful samples. ISI 
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Figure 5.4: Receive resampling approach 1, = 20, Lrx = 16, and nearest neighbor 
interpolation 

is due to the shaping filter truncation and to the distortion introduced by resampling, 
while the useful samples energy decrease is only due to resampling distortion. 
We consider 

• the QPSK modulation 

• the root-raised cosine shaping filter proposed by the UMTS norm with roll-off 
factor r 0 ff = 0.22, approximated by a FIR spanning 10 symbol intervals 

• an IF carrier of 70 MHz 

• the symbol rate (chip rate) proposed by the UMTS-TDD norm, that is R s = 
3.84Mchips/s 

The transmitter and receiver sampling frequency is f s = /if/ (5 + 1/4) = 13.33 MHz, 
yielding f s /R s ~ 3.4722 according to (5.4). In approach 1 the low-pass filter h 
is a linear phase FIR filter designed by using the Remez algorithm [51] in order to 
minimize the distortion in the signal bandwidth. For the sake of notation we denote 
as Ltx and Lrx the upsampling factors for the transmitter and receiver resampling 
system respectively, instead of using L u as done in the previous sections. 

Figures 5.8-5.11 show the performance of the receiver front-end with resampling 
approach 2 using nearest neighbor, linear, cubic, and cubic spline interpolation at the 
receiver, respectively. Each figure shows the performance for transmitter up-sampling 
factors Ltx = 4, 8, 16, 32, 64. With large Ltx and cubic spline interpolation at the 
receiver the system yields negligible SINR degradation over the range [—20, 30] dB of 
nominal £//q- 
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Figure 5.6: Receive resampling approach 1, = 10, Lrx = 1, and linear interpola- 
tion 
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Nearest N. interp.; L Rx =2; N h =6 




Figure 5.7: Receive resampling approach 1, = 6, Lrx = 2, and nearest neighbor 
interpolation 




Figure 5.8: Receive resampling approach 2 and nearest neighbor interpolation 
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Figure 5.9: Receive resampling approach 2 and linear interpolation 




Figure 5.10: Receive resampling approach 2 and cubic interpolation 
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Figure 5.11: Receive resampling approach 2 and cubic spline interpolation 



The best performance with approach 2 are achieved by cubic spline interpolation 
whose complexity is approximately equivalent to that of approach 1 when a filter with 
Nfr = 20 coefficients per phase, Lrx = 16, and nearest neighbor interpolation are 
employed at the receiver (see figure 5.4). Although the complexity of the two ap- 
proaches is almost the same in this case, we remark that with approach 1 up to 320 
coefficients need to be stored in and loaded from memory. Moreover when the num- 
ber of filter coefficients becomes high because of large values of Lrx, the DSP has to 
handle memory accesses and loading data while doing computations and this fact can 
significantly reduce the actual computation speed. 

Approach 1 can achieve better performance than approach 2 at the cost of an in- 
creased complexity, e.g. by using the same Lr x = 16 but employing linear inter- 
polation instead. As shown in Figure 5.5 the performance gain is about 3 dB in the 
high £/7q region while the complexity is doubled with respect to the nearest neighbor 
interpolation and to approach 2 with cubic spline interpolation. 

We may also try to keep the same complexity as approach 2 with cubic spline inter- 
polation by using a filter of length = 10 samples per phase and linear interpolation. 
However this solution yields worse performance as shown in figure 5.6. 

In general, a way to reduce the complexity in approach 1 consists of reducing 
the length of the filter h. Anyway a shorter filter rejects less noise and introduces 
larger signal distortion. So when increasing Lrx this corresponds to increasing the 
temporal resolution of the interpolation and can result in poor performance. Indeed, 
high temporal resolution makes the interpolation follow the signal variations due to 
the high frequency components associated with noise and distortion. For instance, in 
the case of = 10 and linear interpolation we have seen that for Lrx > 1 the 
performance monotonically decreases as Lrx increases. For the same value of 
nearest neighbor interpolation yields the maximum performance for Lrx = 4. This 
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leads to conclude that for a given filter length there exists an optimal trade-off between 
the type of interpolation and the number of filter phases Lrx. 

Another example of this behavior is given in figure 5.7 where the performance is 
shown for Lrx = 2, = 6, and nearest neighbor interpolation. Larger values for 
Lrx worsen the performance for the reasons explained above. The complexity of this 
system is about the double of the one of approach 2 with linear interpolation (see figure 
5.9), but yields worse performance. 

We can conclude that approach 2 is actually an attractive solution for many practi- 
cal scenarios in terms of performance and complexity. 



5.5 Example: training-sequence based channel estimation 

The receiver front-end output signal, sampled at rate multiple integer of the symbol 
rate is suitable to be further processed by the receiver. In this section we test the effec- 
tiveness of the proposed front-end algorithms on a training-sequence based multiuser 
channel estimation procedure for block-synchronous CDMA (e.g., UMTS in TDD 
mode). 

In this scheme users are roughly synchronized to a common time-reference and 
transmit their training sequences at the same time (user timing errors are included as 
effect of the channel and accounted for in the estimation procedure). The maximum 
channel length (including possible timing errors) is C 0 ff chips (see Section 2.1.4) and 
the training sequence sent by each user is built from the same common base sequence 

- -iT 

s% = 5t[0], 5t[l], • • • , sjftr — 1] of length P tT chips, by adding a cyclic exten- 
sion of C 0 q chips. This solution allows a joint estimation of all users' channels if 
-Ptr > C 0 ffU, where U is the number of interfering users (see [18, 19], and with some 
modifications [20], and references therein for more details). 

Under these assumptions we can write the received signal sampled at frequency 
R f s = N C R S during the P tT symbols spanned by the training sequence as 

V = S ig + u (5.16) 

where 

y=[y[0], y[l],..., y[Pt r N c - 1]] T 

is the received signal, 

[T T T ] ^ 

is a vector containing the channel impulse responses of all U users, 



9u = [fti[0], ffu[l], ■ ■ ■ , 9u[C 0 &N c - 1]] 



is the 11 th user channel impulse response, and v is a vector of interference plus noise 
samples, assumed white. The Pt r N c x Pt r N c matrix St is defined as 



St = S t ®I Nc 



(5.18) 
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Figure 5.12: Estimation of the complex pass-band channel response 

where S t * s a circulant matrix containing all the possible cyclic shifts (by columns) of 
the base sequence s t . The matrix S t is also circulant and it is unitary similar [56] to 
the diagonal matrix diag(7), where 




N c times 



and where 7 is the discrete Fourier transform of s t - After some algebra, it is possible 
to show that the least squares estimation of the overall channel impulse response g is 
given by [20] 

f DFT{w} 1 

g = IDFTj (519) 

where DFT and IDFT denote direct and inverse discrete Fourier transforms and the 
ratio of two vectors should be interpreted as the element-by-element division. 

In the case of pass-band signals considered here, figure 5.12 shows qualitatively 
the spectrum Y(e^ 2n ^ Rs ) of y, and the spectrum corresponding to 7 for N c = 4. 
Since the transmitted signal is (approximately) bandlimited and centered around i?' 5 /4, 
the pass-band one-sided (i.e., complex) channel impulse response of all users can be 
estimated by setting to zero all coefficients of DFT{y} but those in the frequency 
range [i?' 5 /4 - /2, i?' 5 /4 + /2] . In this way, the computation of (5.19) is further 
simplified. 

The method applies to pass-band signals as well as to base-band signals, without 
complexity increase. Furthermore, pass-band signals do not need explicit demodula- 
tion. Demodulation can be automatically achieved by down-sampling the output of the 
pass-band symbol matched filter at symbol rate. 



5.5.1 Performance with channel estimation 

We used the training sequence defined in the 3GPP proposal for the UMTS-TDD 
mode [13] and the "Midamble type 2" defined in the Table 2.1. The sequence length 
is P tr = 192 chips and the assumed channel length is C 0 ff = 64 chips, allowing the 
simultaneous estimation of U = 3 users. Due to the oversampling factor N c = 4 we 
need to perform an FFT over 192 x 4 = 256 x 3 (real) samples, which can be efficiently 
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implemented as a mixed radix FFT, (i.e., one radix 4 and one radix 3). We considered 
a channel with MIP given by the set of delays 

r = [0, 0.9, 2.8, 4.7] 

normalized with respect to the chip interval, and by the relative path strengths 

a 2 = [0, -2, -7, -8.5] dB 

(the sum of the path gains is normalized to 1). The output SINR is not useful here 
since this would be dominated by the ISI introduced by the multipath channel. Then, 
as performance measure we use the Matched Filter Bound (MFB) [16], that is, the 
SINR at the output of the matched filter for the transmission of an isolated symbol 
(i.e., without ISI). This can be regarded as an upper bound to the performance of the 
optimal ML sequence estimator. 

Figures 5.13, 5.14, 5.15, and 5.16 show the MFB resulting from the estimated 
matched filter versus the MFB of an ideally matched filter, equal to the nominal £/7 0 , 
the total channel gain being normalized to one. Since in this case we are dealing 
with a CDMA system £ corresponds to the average chip energy. In the light of the 
previous results we consider the transmitter to employ nearest neighbor interpolation 
while the receiver resampling system adopts one of the methods of approach 2. The 
curves denoted by "Ideal interp." refer to the performance when ideal interpolation 
is used at both transmitter and receiver. Notice that in this case this curve does not 
have a horizontal asymptote since ISI is not present. Here, the degradation is due 
to the channel estimation errors (dominating the low £/7 0 region) and to resampling 
distortion (dominating the high £/7 0 region). 

With linear interpolation and Ltx — 32 phases at the transmitter an operating 
range between 0 and 35 dB is achieved. This shows that the degradation of linear inter- 
polation with respect to cubic and cubic spline interpolators evidenced in figures 5.8- 
5.11 is mainly due to ISI, namely inter-chip interference for a CDMA system. The 
effect of the inter-chip interference is negligible when the processing gain is large. 
Because of the presence of the processing gain a CDMA system typically operates in 
the low £/7o region, therefore linear interpolation can be chosen in this case. In fact 
there is no need to insist on more complicated interpolation techniques to increase the 
performance in the high £/7q region (namely higher than 10-20 dB). On the contrary, 
more advanced interpolation might be needed for non-spread spectrum system, espe- 
cially employing modulation scheme with large spectral efficiency (e.g., 16, 64, 256 
QAM modulations) which require to operate in the high £/7 0 region to achieve the 
desired performance. 

5.6 Conclusions 

We presented simple and effective algorithms suited for a flexible implementation of 
transmit and receive front-ends based on programmable DSPs, for SDR applications. 
We compared various solutions in terms of complexity and performance. The resulting 
architectures for transmitter and receiver are both based on the concept of IF sampling 
and resampling. At the transmitter, we can generate directly the IF analog signal via 
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Figure 5.14: Channel estimation with multipath channel and linear interpolation at the 
receiver 
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Figure 5.15: Channel estimation with multipath channel and cubic interpolation at the 
receiver 




Figure 5.16: Channel estimation with multipath channel and cubic spline interpolation 
at the receiver 
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the D/A converter, without using mixers, filters and I-Q adders. At the receiver, the 
IF analog signal is downsampled at a much lower rate than the Nyquist rate in or- 
der to generate a pass-band digital signal with low carrier frequency, suited to further 
processing. Resampling at both transmitter and receiver is needed to support sev- 
eral different symbol (or chip) rates. At the transmitter resampling can be integrated 
with pulse-shaping and can be efficiently implemented by a polyphase filter-bank fol- 
lowed by nearest-neighbor interpolation. At the receiver more accurate interpolation 
techniques are needed while the pulse-shaping matched filter can be more efficiently 
integrated with the channel matched filter, provided by a channel estimator. Our simu- 
lation results show that cubic-spline interpolation with approach 2 provide very small 
SINR degradation over a wide range of nominal £//o- Finally, we tested our front-end 
algorithms with an actual channel estimation technique. In terms of MFB, results are 
very good also for linear interpolation. In the case of CDMA with sufficient processing 
gain, the inter-chip interference due to resampling distortion might have a negligible 
effect and linear interpolation in the receiver might be a good low-complexity solution. 
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Chapter 6 



Architecture of the Software Radio 
Platform 



We describe a software radio architecture developed for providing real-time wide- 
band radio communication capabilities in a form attractive for advanced 3G system 
research. It is currently being used to implement signaling methods and protocol sim- 
ilar, but not limited to, evolving 3G radio standards (e.g. UMTS, CDMA2000). An 
overview of the hardware system is provided along with example software implemen- 
tation on both high performance DSP systems and conventional microprocessors. 

6.1 Introduction and Motivation 

Motivated by the worldwide activity around third generation mobile communication 
systems, Eurecom and EPFL have launched a joint project whose objective is to de- 
sign and implement a real-time software radio communication platform to validate 
advanced mobile communication signal processing algorithms. The right to transmit 
has been granted for one 5 MHz UMTS channel in both France and Switzerland for 
experimental purposes. 

The platform is characterized by the following major features: 

• Flexibility, achievable by a software driven system. 

• Duplex communication. 

• Multiple antennas transmit and receive signal processing (i.e. joint spatio tem- 
poral signal processing). 
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Flexibility remains a key word for a software-defined system. In our case it serves sev- 
eral purposes. For instance one may perform propagation channel measurements and 
transmitter characterization, evaluate the performance of different signal processing 
algorithms for both single user and multi-user systems under different operating con- 
ditions. Duplex communication is also necessary to allow higher-layer protocol testing 
and services, and to analyze more complex system aspects, such as multiple-access, 
power control, and optimize down-link signal processing from up-link measurements. 
In a second phase, the platform will allow multiple-antenna signal processing, or more 
generally, spatio-temporal signal processing (also known as array processing) since 
this a very promising ensemble of techniques able to significantly increase the capac- 
ity of wireless communication systems. 

The main focus of this Chapter is on the description of the basic platform architec- 
ture and signal processing to implement the essential physical layer level procedures 
of the UMTS -TDD mode. 

We have previously learned that one of the keys to the flexibility is to move the 
digital processing toward the antenna. We then choose to extend the digital signal 
processing to the Intermediate Frequency, implementing the required algorithms on a 
general purpose programmable processor. In our case a multi-band Radio Frequency 
section takes care of translating the desired signal from a fixed IF carrier to the required 
RF carrier in the transmitter, and vice-versa in the receiver. 

We are then concerned with the efficient generation of an IF analog signal from a 
base-band digital signal (transmitter front-end) and with the reverse operation (receiver 
front-end). Once the end-to-end transmitter and receiver architectures is defined the 
next step consists of synchronizing the receiver, estimating the channels associated 
with the users, and eventually detecting the transmitted symbols. All these operations 
are performed in real-time on our platform. 

The adopted modulation scheme (BPSK and QPSK), the spread signal bandwidth 
(Wj/, « 5 MHz), the spreading gain (DS-CDMA with allowed spreading gain of L = 
1, 2, 4, 8, 16) and the frequency band (around 2.1 GHz) are the same as those defined 
in the UMTS/TDD specifications [15]. The TDD mode has been chosen because under 
certain circumstances it allows the exploitation of the channel reciprocity between up- 
link and down-link in duplex communication, and also reduces the DSP computational 
load for a given bandwidth. 

6.2 System architecture 

In a first phase, the development of the platform has been based on a single antenna 
architecture for the Mobile Terminal as well as for the Base Station. The retained archi- 
tecture permits to easily enhance the capabilities of the platform without redesigning 
the essential hardware and software parts. 

The BS and the MT are based on the same hardware platform and differ only in 
the software implementation. The different hardware parts of the platform are highly 
partitioned in order to enable the use as many standard cards and components as pos- 
sible. 

The hardware portion of the test-bed consists of 4 elements which are under soft- 
ware control, namely 
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Figure 6.1: System Architecture 



1. a PCI bus based reconfigurable Data Acquisition Card (DAQ) based on a Field 
Programmable Gate Array (FPGA) 

2. an RF front-end 

• a single stage up/down-conversion from/to a 70 MHz IF carrier with Time- 
Division-Duplex multiplexing 

• 1 high-speed 12-bit pass-band sampling A/D converter working at rate 
/ado 

• 1 high-speed 12-bit up-sampling D/A converter working at rate /dac- 

• 8 slow D/A converters for controlling various amplifier gains on the RF 
card 

• control for various switches on RF card 

3. a clock card for generating sampling clocks (fixed frequency) and local oscilla- 
tors (programmable frequency) 

We show a simplified block diagram of the entire system for a single antenna in Fig- 
ure 6.1. 

We have considered two software implementations, the first using a combination 
of commercially available embedded DSP cards and a common PC and the second 
using the DSP units (e.g. MMX) of a standard PC under an operating system proving 
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hard real-time support (e.g. RTLinux). Both run in real-time and are compatible over 
the air. 

6.2.1 RF front-end 

A simplified block diagram of our RF front-end is shown in Figure 6.2. It was designed 
in conjunction with STMicroelectronics in Geneva, Switzerland. 

On reception (Rx), the RF signal is filtered, amplified and down-converted to a 70 
MHz IF. The Local Oscillator (LO) frequency is digitally tunable in steps of 500 Hz. 
The IF signal is amplified by a digitally tunable gain control and directly sampled at 
rate /adc- The samples are transferred via a ribbon cable using high-speed line drivers 
to the Data Acquisition Card. 

On the transmission end (Tx), samples from the DAQ drive a hardware up-sampling 
circuit and a high-speed D/A converter, working at rate /dac> to directly synthesize 
the 70 MHz IF signal. The IF signal is then amplified, up-converted to RF, filtered and 
amplified by a variable-gain power amplifier. 

The mode of the transceiver (i.e. transmission, reception, calibration) is fully con- 
trollable from the software portion of the platform. 

Special low-speed lines control the gains of the Rx variable-gain amplifier and Tx 
Power Amplifier, LO frequency, as well as the antenna switch. Although not included 
in Figure 6.2, automatic wide-band calibration capabilities (for both Tx and Rx) are 
included in the design, which are required for multi-antenna systems (see e.g. [57]). 
These will be used in a later stage of the project where the architecture will be extended 
to implement multiple antenna transceivers. The basic characteristics of the RF front- 
end are summarized in Table 6.1. 

Pass-band A/D conversion 

The main difference of the current implementation with respect to the UMTS-TDD 
3GPP standard is that we choose a chip rate f c = 3.6864 Mchips/s (CDMA2000 chip 
rate) and not 3.84 Mchips/s. Moreover, in order to avoid signal re-sampling 1 we 
suppose the rate /adc to be a multiple integer of the chip rate (i.e. /adc = N c f c ). In 
our implementation we set N c = 4 that yields 

/adc = 4/c = 14.7456 Msamples/s 

So, once the RF signal at the antenna has been down-converted to IF, it is sampled 
by the A/D converter at rate /adc ( see Figure 6.4). Because of the periodicity of the 
discrete-time signal spectrum, the resulting real sampled signal is given by 2 

y[n] =y w (n// AD c) 

lr The re- sampling techniques studied in Chapter 5 will be implemented in the next version of the 
software 

2 Notice that, by using the notation of Chapter 5 (see for example Figures 5.2 and 5.3), and avoiding 
interpolation at the receiver, we have y[n] = r[n] = yw(n/ /adc) 
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Figure 6.2: RF Front-end 



The signal y[n] is pass-band with a spectrum replica centered at /adc/4 as shown 
in Figure 6.4. In practice we applied (5.4) using I = 5 and choosing the sign (— ) in 
the denominator. The sampling frequency /adc satisfies the condition 



/adc>2WV wlOMHz 

required by the theory in Section 5.2.2. 

Although a base-band version of the received signal can be obtained by multiplying 
y[n] by (— j) n and then by low pass filtering, we will show further how the channel 
estimation and the data detection processes can be performed at pass-band. 



Pass-band D/A conversion 

At the transmitter side we send to the D/A converter the pass-band signal x f [n], ob- 
tained by the base-band signal x[n] using (5.5). The signal x f [n] shows a replica of the 
spectrum in / IF and is sampled at rate 

fs = N c f c = 14.7456 Msamples/s 

where we set N c = 4 as for the receiver. Also in this case the sampling rate satisfies 
the inequality f s > 2W^. 

In Section 5.2.3 it is shown that, in order to reduce the attenuation introduced by 
the D/A impulse response up-sampling is required. In our implementation we chose 

/dac = Sf 8 = 117.9648 Msamples/s 

so that the spectrum replica of x'[n] at fjp falls inside the main lobe of the D/A fre- 
quency response (see Figure 5.1(b)). 



6.2.2 Data Acquisition Card (DAQ) 

The data acquisition card is a PCI bus-mastering device permitting high-speed full du- 
plex parallel transfer of digital data from an external device. It contains the necessary 
glue logic which connects the input (A/D) and output (D/A) sample streams as well 
as some control signals to the main CPU/DSP. A PCI architecture was adopted since 
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Figure 6.3: Data Acquisition Unit 



it is the most general purpose bus architecture and is used on most standard PCs, as 
well as DSP systems. It consists of two components, namely a powerful Xilinx Field 
Programmable Gate Array (FPGA) XCV300 [58] and a bus-mastering PCI controller 
PLX9080 [59]. The format is a PCI Mezzanine Card (PMC) to allow for integration 
into both embedded DSP architectures and ordinary PCs. A simplified overview of the 
DAQ is shown in Figure 6.3. 

Firmware on the FPGA for formatting and transferring data via the DMA engines 
of the PLX9080 to the host (DSP, CPU) memory has been developed in VHDL (Very 
High Speed Integrated Circuit Hardware Description Language). After an initial con- 
figuration phase, transfers are continuous and completely transparent to the host, who 
just "sees" a circular buffer containing samples acquired from or to be transferred to 
the external RF front-end. 

The basic components of the DAQ are the following 

1. Line Drivers/Receivers for transfer from external devices via ribbon-cables 

2. A reconfigurable FPGA-based 16-bit bi-directional interface external device (up 
to 30 Msamples/s full-duplex) 

3. A PCI bus-mastering controller for direct transfer of samples to/from memory 
(DMA) on PCI-based signal processing units (e.g. DSP cards, high-performance 
PCs, workstations (SPARC, PowerPC, Alpha), embedded processor cards, etc.) 

4. A Processor Mezzanine Card (PMC) form-factor for maximum flexibility 
The FPGA is programmed via the PCI bus by serial download. 

6.2.3 Texas Instruments TMSC6201 implementation 

The embedded DSP architecture is based on a commercially available dual-DSP card 
(Spectrum Signal Processing Daytona [60]). The basic architecture is shown in Fig- 
ure 6.7 and is centered around 2 Texas Instruments TMSC6201 fixed-point DSPs. 
These DSPs are capable of providing a maximum of 1600 MIPS each. Our DAQ 
is placed on the local PCI bus of the DSP board and transfers samples to/from both of 
the memory buffers on the DSPs busses. These buffers are used as temporary storage 
as the internal (fast) memory of the DSP is rather small. The DSP DMA engines take 
care of automatically transferring data to their internal memory concurrently with the 
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Figure 6.4: Receiver front-end 

signal processing functions. One DSP is used exclusively for transmission functions 
and the other for reception. 

The DSPs are used for the front-end processing as described in the previous sec- 
tions. Symbol rate data is transferred via the PCI bus to the PC which hosts the DSP 
card. This data is processed by the Pentium and handles tasks such as Viterbi decoding, 
carrier frequency offset compensation, higher layer protocol stacks, etc.. 



6.2.4 RTLinux-based PC implementation 

The second implementation does not rely on embedded DSPs. Here the DAQ is placed 
on the master PCI bus of a workstation, possibly multiprocessor, running the hard real- 
time extension to the Linux operating system, RTLinux[61]. The software radio runs 
in kernel space and is integrated into the IPv4 (Internet Protocol) subsystem of Linux 
as a network device. 

This x86 implementation makes use of the MMX (multimedia extensions) SIMD 
(single instruction multiple data) instructions for obtaining maximum processor effi- 
ciency for intensive DSP computations. All DSP routines use fixed-point arithmetic 
and are written in C with embedded assembly macros for time-critical code sections. 
We typically make use of 

1. MMX packed 16-bit arithmetic (multiply, add, MAC, interleaving, etc.) 

2. loop unrolling 

3. software pipelining 

Because of the high-level software structure, this should be portable to other processor 
architectures (e.g. PowerPC, Alpha, etc.). For the same reason, it is easily portable 
to large-scale SMP (symmetric multi-processing) platforms which could be useful for 
advanced base-station implementations. 
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6.3 Physical layer digital signal processing 

This section gives an overview of some theoretical principles on which the platform 
software has been implemented. We have implemented variant of the UMTS-TDD 
3GPP standard and are now implementing a complete subset of layers 1 and 2 of the 
true standard (including 1.28 Mchips/s version). After D/A conversion, a pass-band 
filter centered at /if removes the other replicas, generating the desired IF modulated 
signal. The discrete-time modulation by /dac/4 requires a negligible computational 
cost since it corresponds to change alternatively the signs of x[n] as can be noticed by 
(5.3). 



6.3.1 Receiver 

In this section we give a description of the receiver front-end architecture. Then we 
focus on channel estimation, matched filter synthesis and symbol detection. 



Frame synchronization 

In order to minimize the computational cost, the frame synchronization is achieved 
using a simplified version of the "MAX" timing estimator defined in Section 3.5. In our 
implementation we correlate the input signal y[n] with a filter matched to the primary 
synchronization sequence s tr but not to the chip pulse ift(t). The timing estimation is 
achieved by filtering the pass-band received signal y[n] as 



n} = 



255 



k=0 



255 



k=0 



(6.1) 



and averaging the y ps [^] 2 over several frames. This is typically the most computation- 
ally intensive part of the receiver front-end since it requires a fairly long filter operating 
at the sampling rate. The 3GPP standard uses a hierarchical structure for the primary 
synchronization sequence (see Section 2.1.3) which allows the filter to be implemented 
as a concatenation of 2 FIR filters of length 16. The filter output y ps [n] can then be 
obtained by 



2/psW = 



15 



£W*]*[n + 64fc] 
k=o 

15 



15 



k=0 



^8 psi [Jfc]5[n + 4fc] 

k=0 
15 



*y[n] 



k=0 



(6.2) 



reducing the computational cost by a factor 8. The maximum output of this filter is 
used to adjust the receive signal strength (via a variable gain IF amplifier) and syn- 
chronization is achieved when the maximum is greater than a pre-defined threshold. 
Note that this filtering operation involves purely real quantities. A block diagram of 
the primary synchronization signal processing is shown in Figure 6.5. 



6.3 Physical layer digital signal processing 



93 



y[n\. 







s ps 2 




(■? 







>■ 





^Iframe 



l-a 
V 



^0 




Figure 6.5: Scheme of the implemented primary synchronization signal processing 
Channel estimation 

The training-sequence based multiuser channel estimation procedure for block syn- 
chronous CDMA described in the UMTS/TDD standard has been implemented fol- 
lowing the FFT approach given in Section 5.5. The algorithm provides the estimates g 
of the channel g. A block scheme of the channel estimation procedure is summarized 
in Figure 6.6. 



Matched filter synthesis and data detection 

Given the channel estimates, we are then interested in synthesizing the Matched Fil- 
ters f u , for u = 1, . . . , U, matched to the cascade formed by the user data spreading 
sequence, the chip pulse shape filter and the user channel. This is also shown in Fig- 
ure 6.6. The overall impulse response for the u ih user is given by 

L u -l 

fu(t)= ]T s u [k]g u (t-k/f c ) (6.3) 
k=o 

where L u is the 11 th user spreading gain. Using the sampled channel estimate g we 
synthesize the discrete-time filters 

fu[k}= fu(t)\t=k/f An c 

matched to the overall response as follows 

fu[k] = s u [k]*g u [k] (6.4) 
In order to extract the data symbols we filter the received signal with the MF obtaining 

«»[*] = y[*]*/i[-*] (6-5) 

In a pass-band approach y[k] is real while f u [k] is complex so the product requires two 
real multiplications. In the base-band case, instead, the samples of y[k] would be at 
half the sampling rate but would be complex. So the two complexities are identical. 

Notice that the signal after matched filtering is still pass-band and the symbol esti- 
mates after sub-sampling are given by 



a u [k]=v' u [NN c k] 



(6.6) 
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Data field 1 Midamble Data field 2 Guard 




Figure 6.6: Channel Estimation and Matched Filter Synthesis 



where v f u [k] is the complex envelope of the MF output. But since 

= H) k v u [k] (6.7) 

substituting into (6.6) and for N c = 4 we get 

&u[k] = H) LuNck v u [L u N c k] = v u [±L u k] (6.8) 

Then the symbol estimates are given by sub-sampling the MF output at symbol rate 
without taking care of the demodulation. 

Carrier synchronization and decoding 

The carrier synchronization is done at symbol rate with a classical decision directed 
algorithm [16]. The algorithm then takes a decision on the symbols and recovers the 
data. 

6.4 Validation of the existing platform 

The described platform has been validated by the transmission and the reception of 
two user full-duplex real-time video flows in an indoor environment. Two H263 video 
streams are transmitted in parallel and decoded in real time. For this we use the fol- 
lowing parameters: 

• spreading factor 16 

• bit rate 397 kbit/s (peak) per user 

• TDD configuration: 1 Tx slot followed by 1 Rx slot (transmission is done every 
2 slots). 

• two synchronous full-duplex streams per slot 



• RF band: 5 MHz at 2. 1 GHz 



6.5 Conclusion 
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Figure 6.7: Embebbed DSP Architecture 

6.5 Conclusion 

This first demonstration shows that the architecture of the platform is capable of sus- 
taining real-time communications and is thus promising for future developments. The 
platform is currently being enhanced and opened to both industrial and academic col- 
laboration. 

The enhancements will consist of 

• support for multiple antenna transceivers 

• more sophisticated signal processing algorithms 

• multi user detection 

• layer 2 (RLC, MAC) functionality 

Collaboration has already begun under the label of the RNRT (Reseau National de la 
Recherche en Telecommunications) financing program organized by the French Min- 
istry of Industry and Finance. Three projects have been initiated covering the following 
topics: 
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Frequency Band 


2100-2170 MHz 


Bandwidth 


5 MHz (initially) 


Transmit Power (per antenna) 


1 W 


Receiver Sensitivity 


-100 dBm 


Noise Figure 


< 5dB 
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> -7 dBm 


Duplex mode 


Time Division 


Rx Gain Control 


Digital tuning, ldB steps over 40 dB 


Tx Power Control 


Digital tuning, ldB steps over 80 dB range 


Local Oscillator 


Digital tuning, step: few kHz in each band 


RF Calibration 


Digital control, Tx and Rx 


Direct IF sampling (70 MHz) on Rx 


12-bit A/D @ 14.7456 Msamples/s 


Direct IF synthesis (70 MHz) on Tx 


12-bit D/A @ 117.9648 Msamples/s 


Digital Interface 


High speed: LVDS 

Low speed: 3.3V CMOS line drivers 



Table 6.1: RF Front-end characteristics 



• radio subsystem improvement (flexibility and sensitivity) 

• compliance with the 3GPP UMTS/TDD specification 

• higher-level protocol stacks 

• integration to an IPv6 experimental backbone 



Chapter 7 



General conclusions 



In this thesis we addressed the problem of improving the physical layer reliability of 
W-CDMA third generation transceivers by designing simple and low complexity al- 
gorithms suited to be implemented on a real system. The work was developed in a 
framework of a joint project between Eurecom and EPFL whose objective is to design 
a real-time Software Defined Radio platform able to validate advanced signal process- 
ing algorithms. The UMTS-TDD 3GPP standard, on which the platform has been 
initially targeted, was used throughout this work as a reference system for compar- 
isons. These algorithms, suited to SDR, cover some major issues of the physical layer 
signal processing such as synchronization, and multiuser detection. Since the aim was 
to propose solutions for a practical implementation, particular care has been devoted 
to the algorithm complexity evaluation. Finally, the proposed algorithms have been 
tested in a realistic scenario, including frequency selective channels, fast and slow 
fading conditions, power control, and timing estimation errors. 

The initial synchronization between a Mobile Terminal and a Base Station has 
been covered by Chapter 3 where we derived a Joint Maximum Likelihood slot timing 
estimator for bursty pilot signals in the presence of multipath fading channels, with the 
following conclusions: 

• While the traditional "MAX" slot timing estimator provides good performance 
for short channels it performs very poorly in the presence of sparse channels 

• The proposed "JML" estimator shows good performance in all the considered 
conditions, and is very robust to mismatch. Moreover, at the same time, it pro- 
vides an estimate of the multipath intensity profile and of the interference plus 
noise power spectral density. 

• With respect to a PN sequence, the primary synchronization sequence proposed 
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Chapter 7 General conclusions 



by the UMTS-TDD standard reduces the performance of the slot timing estima- 
tors because of its less peaky aperiodic autocorrelation function. 

Future research in this area may consider the exploitation of the additional informa- 
tions provided by the JML algorithm for improving the communication reliability im- 
mediately after the primary synchronization process (e.g. demodulation of the UMTS- 
TDD secondary synchronization signal and detection of the scrambling sequence). 

In order to increase the system capacity, in Chapter 4 we designed an iterative turbo 
equalizer and multiuser decoder/detector for TD-CDMA systems. Simulation results 
showed that: 

• Simple and low-complexity multiuser detection techniques outperforms the pow- 
erful and high complexity channel coding envisaged by the UMTS-TDD norm. 

• High channel loads can be obtained by using simple convolutional codes and 
Viterbi decoding instead of using the more computationally intensive SISO de- 
coding. 

• Since the proposed algorithm is very simple and it is robust to severe ISI condi- 
tions it is well suited to be implemented on a UMTS-TDD base station. 

The problem of designing a SDR architecture able to support different standards 
with the same hardware has been investigated in Chapter 5 where we proposed some 
resampling techniques for the transceiver front-end. These techniques have been tested 
in realistic scenarios, evaluating the SINR degradation with respect to the ideal inter- 
polation. From the results of this Chapter the following conclusions can be pointed 
out: 

• Resampling techniques are an efficient way to make the baseband signal pro- 
cessing independent on the Intermediate Frequency. 

• Resampling techniques at both transmitter and receiver have a complexity af- 
fordable by the technology available on the market. 

• In general, most part of the computational complexity required by the interpola- 
tion techniques stays in the receiver. 

• With a negligible computational cost the transmitter front-end can directly gen- 
erate the IF signal at the output of the D/A converter. 

• Even if interpolation distorts the signal, a careful choice of the interpolation 
technique allows a wide SINR operating range at the receiver. 

• In the case of CDMA with sufficient processing gain, the inter-chip interference 
due to resampling distortion might have a negligible effect. 

Moreover we have shown that with a common PC, equipped with conventional mi- 
croprocessor and/or with high performance DSPs it is possible to implement real-time 
UMTS-TDD Base Stations and Mobile Terminals, compliant to the 3GPP standard. 
The theoretical research developed throughout this thesis has been of main relevance 
for the platform physical layer implementation and, in some cases, addressed solutions 
to improve the UMTS-TDD standard. 
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